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Enclosed  is  a  copy  of  report  FAA/RD-90/25  entitled  "Rotorwash 
Computer  Model  -  User's  Guide".  This  model  is  based  on 
measured  data  and  it  provides  a  means  for  heliport  and 
airport  designers  to  make  decisions  on  some  of  the  more 
complex  rotorwash-related  design  issues. 


The  majority  of  helicopters  are  light  in  weight  and  it  is 
rare  for  them  to  cause  a  rotorwash-related  accident  or 
incident.  However,  with  heavy  rotorcraft,  rotorwash-related 
mishaps  are  more  of  a  concern.  The  current  FAA  Heliport 
Design  Advisory  Circular  addresses  this  issue  to  a  limited 
degree.  By  and  large,  however,  the  assurance  of  safety  is 
the  responsibility  of  the  pilot. 


Anticipating  the  introduction  of  large  tiltrotor  and  the 
expanded  use  of  larger  helicopters,  it  is  appropriate  to 
consider  whether  protection  against  rotorwash  mishaps  should 
continue  to  depend  so  heavily  on  pilot  judgment.  Would  it 
not  be  better  to  provide  an  additional  safety  margin  by 
addressing  this  issue  via  facility  design  and  operational 
procedures  guidelines?  This  is  the  avenue  that  the  FAA  has 
been  pursuing. 


The  FAA  approach  to  this  task  is  four-fold: 

a.  to  measure  the  rotorwash  of  helicopter  and  tiltrotor 
aircraft. 

b.  to  develop  and  validate  a  rotorwash  computer  model  based 
on  the  available  rotorwash  data. 

c.  to  analyze  rotorwash  induced  mishaps  to  see  if  we  can 
determine  the  thresholds  at  which  rotorwash  becomes  a 
potential  hazard. 

d.  to  apply  the  model  to  an  analysis  of  a  variety  of 
operational  scenarios  using  these  thresholds  and  to 
determine  ways  of  alleviating  rotorwash-related  accidents 
and  incidents. 


This  report  is  the  fifth  in  a  set  of  documents  on  the  issue 
of  rotorwash.  The  subject  is  technically  complex  and 
additional  research  and  development  is  required.  We  want  to 
work  with  industry  in  developing  cost-effective  solutions. 
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Thus,  we  welcome  any  suggestions  that  you  may  wish  to  make  in 
this  regard.  Please  reference  this  report  and  send  your 
suggestions  to: 

Vertical  Flight  Program  Office,  ARD-30 
Federal  Aviation  Administration 
800  Independence  Ave.,  SW 
Washington,  DC  20591 
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SUMMARY 


The  user's  guide  for  the  Rotorwash  (ROTWASH)  Analysis 
program  is  documented  in  this  report.  This  documentation 
provides  step-by-step  descriptions  on  the  use  of  each 
analysis  option  and  a  listing  of  the  IBM  PC/PC-compatible 
based  FORTRAN  77  software.  A  brief  introductory  section  to 
the  user's  guide  describes  the  history  of  the  ROTWASH 
analysis  software.  References  for  the  mathematical  modeling 
utilized  in  the  various  analysis  modules  are  also  included 
in  the  report. 
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1989  and  June  1990.  The  work  was  performed  by  Mr.  Samuel  W. 
Ferguson  at  EMA. 


TABLE  OF  CONTENTS 


Page 

1.0  Introduction  .  1 

2.0  ROTWASH  Program  User's  Guide  .  3 

2.1  Getting  Started  .  3 

2.2  Input  Data  Requirements  .  4 

2.3  Analysis  Problem  Definition  .  6 

2.4  The  Wall  Jet  Option  .  8 

2.5  The  Twin-Rotor  Interaction  Plane  Option  .  14 

2.6  Disk  Edge  Vortex  Analysis  Option  .  22 

2.7  Personnel  Overturning  Force  and  Moment 

Analysis  .  24 

2.8  Particulate  Cloud  Analysis  Option  .  29 

3.0  ROTWASH  Program  Data  Output  File  Formats  .  35 

References  .  41 

Appendix  A  Letter  Report,  Oblique  Impingement  of 

Rotor  Downwash  with  the  Ground  Plane  .  A-l 

Appendix  B  Rotorcraft  Design  Data  in  ROTWASH  Format  . .  B-l 

Appendix  C  ROTWASH  Program  FORTRAN  77  Listings  .  C-l 


i  i  i 


LIST  OF  FIGURES 


Page 

Figure  1  ROTWASH  Program  Header  Output  .  3 

Figure  2  ROTWASH  Program  Input/Output  Control  Menu  ....  4 

Figure  3  Master  Input  Data  Menu  .  5 

Figure  4  ROTWASH  Program  Logic/Comment  Menu  .  7 

Figure  5  Velocity  Analysis  Option  Menu  .  8 

Figure  6  Hazard  Analysis  Option  Menu  .  9 

Figure  7  Rotorwash  Flowfield  Characteristics  of  Single 
and  Twin-Rotor  Configurations  Operating  in 
Close  Proximity  to  the  Ground  .  9 

Figure  8  Vertical  Section  of  the  Nondimensional  Wall 

Jet  Velocity  Profile  .  10 

Figure  9  Velocity  Profile  Status  Menu  .  11 

Figure  10  XV-15  Flight  Test  Data  Measurement  Locations  .  12 

Figure  11  Wall  Jet  Velocity  Profile  Output  Format  .  13 

Figure  12  XV-15  Velocity  Profile  Correlation  .  13 

Figure  13  Interaction  Plane  Velocity  Profile  Output 

Format  .  15 

Figure  14  XV-15  Interaction  Plane  Velocity  Profile 

Correlation  .  16 

Figure  15  Diagram  of  the  Ground  Vortex  Structure  .  17 

Figure  16  CH-53E  Input  Data  for  the  Ground  Vortex 

Example  . 18 

Figure  17  Ground/Disk  Vortex  Input  Data  Menu  .  18 

Figure  18  Horse  Shoe  Vortex  Geometry  for  the  Ground 

Vortex  Analysis  Option  .  20 

Figure  19  Ground  Vortex  Analysis  Intermediate  Output  ...  20 

Figure  20  Ground  Vortex  Circulation  Strength  as  a 

Function  of  Advance  Ratio  .  21 

Figure  21  Ground  Vortex  Velocity  Field  Output  Data  .  21 


IV 


Page 


Figure  22  Disk  Edge  Vortex  Geometry  Definition  .  22 

Figure  23  Disk  Edge  Vortex  Option  Intermediate  Output  . .  23 

Figure  24  Velocity  Field  Output  Data  from  the  Disk 

Edge  Vortex  Option  .  24 

Figure  25  Overturning  Force  and  Moment  Calculation 

Procedures  .  25 

Figure  26  Overturning  Force  and  Moment  Data  Menu  .  26 

Figure  27  Simplified  Output  Format  for  the  Overturning 

Force  and  Moment  Analysis  .  27 

Figure  28  Bell  XV-15  Overturning  Force  as  a  Function 
of  Distance  Along  the  Interaction  Plane 
(0  Degree  Azimuth)  .  28 

Figure  29  Menu  Specification  of  Detailed  Force/Moment 

Output  .  29 

Figure  30  Detailed  Output  Format  for  the  "Large" 

Person  Analysis  Option  .  30 

Figure  31  Detailed  Output  Format  for  the  Second  Screen 

Frame  of  the  "Small"  Person  Analysis  Option  . .  31 

Figure  32  Particulate  Cloud  Analysis  Geometry 

Definition  .  31 

Figure  33  Approximate  Values  for  the  Terrain  Erosion 

Factor  as  a  Function  of  Terrain  Type  .  33 

Figure  34  Particulate  Cloud  Analysis  Option  Output  .  34 

Figure  35  Example  Wall  Jet  Option  Graphics  File  Format  .  37 

Figure  36  Example  Interaction  Plane  Option  Graphics 

File  Format  .  38 

Figure  37  Example  Personnel  Overturning  Force  and 

Moment  Graphics  File  Format  Created  with  the 
Wall  Jet  Analysis  Option  .  39 

Figure  38  Example  Personnel  Overturning  Force  and 

Moment  Graphics  File  Format  Created  with  the 
Interaction  Plane  Analysis  Option  .  40 


v 


LIST  OF  TABLES 


Table  1 


Page 


Graphics  File  /  Source  Code  Reference 

Matrix  .  36 


vi 


LIST  OF  ACRONYMS 


AGL 

DAIP 

DFAC 

DFRC 

DOT 

DOS 

FAA 

GW 

HAGL 

HROTOR 

Hz 

NASA 

PSF 

ROTHAZ 

ROTWASH 

WAGL 


Above  ground  level 

Distance  along  the  interaction  plane,  feet 

Distance  from  aircraft  center,  feet 

Distance  from  the  rotor  center,  feet 

Department  of  Transportation 

Disk  Operating  System 

Federal  Aviation  Administration 

Gioss  weight,  pounds 

Rotor  or  sensor  height  above  ground  level,  feet 
Rotor  height  above  ground  level,  feet 
Hertz,  cycles  per  second 

National  Aeronautics  and  Space  Administration 

Pounds  per  square  foot 

Rotorwash  and  Hazard  Analysis  Program 

Rotorwash  Analysis  Program 

Wheel  height  above  ground  level,  feet 


Vll 


1.0  INTRODUCTION 


Rotorwash  flowfield  characteristics  of  helicopter, 
tiltrotor,  and  tiltwing  aircraft  are  quite  complex  in  close 
proximity  to  the  ground.  All  theoretically  rigorous 
analyses  of  these  characteristics  require  quite  complex 
computer-based  mathematical  models  and  large  high-speed 
computers  for  solution.  Unfortunately,  in  most  situations, 
these  types  of  resources  are  either  unavailable, 
unaffordable,  or  can  not  effectively  be  used  when 
information  is  required  in  a  short  frame  of  time. 

Therefore,  in  response  to  the  need  for  a  simple  yet 
effective  analysis  tool,  the  U.S.  Department  of 
Transportation  (DOT)  sponsored  the  development  of  the 
Rotorwash  Hazard  (ROTHAZ)  Analysis  methodology  and  its 
associated  FORTRAN  77  software.  This  research  effort  was 
completed  in  1986  and  is  documented  in  reference  1.  In 
1989,  the  Federal  Aviation  Administration  identified  the 
need  for  analysis  of  the  rotorwash  characteristics  of  eleven 
tiltrotor  and  tiltwing  aircraft  configurations.  The  results 
of  this  study  are  documented  in  DOT/FAA/RD-90/ 16 , 

"Evaluation  of  Rotorwash  Characteristics  for  Tiltrotor  and 
Tiltwing  Aircraft  in  Hovering  Flight"  (reference  2) .  As  a 
result  of  this  research  effort,  the  ROTHAZ  analysis 
methodology  has  been  improved  and  is  now  referred  to  as  the 
Rotorwash  (ROTWASH)  analysis  methodology. 

Technical  features  of  -he  mathematical  models  used  in 
the  ROTWASH  analysis  continue  to  be  based  on  the  ROTHAZ 
research  effort  with  minor  improvements.  One  of  these 
improvements  allows  for  the  calculation  of  oblique 
impingement  of  rotorwash  with  the  ground  plane.  Details  on 
the  mathematical  modeling  approach  used  for  this  effect  are 
contained  in  Appendix  A.  Each  of  the  analysis  options 
included  in  the  earlier  FORTRAN  77  software  is  still 
contained  in  the  new  software.  However,  extensive  changes 
have  been  developed  for  improvement  of  the  IBM  PC/PC- 
compatible-based  software  user  interface.  This  software 
interface  has  been  completely  restructured  using  a  much  more 
flexible  and  timesaving  menu  structure.  New  options  have 
also  been  developed  for  the  output  of  analysis  results  in 
several  different  formats.  These  options  include  printout 
of  data  to  the  screen,  printer,  and  disk.  Data  files 
printed  out  to  disk  can  also  be  formatted  for  acceptance  by 
PC-based  graphics  programs  for  the  first  time. 

The  purpose  of  this  report  is  to  document  each  feature 
of  the  ROTWASH  analysis  software  in  a  written  user's  guide, 
and  this  documentation  is  provided  beginning  in  Section  2.0. 
However,  it  must  be  stated  that  this  program  and  its 
documentation  are  not  targeted  for  the  novice  user.  Since 
rotorwash  flowfields  are  quite  complex  and  new  data  are  in 
the  process  of  being  acquired,  the  ROTWASH  analysis 


methodology  must  be  considered  an  evolving  engineering  tool 
at  the  present  time.  Therefore,  it  is  important  that  the 
user  be  familiar  with  this  field  of  research  to  effectively 
interpret  and  use  ROTWASH  program  output.  All  input  data 
requirements,  technical  analysis  options,  and  data  output 
options  are  described  in  detail  in  the  user's  guide  and 
example  printouts  are  provided  for  user  reference. 

Commentary  is  also  provided  that  identifies  known  strengths 
and  weaknesses  of  the  various  program  technical  modules. 
Sources  of  flight  test  data  for  correlation  with  output  from 
the  technical  modules  are  referenced  for  researchers 
interested  in  modifying  the  mathematical  models.  Data 
describing  most  modern  types  of  rotorcraft  in  the  ROTWASH 
input  data  format  are  summarized  in  Appendix  B.  Appendix  C 
provides  a  listing  of  the  ROTWASH  program  as  compiled  and 
linked  for  use  on  IBM  PC/PC-compatible  computers  using 
MICROSOFT  FORTRAN,  Version  5.0. 
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2.0  ROTWASH  PROGRAM  USER'S  GUIDE 


The  style  of  the  ROTWASH  user's  guide  presented  in  this 
section  is  primarily  narrative.  This  format  is  designed  to 
guide  the  reader  through  a  step-by-step  explanation  on  the 
use  of  each  program  software  option.  Example  output  from 
each  option  is  presented  exactly  as  it  would  be  viewed  by 
the  user  on  a  video  terminal  for  reference  purposes. 

Examples  of  user  keyboard  input  are  presented  in  <BOLD>  text 
as  an  aid  to  the  reader. 

2.1.  GETTING  STARTED 

The  ROTWASH  program  is  executed  by  typing  the  program 
name  at  the  DOS  system  prompt: 


<ROTWASH> 


To  avoid  possible  system  errors,  the  user  should  execute  the 
program  from  the  directory  containing  the  program  (or  set 
the  appropriate  DOS  system  PATH  command) .  The  user  should 
also  be  aware  that  menu  and  data  printouts  to  the  screen 
will  not  work  correctly  unless  the  device=C:\DOS\ANSI.SYS 
command  is  contained  in  the  CONFIG.SYS  file. 

The  ROTWASH  program  responds  with  the  screen  output  as 
presented  in  figure  1.  This  output  is  the  ROTWASH  program 
header  page. 


ROTWASH  PROGRAM 


ROTORCRAFT  DOWNWASH  HAZARD  ANALYSIS 

EMA/COMPUTATIONAL  METHODOLOGY  ASSOCIATES 
***  PROGRAM  VERSION  2.0,  APRIL  1991  *** 


PRESS  <RETURN> 


FIGURE  1  ROTWASH  PROGRAM  HEADER  OUTPUT 


After  typing  the  carriage  return  <RETURN>  or  <ENTER> 
key,  the  user  is  asked  to  specify  the  path  for  ROTWASH 
program  input  and  output  (I/O).  The  only  option  for  input 
at  the  present  time  is  the  terminal  keyboard.  This  is 
specified  by  typing  <CON>  for  console  or  keyboard  (lowercase 
inputs  such  as  <con>  are  also  permitted) .  Program  output 
may  be  sent  to  one  of  three  different  locations.  These 
output  locations  are: 


OUTPUT  OPTION 


TYPING  COMMAND 


Screen 

Printer 

Disk  Plotting  File 


<CON> 

<PRN> 

<PLT> 


An  example  of  the  video  screen  output  for  this  program  menu 
is  presented  in  figure  2  where  the  <CON>  option  has  been 
chosen  for  both  input  and  output. 


I/O  CAN  BE  DIRECTED  TO  FILES  OR  DEVICES 
VALID  DEVICES  ARE  AS  FOLLOWS: 


<CON> 

<PRN> 

<PLT> 


==>  CONSOLE 
==>  PRINTER 
==>  GRAPHICS  FILE 


ENTER  INPUT  FILE/ DEV  NAME  ==>  CON 
ENTER  OUTPUT  FILE/DEV  NAME  ==>  CON 


FIGURE  2  ROTWASH  PROGRAM  INPUT/OUTPUT  CONTROL  MENU 


2.2  INPUT  DATA  REQUIREMENTS 

Rotorcraft  characteristics  and  atmospheric  conditions 
that  are  common  to  all  program  options  are  input  to  the 
program  using  the  master  input  data  menu.  Four  basic 
configurations  of  rotor  or  propeller  driven  aircraft  can  be 
represented  using  this  menu.  These  configurations  include 
single  and  tandem  rotor  helicopters,  tiltrotors,  and  twin- 
propeller  tiltwings.  This  menu  is  presented  to  the  user  as 
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ROT WASH  USER  INPUT  DATA  MENU 


CODE 

PARAMETER 

VALUE 

UNITS 

A 

NUMBER  OF  ROTORS  (1  OR  2) 

2 

-ND- 

B 

HUB  TO  HUB  ROTOR  SEPARATION 

32.2 

FT 

C 

ROTOR  RADIUS 

12 . 5 

FT 

D 

GROSS  WEIGHT 

13000.0 

LB 

E 

FUSELAGE  DOWNLOAD  FACTOR 

13 . 0 

PCT 

F 

ROTOR  HEIGHT  ABOVE  GROUND 

37 . 0 

FT 

G 

SHAFT  TILT  ANGLE  (<20  DEG) 

.  0 

DEG 

H 

AIR  DENSITY  RATIO 

1.0000 

ND 

I 

AMBIENT  WIND  (-10  TO  10  KT) 

.0 

KT 

ENTER 

DATA  ENTRY  CODE  OR  <RETURN> 

TO 

CONTINUE  ==> 

D 

GROSS 

WEIGHT  =  13000. 

0 

ENTER 

NEW  VALUE  OR  <RETURN>  TO  LEAVE 

AS  IS  ==> 

13500. 

FIGURE  3  MASTER  INPUT  DATA  MENU 


shown  in  figure  3  after  the  I/O  menu  is  completed.  The 
default  values  provided  in  the  menu  define  the  Bell  XV-15 
tiltrotor.  Design  data  describing  most  other  modern  types 
of  rotorcraft  are  provided  in  Appendix  B. 

The  first  input  data  variable  listed  on  the  menu  is  the 
number  of  rotors  (or  propellers) .  This  value  will  always  be 
2  except  for  single  main  rotor  helicopters  (the  tail  rotor 
is  not  considered  a  lifting  rotor  by  the  ROTWASH  analysis) . 
The  next  four  menu  variables  define  the  design  details  of 
the  rotorcraft  configuration.  These  variables  are  the  hub- 
to-hub  separation  distance  for  twin-rotor  configurations  (in 
feet),  rotor  radius  (feet),  rotorcraft  gross  weight 
(pounds) ,  and  the  rotor  on  wing  download  factor  (percent) . 
The  download  factor  represents  the  percent  increase  in  hover 
rotor  thrust  required  to  overcome  the  rotor-induced  vertical 
drag  force  on  the  airframe  (where  thrust  is  initially 
assumed  to  equal  gross  weight) .  For  most  helicopters,  this 
value  is  less  than  5  percent.  Download  on  tiltrotors  can 
generally  be  expected  to  vary  from  8  to  13  percent  because 
of  the  large  wing  area  under  the  rotor. 

The  next  two  input  data  variables  define  the  position 
of  the  rotor  with  respect  to  the  ground  plane.  The  rotor 
height  (feet)  is  defined  as  the  distance  from  the  ground 
surface  to  the  plane  of  the  rotor.  The  mast  angle  (degrees) 
is  defined  as  the  tilt  of  the  rotor  plane  with  respect  to 
the  ground  plane  (this  improvement  to  the  ROTWASH 
methodology  is  described  in  detail  in  Appendix  A) .  The  mast 
angle  is  defined  as  0  degrees  when  the  plane  of  the  rotor  is 
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parallel  to  the  ground  plane  and  a  positive  angle  is  a 
forward  tilt  of  the  mast  (which  is  presently  limited  in  the 
program  to  20  degrees) .  Since  most  rotorcraft  hover  with 
the  plane  of  the  rotor  parallel  to  the  ground,  it  is 
recommended  that  caution  be  exercised  when  non-zero  values 
are  used  for  mast  angle.  Non-zero  uses  of  the  variable 
might  involve  hover  investigations  for  tiltwing  aircraft 
with  the  wing  tilted  forward  (where  a  pitch  fan  provides  the 
trim  pitching  moment).  Also,  limited  rotorwash 
investigations  for  takeoff  and  landing  maneuvers  might  be 
attempted  if  approximate  rotor  thrust  levels  are  known  for 
various  segments  of  the  maneuvers. 

The  last  two  input  data  variables  on  the  menu  define 
the  atmospheric  conditions.  The  air  density  ratio  is 
defined  as  the  ratio  of  the  desired  air  density  to  the  sea 
level  standard  air  density  (0.0023769  slugs/ feet 3 ) .  Ambient 
wind  speed  (knots)  is  specified  by  the  user  up  to  a  limit  of 
10  knots.  Values  greater  than  this  limit  are  believed  to 
invalidate  several  empirically  determined  mathematical 
modeling  assumptions  (limits  are  discussed  in  reference  1). 

The  mechanics  of  using  the  menu  are  quite  simple.  The 
user  types  in  the  code  value  for  the  variable  to  be  changed 
and  then  types  <RETURN> .  The  next  prompt  asks  for  the  new 
parameter  value.  After  this  value  and  another  <RETURN>  are 
typed,  the  menu  is  rewritten  to  the  screen  with  the  new 
value.  This  simple  process  is  continued  until  the  user 
specifies  each  variable  to  its  desired  value.  At  this 
point,  the  <RETURN>  key  is  typed  by  itself.  This  menu  can 
also  be  reached  from  most  of  the  other  menus  in  the  program 
whenever  the  user  decides  that  the  basic  configuration  needs 
to  be  modified.  This  is  accomplished  by  typing  <N>  for  NEW 
CASE  when  the  option  is  offered. 

2.3  ANALYSIS  PROBLEM  DEFINITION 

After  the  master  input  data  menu  is  completed,  the  user 
specifies  the  desired  type  of  analysis  option.  Figure  4 
presents  the  program  logic/comment  menu  and  the  associated 
list  of  default  values.  This  menu  has  two  groupings  of 
parameters  which  need  to  be  specified. 

The  first  parameter  on  the  menu  specif  ion  the  choice  of 
either  the  velocity  calculation  analysis  option  or  the 
hazard  analysis  option.  The  velocity  analysis  opt  ion  is  the 
default  option  on  the  menu.  This  option  is  otherwise  chosen 
by  typing  the  code  <A>,  then  <RETURN> ,  then  <V>,  and  finally 
<RETURN>.  The  same  process  is  used  to  spec  i  .'  y  the  hazard 
analysis  option  except  that  <H>  is  substitute  I  : or  <V>  as 
the  parameter  value .  The  second  menu  parameter  provides  an 
interactive  toggle  switch  for  the  option  which  writes  out 
data  files  to  disk  for  graphics  programs  ('I  his  parameter  is 
also  offered  as  an  option  on  the  initial  ROli.’ASH  menu  by 


ROTWASH  PROGRAM  LOGIC/ COMMENT  MENU 

CODE  PARAMETER  VALUE 

A  ANALYSIS  TYPE,  <V>  OR  <H>  V 

B  GRAPHICS  FILE,  <Y>  OR  <N>  N 

USER  INPUT  COMMENTS  (FOR  " PRN"  AND  "PLT"  OUTPUT) 

< -  50  SPACES  - > 

C  BELL  XV- 15  CHARACTERISTICS  ARE  USED  AS  INPUT  DATA 
D  GROSS  WEIGHT  MIGHT  BE  ONE  OF  THE  COMMENT  STRINGS 

ENTER  CODE  FOR  DATA  INPUT  OR  <RETURN>  TO  CONTINUE  ==> 


FIGURE  4  ROTWASH  PROGRAM  LOGIC/COMMENT  MENU 


typing  <PLT>) .  There  is  no  limit  to  the  number  of  times 
this  switch  can  be  toggled.  As  long  as  the  parameter  has  a 
<Y>  or  "yes"  value,  the  user  must  specify  output  filenames 
before  data  files  are  written  to  disk.  The  user  is  not 
allowed  to  write  over  files  previously  written  to  disk  by 
specification  of  the  same  filename  twice. 

The  last  two  lines  in  the  menu  are  used  to  specify  user 
comments  in  all  data  files  that  are  written  to  disk.  These 
comments  are  also  written  out  as  header  information  on 
screen  output  sent  directly  to  the  printer.  Both  of  the 
comment  lines  can  be  changed  at  any  time  during  program 
execution.  The  only  restriction  is  that  the  character 
strings  on  both  lines  be  less  than  or  egual  to  50 
characters.  The  arrowhead  symbols  above  the  comment  lines 
in  the  menu  define  a  50-space  line  width. 

If  the  velocity  analysis  option  is  specified,  the  user 
is  then  required  to  choose  one  of  the  four  analysis  options 
presented  in  figure  5.  Velocity  analysis  options  reached 
through  this  menu  are  the: 


1.  simple  wall  jet  (for  both  single  and 
twin-rotor  configurations), 

2.  interaction  plane  (twin-rotor  only), 

3.  ground  vortex  (single  rotor  only),  and 

4.  disk  edge  vortex  (single  rotor  only). 
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SELECT  TYPE  OF  FLOW  TO  BE  ESTIMATED 


WALL  JET  PROFILE 

TYPE 

<W> 

INTERACTION  PLANE  PROFILE 

TYPE 

<I> 

GROUND  VORTEX 

TYPE 

<G> 

DISK  VORTEX 

TYPE 

<D> 

TO  EXIT  PROGRAM 

TYPE 

<X> 

ENTER  DATA  ENTRY  CODE 

ll 

ll 

v 

32 

FIGURE  5  VELOCITY  ANALYSIS  OPTION  MENU 


Each  of  these  four  options  is  described  with  flowfield 
characteristics  sketches  in  subsequent  sections.  For 
technical  details  and  discussion  on  practical  limitations  of 
these  options,  the  user  is  referred  to  Sections  III  and  IV 
of  reference  1.  The  choice  of  one  of  these  menu  options  is 
made  by  typing  the  appropriate  code  and  <RETURN> .  If  one  of 
the  five  allowable  characters  is  not  chosen,  the  menu  will 
reappear  and  the  user  will  be  forced  to  choose  an  acceptable 
option . 

If  the  hazard  analysis  option  is  specified,  the  user  is 
presented  with  the  figure  6  menu.  This  menu  allows  the  user 
to  choose  either: 


1.  human  overturning  force/moment  analysis,  or 

2.  particulate  cloud  analysis. 


Both  of  these  analyses  can  be  applied  to  either  single  or 
twin  rotor  configurations.  The  mechanics  of  this  menu 
operate  exactly  like  those  of  the  velocity  analysis  option 
menu . 

2.4  THE  WALL  JET  OPTION 

Rotorwash  velocity  profiles  are  calculated  for  single 
main  rotor  helicopter  configurations  using  the  wall  jet 
option.  Velocity  profiles  along  the  90  and  270  degree 
azimuths  for  twin  rotor  configurations  are  also  calculated 
using  this  option  (90  degrees  is  out  the  right  wing  on  a 
tiltrotor  and  along  the  centerline  of  the  fuselage  for 
tandem  rotor  helicopters).  Figure  7  provides  a  three- 
dimensional  view  of  the  rotorwash  flowfields  associated  with 


SELECT  TYPE  OF  HAZARD 


OVERTURNING  FORCE/MOMENT 
PARTICULATE  CLOUDS 
TO  EXIT  PROGRAM 


TYPE  <M> 
TYPE  <C> 
TYPE  <X> 


ENTER  HAZARD  CODE  ==> 


FIGURE  6  HAZARD  ANALYSIS  OPTION  MENU 


both  rotor  configurations.  Figure  8  provides  a  cross 
sectional  view  of  the  non-dimensionalized  ROTWASH  wall  jet 
velocity  profile  model.  The  program  menu  associated  with 
the  use  of  this  option  is  presented  in  figure  9. 


FIGURE  7  ROTORWASH  FLOWFIELD  CHARACTERISTICS  OF  SINGLE 
AND  TWIN  ROTOR  CONFIGURATIONS  OPERATING  IN 
CLOSE  PROXIMITY  TO  THE  GROUND 
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VELOCITY  PROFILE  STATUS  MENU 


CODE 

PARAMETER 

VALUE 

UNITS 

A 

PROFILE  STATION  POSITION 

50.00 

FT 

B 

VERTICAL  INCREMENT 

1.00 

FT 

C 

MAXIMUM  PROFILE  HEIGHT 

10.00 

FT 

D 

DATA  OUTPUT  FILENAME 

DFRC . PTS 

ENTER  DATA  ENTRY  CODE  OR  <RETURN>  TO  CONTINUE  ==>  A 

PROFILE  STATION  POSITION  =  50.0000 


ENTER  NEW  VALUE  OR  <RETURN>  TO  LEAVE  AS  IS  ==>  60.0 


FIGURE  9  VELOCITY  PROFILE  STATUS  MENU 


The  velocity  profile  status  menu  provides  the  user  the 
option  to  specify  three  parameters  before  calculations  are 
conducted.  The  horizontal  distance  on  the  ground  from  the 
center  of  the  rotor  where  the  velocity  profile  should  be 
calculated  is  the  first  parameter  specified  on  the  menu. 
Figure  10  is  provided  to  graphically  illustrate  this 
geometry  using  the  Bell  XV-15,  which  is  the  more  complex 
example  (flight  test  data  results  associated  with  this 
figure  are  documented  in  reference  3).  To  specify  the 
profile  station  position  for  the  wall  jet  with  the  XV-15 
(270  degree  radial),  the  user  would  measure  the  distance 
from  the  aircraft  centerline  (DFAC)  and  subtract  16.1  feet 
(the  distance  from  the  centerline  to  the  center  of  the 
rotor) .  The  remaining  parameters  to  be  specified  are  the 
vertical  calculation  increment  and  the  maximum  height  above 
ground  level  (AGL)  to  which  the  profile  should  be 
calculated.  The  default  values  for  these  three  parameters 
are  50  feet,  1  foot,  and  10  feet  respectively.  Each  of 
these  parameters  is  input  by  typing  the  appropriate  code, 
<RETURN>,  the  new  input  data  value,  and  <RETURN>  to  end  the 
sequence.  The  last  parameter  on  the  menu  is  the  filename 
for  data  that  is  written  to  disk  if  the  graphics  file  toggle 
switch  is  set  to  <Y> . 

The  wall  jet  velocity  profile,  as  calculated  by  the 
program  using  default  inputs,  is  presented  in  figure  11  when 
the  <RETURN>  key  is  typed  by  itself.  Output  from  the 
analysis  describes  the  shape  of  both  the  mean  and  peak 
velocity  profiles;  an  example  being  presented  in  figure  12 
as  correlated  with  Bell  XV-15  flight  test  data  (this 
particular  example  correlates  to  an  HROTOR  height  of  37.5 
feet  and  a  DFAC  value  of  66.1  feet  on  figure  10). 
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SINGLE  ROTOR  VELOCITY  PROFILE  AT  RADIUS 


50.0  FT 


PROFILE 

BOUNDARY 

HEIGHT 

=  10.35 

FT 

HALF-VEL. 

HEIGHT 

=  3.70 

FT 

MAX-VEL 

HEIGHT 

=  .72 

FT 

EIGHT 

MEAN  VELOCITY 

PEAK  VELOCITY 

MEAN  Q 

(FT) 

(FPS) 

(KN) 

(FPS) 

(KN) 

(PSF) 

.  00 

.  000 

.000 

.  000 

.  000 

.  000 

1.00 

41 . 424 

24 . 555 

77.788 

46.110 

2 . 039 

2 . 00 

33 . 642 

19.942 

70.006 

41 . 497 

1.345 

3.00 

26.396 

15.647 

62 . 760 

37.202 

.828 

4.00 

19.903 

11.798 

56 . 267 

33 . 353 

.471 

5.00 

14 .252 

8 .448 

50.615 

30 . 003 

.241 

6.00 

9.495 

5.629 

45.859 

27 . 184 

.  107 

7 . 00 

5.672 

3.362 

42 .036 

24.918 

.  038 

8.00 

2.811 

1.666 

39.175 

23 . 222 

.  009 

9.00 

.  936 

.  555 

37 . 299 

22 . 110 

.  001 

10.00 

.  064 

.  038 

36.428 

21 . 593 

.  000 

TYPE  <C>ONTINUE,  NEXT  <P>0INT,  <N>EW  CASE,  E<X>IT  ==> 


FIGURE  11  WALL  JET  VELOCITY  PROFILE  OUTPUT  FORMAT 
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BELL  XV- 15  TILTROTOR 
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FIGURE  12  XV- 1 5  VELOCITY  PROFILE  CORRELATION 


The  output  format  provides  velocity  profile  data  in  units  of 
either  feet  per  second,  knots,  or  pounds  per  square  foot 
(also  referred  to  as  dynamic  pressure) . 

If  the  specified  increment  in  vertical  height  is  small 
or  the  maximum  calculated  height  is  large,  the  quantity  of 
data  to  be  output  to  the  screen  may  exceed  the  ten  line 
limit  for  one  screen  frame.  When  this  situation  occurs,  the 
typing  of  <C>  at  the  prompt  (at  the  bottom  of  the  screen) 
results  in  the  next  10  lines  of  data  being  written  to  the 
screen.  If  the  code  value  <P>  is  typed,  the  program  returns 
to  the  wall  jet  analysis  menu.  The  typing  of  code  <N> 
results  in  the  program  returning  to  the  master  input  data 
menu.  If  the  code  <X>  is  typed,  the  program  returns  to  the 
DOS  system  prompt. 

Flight  test  data,  as  correlated  with  output  from  the 
wall  jet  option,  are  presented  in  references  1  and  2  for  the 
Bell  XV- 15 ,  Sikorsky  CH-53E,  and  Canadair  CL-84.  Original 
sources  for  these  data  are  references  3  through  5.  As  based 
on  the  correlation  conducted  for  references  1  and  2,  it  is 
generally  recommended  that  the  wall  jet  option  be  used  for 
calculation  of  velocity  profiles  at  distances  greater  than 
1.5  times  the  rotor  radius  from  the  center  of  the  rotor.  At 
distances  less  than  this  value,  the  mathematical  model  is 
not  detailed  enough  to  accurately  predict  rotorwash 
flowfield  characteristics.  This  1  '  -.itation  is  not  serious 
because  distances  closer  to  th->  rotor  tip  than  1.5  times  the 
rotor  radius  have  litt1°  practical  reason  for  being  analyzed 
for  rotorwash  effects  on  the  environment.  Collision 
avoidance  with  respect  to  objects  in  close  proximity  to  the 
rotorcraft  is  the  critical  i^sne  at  this  close  a  distance. 

2.4  THE  TWIN-ROTOR  INTERACTION  PLANE  OPTION 

The  twin-rotor  interaction  plane  option  calculates  the 
velocity  profile  contai  ied  in  the  plane  which  is  oriented 
perpendicular  to  the  ground  and  to  the  line  segment  which 
connects  the  center  of  both  rotor  hubs  (refer  to  figure  7) . 
After  choosing  the  interaction  plane  option  on  the  velocity 
analysis  option  menu,  the  user  must  specify  the  same 
parameters  on  the  velocity  profile  status  menu  that  are 
required  with  the  wall  jet  option.  The  only  difference 
between  the  parameters  is  the  reference  position  for 
specificatio  of  the  horizontal  location  of  the  velocity 
profile  with  respect  to  the  rotor  (the  first  menu  option) . 
For  the  interaction  plane  analysis,  this  distance  is 
referenced  as  zero  at  the  intersection  of  the  interaction 
plane  and  the  line  connecting  the  rotors  and  not  directly  to 
the  center  of  one  of  the  two  rotors  (i.e.  the  input  value 
for  distance  along  the  interaction  plane  (DAIP)  for  a 
tiltrotor  is  along  a  line  that  is  an  extension  of  the 
fuselage  centerline  as  seen  in  figure  10) .  No  velocity 
profile  differences  are  assumed  to  exist  by  the  mathematical 
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model  for  points  equidistant  along  the  interaction  plane  but 
on  opposite  sides  of  the  line  connecting  the  rotors.  For  a 
tiltrotor,  this  means  that  the  calculated  velocity  profiles 
both  directly  in  front  of  and  directly  aft  of  the  aircraft 
are  the  same  when  equidistant  along  the  interaction  plane. 

Example  output  from  the  interaction  plane  option  is 
presented  in  figure  13.  This  data  format  closely  resembles 
that  of  the  wall  jet  option  except  that  both  horizontal  and 
vertical  velocity  profile  components  are  calculated.  The 
horizontal  velocity  component  is  identified  by  the  "H"  in 
the  column  following  the  height  column  and  the  vertical 
component  by  the  "V" .  The  mechanics  for  viewing  data  on  the 
screen  (if  more  data  exist  than  will  fi*-  on  one  screen 
frame)  and  for  transferring  to  other  me  .  s  are  exactly  as 
are  described  for  the  wall  jet  option. 


TWIN  ROTOR  INTERACTION  PLANE  VELOCITY  PROFILE  AT  STATION  = 


HEIGHT 

MEAN 

VELOCITY 

(FT) 

(FPS) 

(KN) 

.  00 

H 

.  000 

.000 

V 

.  000 

.000 

1.00 

H 

59.870 

35.489 

V 

20.476 

12.137 

2 . 00 

H 

59.072 

35.016 

V 

21.384 

12.676 

3 . 00 

H 

58.251 

34.529 

V 

22 .252 

13.190 

4 . 00 

H 

57.409 

34.030 

V 

23 . 078 

13 . 680 

PEAK 

VELOCITY 

MEAN  Q 

PEAK  Q 

(FPS) 

(KN) 

(PSF) 

(PSF) 

.000 

.  000 

.000 

.000 

.000 

.000 

.000 

.000 

106. 42C 

63 . 086 

4.260 

13.461 

36.398 

21.575 

.498 

1.574 

105.007 

62.245 

4.147 

13.104 

38.012 

22.533 

.543 

1.717 

103 . 546 

61.379 

4.033 

12.742 

39 . 555 

23.447 

.588 

1.859 

102.049 

60.492 

3.917 

12 . 377 

41.024 

24 .318 

.633 

2.000 

TYPE  cC^ONTINUE ,  NEXT  <P>OINT,  <N>EW  CASE,  E<X>IT  ==> 


FIGURE  13  INTERACTION  PLANE  VELOCITY  PROFILE  OUTPUT  FORMAT 


Output  from  the  interaction  plane  option,  as  correlated 
with  XV-15  flight  test  data,  is  presented  in  figure  14  as  an 
example.  These  data  are  excerpted  from  reference  2.  Other 
original  flight  test  data  for  both  the  Bell  XV-15  and 
Canadair  CL-84  (tiltwing)  are  documented  in  references  3  and 
5 . 
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BELL  XV- 15  TILTROTOR 
GW  =  12475  LB,  WAGL  =  25  FT 
AZIMUTH  =  0  DEG,  DAIP  =  34.8  FT 


FIGURE  14  XV- 15  INTERACTION  PLANE  VELOCITY 
PROFILE  CORRELATION 


Several  practical  recommendations  need  to  be  noted  for 
users  of  the  interaction  plane  option.  When  analyzing 
tiltrotor  and  tiltwing  configurations,  the  user  must  be 
careful  to  avoid  choosing  analysis  locations  that  are 
coincident  with  components  of  the  nose  or  tail  structure  on 
the  aircraft.  Also,  locations  in  close  proximity  to  the 
nose  or  tail  of  a  real  aircraft  should  not  be  expected  to 
have  rotorwash  flowfield  characteristics  identical  to 
calculated  velocity  profiles.  At  these  locations  unmodeled 
airframe  aerodynamic  interferences  significantly  influence 
the  rotorwash  flowfield  structure.  Flight  test  data 
obtained  from  the  XV-15  and  CL-84  indicate  that  measured 
mean  and  peak  velocity  profiles  at  points  equidistant  along 
the  interaction  plane  (both  in  front  of  and  aft  of  an  actual 
aircraft)  often  do  not  yield  identical  results  as  might  be 
expected.  Therefore,  output  from  this  analysis  option 
should  be  calibrated  with  flight  test  data  whenever  possible 
in  an  attempt  to  determine  whether  positions  forward  or  aft 
of  the  aircraft  may  be  more  critical  for  analysis. 

2.5  GROUND  VORTEX  ANALYSIS  OPTION 

A  ground  vortex  structure  is  formed  when  ambient  wind 
and/or  rotorcraft  translational  velocity  overcome  the  rotor- 
induced  wall  jet  flowfield.  A  diagram  of  the  ground  vortex 
is  presented  in  figure  15.  Due  to  the  elementary  nature  of 
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GROUND 

VORTEX 


FIGURE  15  DIAGRAM  OF  THE  GROUND  VORTEX  STRUCTURE 


the  mathematical  model  formulation  used  in  the  ROTWASH 
program,  the  ground  vortex  option  should  be  used  with 
caution.  As  discussed  in  reference  1,  almost  no  test  data 
exist  to  validate  the  mathematical  model.  Also,  the  single 
main  rotor  helicopter  is  the  only  configuration  that  can  be 
analyzed  with  the  model  as  presently  formulated.  Since  all 
examples  presented  up  to  this  point  in  the  user's  guide  have 
been  for  the  XV-15,  it  is  necessary  to  define  input  data  for 
a  single  main  rotor  helicopter  before  the  ground  vortex 
option  is  explained  in  detail.  The  Sikorsky  CH-53E 
configuration  has  been  chosen  for  this  task.  The  main  input 
data  menu  for  this  configuration,  as  typed  into  the  program, 
is  presented  in  figure  16. 

The  ground  vortex  analysis  option  is  specified  by 
selection  of  the  character  <G>  on  the  velocity  analysis 
option  menu  as  shown  in  figure  5.  The  user  is  then  required 
to  complete  the  ground/disk  vortex  input  data  menu  which  is 
presented  in  figure  17.  Two  of  the  parameters  specified  on 
this  menu  are  rotorcraft  configuration  parameters.  These 
parameters  are  the  rotor  tip  speed  (ft/sec)  and  the  number 
of  rotor  blades.  The  next  parameter  to  be  specified  is  the 
rotorcraft  translational  velocity  with  respect  to  the 
surrounding  air  mass  (i.e.,  an  input  of  15  knots  can  be 
either  15  knots  ground  speed  on  a  no-wind  day  or  0  knots 
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ROTWASH  USER  INPUT  DATA  MENU 


CODE  PARAMETER  VALUE  UNITS 


A 

NUMBER  OF  ROTORS  (1  OR  2) 

1 

-ND- 

B 

HUB  TO  HUB  ROTOR  SEPARATION 

.0 

FT 

C 

ROTOR  RADIUS 

39.5 

FT 

D 

GROSS  WEIGHT 

56000.0 

LB 

E 

FUSELAGE  DOWNLOAD  FACTOR 

5.0 

PCT 

F 

ROTOR  HEIGHT  ABOVE  GROUND 

30.0 

FT 

G 

SHAFT  TILT  ANGLE  (<20  DEG) 

.0 

DEG 

H 

AIR  DENSITY  RATIO 

1.0000 

ND 

I 

AMBIENT  WIND  (-10  TO  10  KT) 

.  0 

KT 

ENTER  DATA  ENTRY  CODE  OR  <RETURN>  TO  CONTINUE  ==> 


FIGURE  16  CH-53E  INPUT  DATA  FOR  THE  GROUND  VORTEX  EXAMPLE 


GROUND/DISK  VORTEX  INPUT  DATA  MENU 
(FOR  SINGLE  MAIN  ROTOR  HELICOPTERS  ONLY) 


CODE 

PARAMETER 

VALUE 

UNITS 

A 

ROTOR  TIP  SPEED 

600.00 

FPS 

B 

NUMBER  OF  ROTOP.  BLADES 

7.00 

-ND- 

C 

TRANSLATIONAL  SPEED 

16.00 

FPS 

D 

XT  POSITION 

-80.00 

FT 

E 

YT  POSITION 

.  00 

FT 

F 

ZT  CALCULATION  INCREMENT 

2.00 

FT 

G 

MAXIMUM  CALCULATION  HEIGHT 

10.00 

FT 

ENTER 

DATA  ENTRY  CODE  OR  <RETURN>  TO 

CONTINUE  = 

=> 

FIGURE  17  GROUND/DISK  VORTEX  INPUT  DATA  MENU 
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ground  speed  on  a  day  with  a  15  knot  headwind) .  Each  of 
these  values  is  input  with  the  keyboard  using  the  same 
techniques  that  have  been  previously  discusssed. 

The  next  four  menu  parameters  define  the  position  in 
three-dimensional  space  (feet)  where  the  velocity  profile 
will  be  calculated  (see  figure  18) .  The  positive  directions 
for  the  coordinate  system  are  aft  and  right  from  the  center 
of  the  rotor.  Therefore,  in  order  to  calculate  a  slice  of 
the  ground  vortex  directly  in  front  of  the  rotor,  a  negative 
X-value  (longitudinal  position)  is  input  along  with  a  zero 
Y-value  (lateral  position).  The  Z-axis  increment  (feet)  and 
maximum  calculation  height  (feet)  parameters  define  the 
number  of  points  that  are  calculated  between  the  ground 
level  and  the  maximum  height  of  interest. 

After  each  input  data  parameter  has  been  specified  the 
user  types  the  <RETURN>  key  to  initiate  the  analysis.  The 
response  of  the  program  is  to  list  three  calculated 
parameters  which  are  followed  by  a  prompt.  An  example  of 
this  response  is  provided  in  figure  19.  The  calculated 
values  are  the  nondimensionalized  rotor  height  above  ground 
and  two  advance  ratio  parameters.  The  program  prompt 
following  the  screen  output  requires  the  user  to  input  the 
ground  vortex  strength  ratio  which  is  obtained  from  the 
graph  in  figure  20  using  the  three  calculated  parameters. 
Background  on  the  use  of  this  figure  is  discussed  in  Section 
III  of  reference  1.  The  limits  presented  on  the  graph  in 
figure  20  define  the  advance  ratio  range  wherein  the  ground 
vortex  would  be  expected  to  occur.  At  advance  ratios  much 
less  than  0.035,  the  ground  vortex  does  not  have  favorable 
conditions  for  formation.  At  advance  ratios  slightly 
greater  than  0.055,  the  ground  vortex  is  dispersed  by  the 
rotor  because  the  translational  velocity  relative  to  the  air 
mass  is  too  high  for  the  vortex  to  maintain  position.  After 
the  user  has  entered  the  ground  vortex  strength  ratio,  the 
program  calculates  the  ground  vortex  circulation  and  the 
position  of  the  ground  vortex  core  with  respect  to  the  axis 
system  presented  in  figure  18.  This  output  is  also 
presented  in  figure  19  below  the  prompt  for  the  ground 
vortex  strength  ratio. 

Engineering  data  from  the  ground  vortex  analysis  option 
is  obtained  by  typing  the  <RETURN>  key  after  the  vortex 
position  parameters  are  displayed.  An  example  of  the  output 
format  is  presented  in  figure  21.  Calculated  field 
velocities  at  the  various  points  along  the  profile  Z-axis 
are  presented  in  both  a  vectorial  XYZ  component  format  and 
as  a  total  resolved  magnitude  in  both  units  of  feet  per 
second  and  knots.  These  same  data  are  also  provided  to  the 
user  as  dynamic  pressures  in  units  of  pounds  per  square 
feet. 
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GROUND  VORTEX 


RELATIVE  WIND 


FIGURE  18  HORSE  SHOE  VORTEX  GEOMETRY  FOR  THE 
GROUND  VORTEX  ANALYSIS  OPTION 


ROTOR  HEIGHT  ABOVE  GROUND  H/D  .3797 
ADVANCE  RATIO  MU-STAR  .5373 
ADVANCE  RATIO  MU  .0450 


ENTER  GROUND  VORTEX  STRENGTH  RATIO 

(SEE  FIGURE  18)  ==>  4.0 


GROUND  VORTEX  CORE  POSITION 

X-LOCATION  (XXGV)  =  -69.51  FT 

Y- LOCATION  ( ZZGV)  =  5.93  FT 

GROUND  VORTEX  CIRCULATION  =  1192.92  FT**2/SEC 


PRESS  <RETURN>  TO  CONTINUE 


FIGURE  19  GROUND  VORTEX  ANALYSIS  INTERMEDIATE  OUTPUT 


C-O 


on 
f — • 

l/> 


X 


on 

o 


o 


o 

on 

o 


ADVANCE  RATIO,  JJ 


FIGURE  20  GROUND  VORTEX  CIRCULATION  STRENGTH 
AS  A  FUNCTION  OF  ADVANCE  RATIO 


HEIGHT 

MEAN 

VELOCITY 

MEAN  Q 

(FT) 

(FPS) 

(KN) 

(PSF) 

.  00 

X 

-14 .839 

-8 . 796 

.262 

Y 

.  000 

.  000 

.000 

Z 

.000 

.  000 

.  000 

T 

14 .839 

8 . 796 

.262 

2.00 

X 

-13 . 988 

-8.291 

.233 

Y 

.000 

.  000 

.  000 

Z 

4 .257 

2.523 

.  022 

T 

14.621 

8 . 667 

.254 

4 . 00 

X 

-11.601 

-6 .877 

.  160 

Y 

.  000 

.  000 

.  000 

Z 

7.768 

4.605 

.  072 

T 

13.961 

8 .276 

.232 

TYPE  <C>0NTINUE,  NEXT  <P>0INT,  <N>EW  CASE,  E<X>IT  ==> 


FIGURE  21  GROUND  VORTEX  VELOCITY  FIELD  OUTPUT  DATA 
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2.6  DISK  EDGE  VORTEX  ANALYSIS  OPTION 


The  disk  edge  vortex  analysis  option  was  developed  to 
provide  a  capability  to  estimate  the  strength  of  trailing 
vortices  behind  helicopters  in  forward  flight  as  described 
by  figure  22.  Like  the  ground  vortex  option,  this  option  is 
limited  to  use  with  the  single  main  rotor  helicopter 
configuration.  Approximations  of  vortex  core  size  are  not 
calculated  by  the  mathematical  model  and  must  be  estimated 
using  flight  test  data.  Limited  data  are  presented  in 
section  IV  of  reference  1.  However,  the  FAA  has  recently 
conducted  flight  tests  to  measure  additional  data. 

Therefore,  the  interested  researcher  should  be  alert  for 
future  FAA  publications  documenting  these  efforts.  All 
examples  presented  in  the  user's  guide  for  this  analysis 
option  utilize  the  same  CH-53E  input  data  array  that  was 
described  previously. 


BOUND  VORTEX 


FIGURE  22  DISK  EDGE  VORTEX  GEOMETRY  DEFINITION 


Execution  of  the  disk  edge  vortex  option  ls  initiated 
by  choosing  <D>  on  the  velocity  analysis  opt  ion  menu  shown 
in  figure  5.  The  ground/disk  vortex  input  data  menu,  as 
shown  in  figure  17,  must  then  be  completed  as  was  discussed 
in  the  previous  section.  All  sign  conventions  used  in  the 
specification  of  locations  for  the  calculation  of  field 
velocities  are  the  same  as  for  the  ground  vortex  option. 
After  the  input  data  menu  is  completed,  the.  analysis  option 
is  executed  by  typing  the  <RETURN>  key. 


DISK  VORTEX  VELOCITY  PROFILE  DATA 


X-LOCATION  (XT)  =  150.00  FT 

Y-LOCATION  (YT)  =  . 00  FT 


VORTEX  CIRCULATION 
SETTLING  ANGLE 


3712.40  FT**2/SEC 
9.29  DEG 


PRESS  <RETURN>  TO  CONTINUE 


FIGURE  23  DISK  EDGE  VORTEX  OPTION  INTERMEDIATE  OUTPUT 


The  initial  program  response,  figure  23,  is  to  write  to 
the  screen  the  calculated  values  for  the  vortex  circulation 
and  the  settling  angle  of  the  trailing  vortex  components  as 
defined  in  figure  22.  When  executing  this  option,  it  is 
important  that  the  user  confirm  that  the  calculated  settling 
angle  is  less  than  approximately  45  degrees.  At  settling 
angles  larger  than  this  value,  the  airspeed  of  the 
helicopter  is  probably  too  slow  to  sustain  the  formation  of 
the  trailing  edge  vortex  system  which  is  predicted  using 
this  mathematical  model,  and  results  should  be  considered 
suspect.  User  specified  velocity  field  calculations  are 
presented  after  the  <RETURN>  key  is  typed  in  the  same  format 
as  was  discussed  with  the  ground  vortex  option.  An  example 
output  for  the  CH-53E  is  presented  in  figure  24. 

Other  features  of  the  mathematical  model  which  should 
be  noted  include: 


1.  The  calculated  trailing  vortex  strength  does 
not  decay  as  a  function  of  increasing  distance 
from  the  helicopter  or  with  increasing  time. 

This  weakness  will  affect  correlation  with 
flight  test  data  and  could  result  in  greater- 
than-predicted  values  for  field  velocities. 

2.  The  decay  of  the  trailing  vortex  structure  after 
impingement  with  the  ground  is  not  modeled.  Any 
prediction  of  field  velocities  signif icantly 
behind  the  initial  impingement  point  should  not 
be  considered  valid.  The  location  of  the 
impingement  point  with  respect  to  the  ground  and 
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HEIGHT 

MEAN 

VELOCITY 

MEAN  Q 

(FT) 

(FPS) 

(KN) 

(PSF) 

.  00 

X 

-.996 

-.590 

.  001 

Y 

.  000 

.  000 

.  000 

Z 

.000 

.  000 

.  000 

T 

.996 

.  590 

.  001 

50.00 

X 

-1.113 

-.660 

.001 

Y 

.  000 

.000 

.000 

Z 

-1.933 

-1.146 

.004 

T 

2.230 

1.322 

.006 

100.00 

X 

-1.651 

-.979 

.003 

Y 

.000 

.  000 

.  000 

Z 

-6.133 

-3.636 

.045 

T 

6.352 

3.765 

.048 

TYPE  COONTINUE,  NEXT  <P>OINT,  <N>EW  CASE,  E<X>IT  ==> 


FIGURE  24  VELOCITY  FIELD  OUTPUT  DATA  FROM 
THE  DISK  EDGE  VORTEX  OPTION 


the  location  of  points  specified  for  velocity  field 
calculations  must  be  checked  by  hand  calculation  to 
ensure  that  geometry  constraints  are  not  violated 
This  is  accomplished  by  using  the  rotor  height  and 
the  settling  angle  to  calculate  the  horizontal 
distance  behind  the  helicopter  where  the  impingement 
occurs . 


2.7  PERSONNEL  OVERTURNING  FORCE  AND  MOMENT  ANALYSIS 

The  personnel  overturning  force  and  moment  analysis 
model  is  formulated  for  use  with  both  the  single  main  rotor 
and  twin-rotor  configurations.  The  initial  task  of  the 
model  is  to  calculate  the  velocity  profile  for  a  specified 
location.  The  calculated  velocity  profile  is  then 
integrated  over  the  projected  area  of  a  human  body  to  obtain 
estimates  of  the  applied  aerodynamic  force  and  moment.  This 
analysis  technique  is  summarized  by  figure  25. 

Use  of  the  overturning  force  and  moment  option  is 
initiated  by  choosing  <H>  on  the  program  logic/comment  menu. 
This  is  followed  by  the  choice  of  <M>  on  the  hazard  analysis 
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STEPS: 


(1)  Calculate  the  dynamic  pressure  at  twelve  vertical  stations. 

(2)  Use  the  dynamic  pressure  to  calculate  twelve  force  values  (multiply 
each  Qx  value  by  the  projected  body  surface  area  at  that  7^  value, 
In  this  case  0.5  ft  by  1.1  ft). 

(3)  Sum  the  twelve  Individual  Fx  values  in  order  to  obtain  the  total 
force  (Fpeak). 

(A)  Multiply  each  Fx  value  by  its  corresponding  Zx  value,  and  sum  the 
incremental  overturning  moments  to  obtain  the  total  peak  moment 

<Mpeak>* 


FIGURE  25  OVERTURNING  FORCE  AND  MOMENT 
CALCULATION  PROCEDURES 
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option  menu  as  shown  in  figure  6.  The  user  then  specifies 
the  parameters  listed  on  the  overturning  force  and  moment 
data  menu  as  presented  in  figure  26.  The  first  parameter  on 
this  menu  specifies  the  use  of  either  the  wall  jet  or  the 
interaction  plane  analysis  for  creation  of  velocity  profile 
data.  The  second  option  specifies  the  use  of  either  the 
"large"  (six  feet  in  height)  or  "small"  (four  feet  in 
height)  human  body  mathematical  model.  The  third  parameter 
provides  the  user  the  capability  to  specify  a  graphics 
output  filename  (assuming  this  option  has  been  toggled  ON 
using  the  program  logic/comment  menu) .  If  the  user  executes 
the  option  without  changing  the  filename  and  a  file  already 
exists  with  the  same  filename,  the  user  is  notified  and 
required  to  change  the  filename.  The  last  three  menu 
variables  define  the  locations  that  are  to  be  analyzed  using 
the  option.  These  variables,  all  in  units  of  feet,  are  the 
initial  station  position  for  analysis,  the  increment  in 
station  position,  and  the  final  station  position 
respectively.  The  mechanics  for  input  of  the  desired  values 
using  this  menu  are  as  described  for  previous  menus. 


CODE 


OVERTURNING  FORCE/MOMENT  DATA  MENU 


PARAMETER 


VALUE  UNITS 


A 

B 

C 


<W>ALL  JET  OR  < INTERACTION  PLANE 
<L>ARGE  OR  <S>MALL  PERSON 
DATA  OUTPUT  FILENAME 


W 

L 

OTDFRC.PTS 


D  INITIAL  STATION  POSITION 

E  HORIZONTAL  INCREMENT 

F  MAXIMUM  STATION  POSITION 


50.00  FT 

10.00  FT 

100.00  FT 


ENTER  DATA  ENTRY  CODE  OR  <RETURN>  TO  CONTINUE  ==> 


FIGURE  26  OVERTURNING  FORCE  AND  MOMENT  DATA  MENU 


After  the  analysis  is  executed  by  typing  the  <RETURN> 
key,  the  calculated  results  are  written  out  in  the  format  as 
presented  in  figure  27.  Three  columns  of  data  are  written 
using  this  format.  The  first  column  identifies  either  the 
distance  from  rotor  center  (DFRC)  for  the  wall  jet  option  or 
the  distance  along  the  interaction  plane  (DAIP)  for  the 
interaction  plane  option.  The  second  and  third  columns  are 
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SUMMARY  OF  OVERTURNING  FORCES  AND  MOMENTS 


RADIUS 

TOTF 

TOTM 

(FT) 

(LB) 

(FT-LB) 

50.00 

35.185 

86.233 

60.00 

32.114 

83 . 975 

70.00 

26.383 

71.428 

80 . 00 

21.770 

60.359 

90 . 00 

18 . 094 

51.075 

100.00 

15.168 

43.438 

TYPE  cOONTINUE,  NEXT 

<P>OINT , 

<N>EW  CASE 

FIGURE  27  SIMPLIFIED  OUTPUT  FORMAT  FOR  THE  OVERTURNING 
FORCE  AND  MOMENT  ANALYSIS 


the  associated  total  force  and  total  moment  values  as 
calculated  for  the  projected  area  of  a  human  body.  Example 
data  for  the  XV-15  using  this  option  are  presented  in  figure 
28  for  reference.  At  the  bottom  of  the  screen,  the  user  is 
required  to  return  to  the  previous  menu  by  typing  <P>,  the 
master  input  data  menu  by  typing  <N> ,  or  to  exit  the  program 
by  typing  <X> .  By  typing  <C>,  the  user  can  view  the  next 
video  screen  of  data  if  more  than  one  screen  of  data  are 
generated.  Otherwise,  the  <C>  option  works  identically  to 
the  <P>  option. 

If  the  user  needs  to  examine  the  detailed  calculations 
used  to  create  the  summary  output,  he  specifies  the  desired 
analysis  location  as  the  initial  station  position.  This 
input  is  followed  by  specification  of  the  increment  value  as 
0.0  (as  shown  in  figure  29).  The  resulting  "large"  person 
output  using  this  option  is  presented  on  two  screen  frames 
as  shown  in  figure  30.  The  first  output  frame  presents  a 
summary  of  the  velocity  profile  calculations  as  a  function 
of  height  above  ground  at  the  specified  station  position. 

The  second  output  frame  presents  the  associated  calcultions 
for  dynamic  pressure,  overturning  force,  and  overturning 
moment.  The  last  two  columns  in  the  second  table  are  values 
of  total  force  and  moment  as  summed  for  the  incremental 
increase  in  height. 


OVERTURNING  FORCE,  POUNDS 


BELL  XV-15 


DISTANCE  ALONG  INTERACTION  PLANE  (DAIP),  FEET 


FIGURE  28  BELL  XV-15  OVERTURNING  FORCE  AS  A  FUNCTION 
OF  DISTANCE  ALONG  THE  INTERACTION  PLANE 
(0  DEGREE  AZIMUTH) 
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OVERTURNING  FORCE/MOMENT  DATA  MENU 


CODE  PARAMETER 


VALUE  UNITS 


A 

B 

C 


<W>ALL  JET  OR  < INTERACTION  PLANE 
<L>ARGE  OR  <S>MALL  PERSON 
DATA  OUTPUT  FILENAME 


W 

L 

OTDFRC . PTS 


D  INITIAL  STATION  POSITION 

E  HORIZONTAL  INCREMENT 

F  MAXIMUM  STATION  POSITION 


50.00  FT 

.  00  FT 

110.00  FT 


ENTER  DATA  ENTRY  CODE  OR  <RETURN>  TO  CONTINUE  ==> 


FIGURE  29  MENU  SPECIFICATION  OF  DETAILED 
FORCE/MOMENT  OUTPUT 


In  this  example,  the  force  and  moment  values  of  35.2 
pounds  and  86.2  foot-pounds  respectively  at  5.75  feet  are 
the  total  force  and  moment  values  as  would  normally  be 
printed  out  in  the  summary  output  (figure  27).  These  values 
are  checked  by  totaling  the  individual  height  related  values 
in  the  overturning  force  and  moment  columns  (second  and 
third  columns) .  An  example  output  for  the  second  screen  of 
the  "small"  person  option  is  presented  in  figure  31  for 
reference . 

Both  qualitative  and  quantitative  overturning  force  and 
moment  data  are  presented  in  reference  1  as  correlated  with 
ROTWASH  program  output  for  the  Bell  XV-15,  Sikorsky  CH-53E, 
and  the  Sikorsky  S-61.  These  calculated  data  all  assume  a 
coefficient  of  drag  for  a  human  body  of  1.1  (which  according 
to  Hoerner,  reference  6,  can  vary  from  1.0  to  approximately 
1.3)  . 

2.8  PARTICULATE  CLOUD  ANALYSIS  OPTION 

The  methodology  used  in  the  calculation  of  particulate 
cloud  size  is  presented  in  Section  V  of  reference  1  along 
with  a  very  limited  amount  of  flight  test  data.  This  option 
is  applicable  to  both  single  main  rotor  and  twin-rotor 
configurations.  The  particulate  cloud  geometry  utilized  in 
the  analysis  option  is  presented  in  figure  32. 
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SINGLE  ROTOR  VELOCITY  PROFILE  AT  RADIUS 


50.0  FT 


MEAN  VELOCITY  PEAK  VELOCITY 

(FPS)  (KN)  (FPS)  (KN) 


HEIGHT 

(FT) 

.25 

.75 

1.25 

1.75 

2.25 

2.75 

3.25 

3.75 

4.25 

4.75 

5.25 

5.75 


37.456 
43.341 

39.457 
35.553 
31.768 
28.142 
24.698 
21.450 
18 .409 
15.582 
12.977 
10.598 


22.203 
25.691 
23.389 
21.075 
18.831 
16.682 
14 . 640 
12.715 
10.912 
9.237 
7.692 
6.282 


73.819 
79 . 704 

75.820 
71.917 
68 . 132 
64 . 506 
61.062 
57 .814 
54.773 
51.946 
49.341 
46.962 


43 . 758 
47.246 
44 . 944 
42 . 630 
40.386 
38.237 
36 . 195 
34.270 
32.467 
30.792 
29.248 
27 . 837 


TYPE  <RETURN>  TO  CONTINUE 


MEAN  Q  PEAK  Q 
( PSF)  (PSF) 

1.667  6.476 

2.232  7.550 

1.850  6.832 

1.502  6.147 

1.199  5.517 

•941  4.945 

.725  4.431 

.547  3.972 

•403  3.565 

.289  3.207 

.200  2.893 

.133  2.621 


SINGLE  ROTOR  FORCE  PROFILE  AT  RADIUS 


50.0  FT 


HEIGHT 

(FT) 


.  4  J 

.75 

1.25 

1.75 

2.25 

2.75 

3.25 

3.75 

4 . 25 

4.75 

5.25 

5.75 


PEAK  Q 
(PSF) 

6.476 
7.550 
6.832 
6.147 
5.517 
4.945 
4.431 
3 . 972 
3 . 565 
3.207 
2.893 
2 . 621 


FOVER 

(LB) 

3.918 
4 . 568 
4.133 
3.719 
3 .338 
2.992 
2 . 681 
2.403 
2  .  1  57 
1.940 
1.750 
1.586 


OVERM 

(FT-LB) 

.980 

3.426 

5.167 
6.508 
7 . 510 
8 .227 
8.713 
9 .012 

9 . 167 
9.216 
9.190 
9 . 118 


TOT  F 
(LB) 

3.918 
8 .486 
12.619 
16.338 
19.676 
22 . 667 
25.348 
27.752 
29 . 909 
31.849 
33.599 
35.185 


TOT  M 
(FT-LB) 

.980 
4.405 
9.572 
16.080 
23.589 
31.817 
40.530 
49 . 542 
58 . 709 
67.925 
77.115 
86.233 


TYPE  <C>ONTINUE ,  NEXT  <P>OINT,  <N>EW  CASE,  E<X>IT  ==> 


FIGURE  30  DETAILED  OUTPUT  FORMAT  FOR  THE  "LARGE" 
PERSON  ANALYSIS  OPTION 
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SINGLE  ROTOR  FORCE  PROFILE  AT  RADIUS 


50.0  FT 


HEIGHT 

PEAK  Q 

FOVER 

OVERM 

TOT  F 

TOT  M 

(FT) 

(  PSF) 

(LB) 

(FT-LB) 

(LB) 

(FT-LB) 

.  25 

6.476 

2.493 

.  623 

2.493 

.  623 

.75 

7.550 

2.907 

2 . 180 

5.400 

2 .803 

1.25 

6.832 

2.630 

3 .288 

8.030 

6.091 

1.75 

6.147 

2.366 

4 . 141 

10.397 

10.233 

2.25 

5.517 

2.124 

4 . 779 

12 . 521 

15.011 

2.75 

4.945 

1.904 

5.236 

14.425 

20.247 

3.25 

4.431 

1.706 

5 . 545 

16.131 

25.792 

3.75 

3.972 

1.529 

5.735 

17.660 

31.527 

TYPE 

<OONTINUE , 

NEXT  <P>OINT 

,  <N>EW 

CASE,  E<X>IT 

==> 

FIGURE  31  DETAILED  OUTPUT  FORMAT  FOR  THE  SECOND 
SCREEN  FRAME  OF  THE  "SMALL"  PERSON 
ANALYSIS  OPTION 


FIGURE  32  PARTICULATE  CLOUD  ANALYSIS  GEOMETRY  DEFINITION 
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The  particulate  cloud  option  is  initiated  with  the 
typing  of  <C>  on  the  hazard  analysis  option  menu  as  shown  in 
figure  6.  The  screen  that  is  subsequently  written  presents 
the  user  with  a  prompt  for  input  of  the  terrain  erosion 
factor.  The  value  for  this  factor  is  chosen  from  the  graph 
in  figure  33. 

The  output  format  for  the  particulate  cloud  option  is 
designed  for  both  single  main  rotor  and  twin-rotor  aircraft 
as  shown  in  figure  34.  In  this  example,  using  the  XV-15, 
the  particulate  cloud  boundaries  located  at  the  90  and  270 
degree  azimuths  (out  the  span  of  the  wing)  are  specified  by 
the  single  rotor  or  "SR"  row  of  output  values.  The 
interaction  plane  or  "IP"  row  defines  the  cloud  boundaries 
which  exist  straight  out  in  front  of  and  directly  aft  of  the 
aircraft  along  the  aircraft  centerline.  If  a  single  main 
rotor  configuration  is  being  evaluated,  the  "IP"  row  in  the 
printout  will  contain  all  zeros. 
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Diameter.,  (p) 
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General  Terrain 


ENTER  TERRAIN  EROSION  FACTOR  (-ND-)  ==>  0.4 


SUMMARY  OF  CLOUD  BOUNDARIES 


RC  AND 

RV  ARE  FROM 

ROTOR  CENTER 

(FT) 

RC 

RV 

ZV 

HC 

SR 

87.0 

68 . 3 

28.6 

40.9 

IP 

122.8 

96.4 

40.4 

57.7 

QSMAX  =  12.7  PSF 


TYPE  <C>ONTINUE,  NEXT  <P>OINT,  <N>EW  CASE,  E<X>IT  ==> 


FIGURE  34  PARTICULATE  CLOUD  ANALYSIS  OPTION  OUTPUT 
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3.0  ROTWASH  PROGRAM  DATA  OUTPUT  FILE  FORMATS 


Four  data  output  file  formats  can  be  specified  from  the 
ROTWASH  program  for  use  with  computer  graphics  programs. 

Two  of  these  output  file  formats  are  generated  by  the  wall 
jet  and  interaction  plane  velocity  profile  analysis  options. 
The  other  two  formats  are  generated  by  the  personnel 
overturning  force  and  moment  option.  These  file  formats 
save  the  summary  force  and  moment  data  for  both  the  wall  jet 
and  interaction  plane  cases. 

The  first  two  lines  in  each  of  the  four  file  formats 
are  user-specified  comments.  These  two  comment  lines  are 
typed  in  through  use  of  the  program  logic/comment  menu.  The 
rest  of  the  data  in  each  of  the  file  formats  is  either 
header  information  or  engineering  data.  The  example  files 
presented  in  this  section  are  written  for  direct  input  to 
the  TECPLOT  Graphics  Program  which  is  written  by  AMTEC 
Engineering.  This  graphics  program  is  one  of  several  IBM 
PC/PC-compatible  engineering  graphics  programs  presently  on 
the  market.  If  the  user  does  not  have  access  to  this 
program,  he  can  easily  modify  the  ROTWASH  FORTRAN  code  for 
other  types  of  graphics  programs  by  modifying  the 
appropriate  write  statements.  Table  1  provides  the  user 
with  a  cross  reference  of  the  figure  number  for  each  of  the 
four  output  types  and  the  source  location  for  the  associated 
FORTRAN  code  that  can  be  modified  (see  program  listings  in 
Appendix  C) . 
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TABLE  1  GRAPHICS  FILE  /  SOURCE  CODE  REFERENCE  MATRIX 


FILE  OUTPUT  TYPE 

FIGURE 

SOURCE  CODE 

LOCATION 

Wall  Jet  Velocity 

Profile  Output 

35 

Subroutine 

WJVEL 

Interaction  Plane 

Velocity  Profile  Output 

36 

Subroutine 

IPVEL 

Overturning  Force/Moment 
Summary  (Wall  Jet) 

37 

Subroutine 

HWJVEL 

Overturning  Force/Moment 
Summary  (Interaction  Plane) 

38 

Subroutine 

HIPVEL 

36 


BELL  XV- 15  CHARACTERISTICS  ARE  USED  AS  INPUT  DATA 
GROSS  WEIGHT  MIGHT  BE  ONE  OF  THE  COMMENT  STRINGS 


TITLE=" SINGLE  ROTOR,  DFRC  =  50.00  FT,  GW  =  13500  LB,  RAGL  =  25  FT" 
VARIABLES  =  X,HT 

ZONE  T  =  "MEAN  PROFILE,  KTS",  1=11,  F=POINT 

.0  .0 

24.6  1.0 

19.9  2.0 

15.6  3.0 

11.8  4.0 

8.4  5.0 

5.6  6.0 

3.4  7.0 

1.7  8.0 

.6  9.0 

.0  10.0 

ZONE  T  =  "PEAK  PROFILE,  KTS",  1=11,  F=POINT 

.0  .0 

46.1  1.0 

41.5  2.0 

37.2  3.0 

33.4  4.0 

30.0  5.0 

27.2  6.0 

24.9  7.0 

23.2  8.0 

22.1  9.0 

21.6  10.0 

ZONE  T  =  "PEAK  Q,  PSF",  1=11,  F=POINT 

.0  .0 

7.2  1.0 

5.8  2.0 

4.7  3.0 

3.8  4.0 

3.0  5.0 

2.5  6.0 

2.1  7.0 

1.8  8.0 

1.7  9.0 

1.6  10.0 


FIGURE  35  EXAMPLE  WALL  JET  OPTION  GRAPHICS  FILE  FORMAT 
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BELL  XV- 15  CHARACTERISTICS  ARE  USED  AS  INPUT  DATA 
GROSS  WEIGHT  MIGHT  BE  ONE  OF  THE  COMMENT  STRINGS 


TITLE="TWIN  ROTOR,  DAIP 
VARIABLES  =  X,HT 


ZONE  T 
.  0 

35.5 
35.0 

34.5 
34.0 

33 . 5 
33 . 0 

32 . 5 

31.9 

31.4 

30.9 
ZONE  T 

.0 

63 . 1 

62.2 

61.4 

60.5 

59.6 

58.7 

57 . 7 

56.8 

55.8 

54 . 9 
ZONE  T 

.  0 

13 . 5 

13.1 
12 . 7 

12.4 
12 . 0 

11 . 6 
11.3 

10.9 

10.5 

10.2 


=  "MEAN 
.0 
1.0 
2 . 0 
3.0 
4.0 
5.0 
6.0 
7 . 0 
8 . 0 
9.0 
10.0 
=  "PEAK 
.  0 
1.0 
2 . 0 
3 . 0 
4 . 0 
5.0 
6.0 
7.0 
8 . 0 
9.0 
10.0 
=  "PEAK 
.0 
1.0 
2 . 0 
3 . 0 
4 . 0 
5.0 
6.0 
7.0 
8.0 
9.0 
10.0 


PROFILE, 


PROFILE, 


Q,  PSF" , 


=  50.00  FT,  GW  =  13500  LB, 
KTS " ,  1=11,  F=POINT 


KTS",  1=11,  F= POINT 


1=11,  F=POINT 


RAGL 


2  5  FT" 


FIGURE  36  EXAMPLE  INTERACTION  PLANE  OPTION  GRAPHICS 
FILE  FORMAT 
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BELL  XV-15  CHARACTERISTICS  ARE  USED  AS  INPUT  DATA 
GROSS  WEIGHT  MIGHT  BE  ONE  OF  THE  COMMENT  STRINGS 


TITLE=" SINGLE  ROTOR  DFRC  DATA” 
VARIABLES  =  DFRC , TOTF , TOTM 
ZONE  T  =  "GW  =  13500  LB,  RAGL  =  25  FT", 
50.00  35.18  86.23 

60.00  32.11  83.97 

70.00  26.38  71.43 

80.00  21.77  60.36 

90.00  18.09  51.08 

100.00  15.17  43.44 

110.00  12.84  37.20 


1=7,  F=POINT 


FIGURE  37  EXAMPLE  PERSONNEL  OVERTURNING  FORCE  AND 
MOMENT  GRAPHICS  FILE  FORMAT  CREATED  WITH 
THE  WALL  JET  ANALYSIS  OPTION 
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BELL  XV-15  CHARACTERISTICS  ARE  USED  AS  INPUT  DATA 
GROSS  WEIGHT  MIGHT  BE  ONE  OF  THE  COMMENT  STRINGS 


TITLE="TWIN  ROTOR  DAIP  DATA" 

VARIABLES  =  DAIP, TOTF, TOTM 

ZONE  T  =  "GW  =  13500  LB,  RAGL  =  25  FT",  1=7,  F=POINT 


50. 

.00 

90. 

.  16 

268. 

.97 

60. 

.00 

79. 

.31 

239. 

,89 

70. 

,  00 

70. 

.66 

215. 

.84 

80. 

,00 

63. 

,81 

196. 

,36 

90. 

,00 

55. 

.54 

172. 

.08 

100. 

.00 

44. 

.  15 

137. 

,68 

110. 

,  00 

35. 

,78 

112. 

.18 

FIGURE  38  EXAMPLE  PERSONNEL  OVERTURNING  FORCE  AND 
MOMENT  GRAPHICS  FILE  FORMAT  CREATED  WITH 
THE  INTERACTION  PLANE  ANALYSIS  OPTION 
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SUBJECT:  Letter  Report,  Oblique  Impingement  of  Rotor 

Downwash  with  the  Ground  Plane 

REFERENCES:  (1)  Ferguson,  S.W.,  "Analysis  and  Recommendation 

of  Separation  Requirements  for  Rotoreraft 
Operation  at  Airports  and  Heliports",  Systems 
Technology,  Inc.,  STI  Report  TR-1224-1, 
September  1986. 
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(3)  Watts,  A.,  "V/STOL  Downwash  Impingement  Study, 
Velocity  Estimate",  Canadair  Limited, 
Memorandum  MAA-284-001,  April  1971. 


BACKGROUND 


A  methodology  which  predicts  the  rotorwash 
characteristics  of  hovering  helicopter,  tiltrotor,  and 
tiltwing  aircraft  was  developed  in  1986  in  reference  1.  This 
analysis  methodology  is  fundamentally  based  upon  momentum 
theory  and  the  assumption  (which  is  acceptable  for  most 
scenarios)  that  the  hovering  rotor  is  contained  within  a  plane 
parallel  to  the  surface  of  the  ground.  Two  scenarios  where 
this  assumption  does  not  apply  are  the  initial  segment  of  a 
departure  from  hover  and  the  final  segment  of  a  deceleration 
to  hover  where  zero  forward  velocity  is  obtained.  During 
these  very  short  intervals  of  time,  the  pilot  is  commanding  a 
rotor  thrust  level  in  excess  of  that  required  to  hover.  This 
is  because  the  tip  path  plane  of  the  rotor  is  tilted  in  such  a 
way  as  to  either  command  an  increase  in  forward  airspeed  (a 
departure  from  hover  to  forward  flight)  or  a  decrease  in 
airspeed  (a  flare  to  a  hover  from  forward  flight).  A  third 
scenario  where  the  assumption  does  not  apply  involves  analyses 
of  hovering  tiltwing  aircraft  where  the  wing  is  tilted  either 


forward  or  aft  of  the  vertical.  In  these  situations,  the 
pitch  propeller  or  fan  is  used  to  trim  out  the  pitching  moment 
generated  by  the  tilted  wing  (In  contrast,  hovering  tiltrotor 
aircraft  always  have  the  plane  of  the  rotor  parallel  to  the 
ground  because  pitching  moments  can  not  be  trimmed  out.  If 
the  nacelle  angle  is  varied,  the  pitch  attitude  must  also 
vary ) . 

METHODOLOGY  IMPROVEMENTS 

Since  the  developed  methodology  was  incapable  of 
analyzing  scenarios  involving  tilted  rotor  tip  path  planes  in 
hover,  an  effort  was  initiated  in  1989  to  develop  a  solution 
for  this  problem.  Technical  aspects  of  this  effort  are 
documented  in  this  letter  report.  Modifications  to  the  user 
interface  of  the  rotorwash  analysis  computer  software  (ROTWASH 
program)  are  documented  in  reference  2. 

The  technical  approach  used  to  modify  the  reference  1 
mathematical  model  to  include  tilted  tip  path  plane  effects  is 
generally  based  upon  work  documented  in  reference  3.  This 
technical  approach  was  developed  us  ng  experimental  results 
which  are  referenced  in  the  memorandum.  Unfortunately,  no 
effort  was  made  to  document  any  correlation  of  results  with 
the  mentioned  experimental  data.  This  mathematical  modeling 
approach  is  graphically  described  using  figures  1  through  3. 

A  hovering  helicopter  at  a  given  gross  weight,  rotor  radius  of 
R  feet,  and  rotor  height  of  H  feet  (figure  1)  generates  a  mean 
induced  velocity  v^  or  dynamic  pressure  q*  across  the  rotor 
disk-  If  the  plane  of  the  rotor  is  parallel  to  the  ground, 
the  stream  tube  of  air  is  turned  90  degrees  upon  impact  with 
the  ground  and  a  circular  wall  jet  flowfield  forms  at  position 
Xp.  A  three-dimensional  view  of  this  process  is  presented  in 
figure  2.  As  the  horizontal  position  of  the  wall  jet  is 
increased  to  a  value  of  xm,  the  peak  velocity  Vjj,  in  the  wall 
jet  slowly  decays.  If  the  rotor  tip  path  plane  is  tilted 
slightly,  the  effect  depicted  in  figure  3  must  be  modeled.  In 
this  scenario,  assuming  the  center  of  the  rotor  is  in  the  same 
position,  the  impingement  point  of  the  stream  tube  moves 
further  aft  as  the  rotor  tip  path  plane  is  tilted  further 
forward.  This  shift  of  the  impingement  point  also  increases 
the  distance  xQ,  the  location  where  the  wall  jet  begins. 
Therefore,  for  any  fixed  distance  Xj^,  the  maximum  profile 
velocity  Vm  is  increased  when  the  tilt  angle  is  increased. 

Modifications  required  to  add  the  improvements  to  the 
reference  1  mathematical  model  are  easy  to  implement.  The 
user  first  defines  the  angle  that  the  tip  path  plane  will  be 
tilted.  This  angle,  0  ,  is  then  used  to  calculate  the  shift 
in  the  impingement  point  of  the  mean  value  of  the  rotor 
downwash.  The  corrected  length  to  the  impingement  point,  , 
is : 


=  H  /  cos  (p 


A -2 


The  horizontal  shift  in  the  location  of  the  point  is: 

Ax0  =  H  *  tan  0 

Analytically,  the  effect  of  the  parameter  H0  is  the  same  as  a 
change  in  the  altitude  of  the  rotor  above  ground  level  when 
the  tip  path  plane  is  not  tilted.  The  effect  of  the  parameter 
Ax0is  to  shift  the  location  of  the  formation  (and  strength)  of 
the  wall  jet  further  away  from  the  rotor. 

Two  practical  constraints  exist  with  respect  to  the  use 
of  the  tilted  tip  path  plane  mathematical  model  improvements . 
The  first  of  these  constraints  is  related  to  a  lack  of 
available  test  data  for  correlation  with  the  model.  At  some 
as  yet  undetermined  impingement  angle,  the  wall  jet  flowfield 
will  cease  to  form  as  a  circular  three-dimensional  flowfield. 
When  this  angle  is  exceeded,  the  mathematical  model  described 
reference  1  will  cease  to  be  valid.  Since  the  critical  angle 
is  presently  unknown,  it  is  recommended  that  tilt  angles 
greater  than  30  degrees  be  avoided  until  flight  test  data 
become  available  for  correlation.  This  recommendation  is  not 
based  on  any  particular  body  of  data,  only  experience  in  the 
analysis  of  rotorcraft. 

A  second  model  constraint  exists  due  to  the  fact  that  the 
model  is  not  designed  for  use  in  the  analysis  of  dynamic 
maneuvers  such  a3  a  takeoff.  The  model  can  only  be  used 
during  very  specific  segments  of  these  maneuvers  where 
airspeed  is  essentially  zero.  Steady  flow  assumptions 
inherent  to  the  model  are  "assumed"  to  be  valid  if  the  output 
data  are  viewed  as  snapshots  of  conditions  generated  during 
the  maneuver.  The  model  must  also  be  used  with  caution  since 
flight  test  data  are  not  yet  available  for  correlation.  The 
segment  of  the  takeoff  maneuver  which  can  be  analyzed  is  the 
initial  pushover  from  hover.  At  this  point  in  the  maneuver, 
airspeed  is  still  essentially  zero  and  the  tip  path  plane  is 
tilted  forward  to  accelerate  the  rotorcraft.  This  phase  of 
the  maneuver  also  requires  an  increase  in  rotor  thrust  above 
that  required  to  hover  in  order  to  avoid  ground  contact. 
Therefore,  the  user  must  know  something  about  the  gross  weight 
of  the  rotorcraft  and  the  thrust-power  relationships  involved 
to  make  practical  use  of  the  model  for  this  segment  of  the 
takeoff  maneuver.  During  a  landing,  the  final  phase  of  the 
flare  manuever  which  is  near  zero  airspeed  can  also  be 
analyzed.  During  this  phase  of  the  maneuver,  the  rotor  is 
tilted  back  very  close  to  its  maximum  position  and  the  thrust 
level  used  to  maintain  altitude  while  decelerating  is  at  its 
maximum. 

EXAMPLE  CASE 


The  effect  of  adding  the  rotor  tip  path  plane  feature  to 
the  reference  1  mathematical  model  is  demonstrated  in  this 
section  through  an  example  using  a  Bell  UH-1H  helicopter.  A 


A-3 


three-view  of  the  UH-1H  is  presented  in  figure  4.  The  gross 
weight  of  the  UH-1H  in  the  example  is  9500  pounds  and  the 
height  above  ground  level  of  the  rotor  is  40  feet. 

Results  from  calculations  for  hover  and  the  initiation  of 
a  takeoff  with  a  20  degree  forward  tilt  of  the  tip  path  plane 
are  presented  in  figure  5.  The  hover  rotorwash  profile 
characteristics  are  depicted  by  the  solid  lines  and  are 
provided  for  distances  of  50,  100,  and  150  feet  from  the 
center  of  the  rotor.  Peak  profile  velocities  vary  from 
approximately  48  knots  at  0.5  feet  above  the  ground  to  22 
knots  at  2  feet  for  the  50  and  150  feet  radial  positions 
respectively.  If  the  UH-1H  accelerates  from  hover  at  a  tip 
path  plane  of  20  degrees,  then  the  pilot  must  increase  thrust 
(through  the  collective  stick)  to  approximately  10500  pounds 
to  maintain  at  least  a  constant  altitude  above  the  ground. 
Calculated  results  at  the  instant  the  rotor  stabilizes  at  the 
20  degree  position  are  represented  by  the  dashed  lines.  Peak 
profile  velocities  at  the  50  and  150  feet  radial  positions  are 
increased  by  approximately  9  and  4  knots  respectively. 

CONCLUSIONS 


A  simple  model  has  been  developed  to  approximate  the 
effects  of  a  tilt  in  the  rotor  tip  path  plane  on  rotorwash 
characteristics  in  close  proximity  to  the  ground.  This  model 
has  been  integrated  into  the  ROTWASH  computer  program  and 
tested  with  example  cases.  Unfortunately,  at  this  time 
appropriate  flight  test  data  do  not  exist  which  can  be  used 
for  correlation  purposes.  These  data  will  have  to  be  acquired 
and  correlated  with  the  improved  mathematical  model  before  any 
further  improvements  can  be  made. 


Samuel  W.  Ferguson,  EMA 
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SINGLE  ROTOR  CONFIGURATION 


FIGURE  2  THREE-DIMENSIONAL  ROTORWASH  FLOW  FIELD 

CHARACTERISTICS  FOR  A  ROTOR  TIP  PATH  PLANE 
PARALLEL  TO  THE  GROUND 
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ROTORWASH  FLOW  FIELD  CHARACTERISTICS  FOR  A 
ROTOR  TIP  PATH  PLANE  WHICH  IS  TILTED  FORWARD 
WITH  RESPECT  TO  A  PLANE  PARALLEL  TO  THE  GROU 


7'-8  7/16" 


FIGURE  4  THREE-VIEW  OF  A  BELL  UH-1H  UTILITY  HELICOPTER 
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10“  2  1/2' 


BELL  UH-1H 

GROSS  WEIGHT  =  9500  POUNDS 
ROTOR  HEIGHT  =  40  FEET 
AIRSPEED  <  3  KNOTS 


133d  ‘33A31  QNnOHO  3A09V  1H9I3H 


FIGURE  5  UH-1H  ROTORWASH  CHARACTERISTICS  AT  RADIAL 

DISTANCES  OF  50,  100,  AND  150  FEET  FOR  TIP  PATH 
PLANE  ANGLES  OF  0  AND  20  DEGREES 
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PEAK  PROFILE  VELOCITY,  KNOTS 


TABLE  B-l  ROTORCRAFT  DATA  S 


APPENDIX  B 


ROTORCRAFT  DESIGN  DATA  IN  ROTWASH  FORMAT 
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W-30-300  15,500  21.8  2.0  10.38  1/4 

EH-101  31,500  30.5  5.0  10.78  1/5 


APPENDIX  C 


ROTWASH  PROGRAM  FORTRAN  77  LISTINGS 


ROTWASH  program  listings  are  presented  in  this 
appendix  for  the  ROTWASH  main  program  and  its  24 
subroutines.  The  listings  are  for  a  version  of  the  program 
that  is  run  on  IBM  PC  /  PC-compatible  computers  using 
MICROSOFT  FORTRAN  77,  Version  5.0.  The  tabular  listing 
below  indexes  subroutine  names  and  briefly  describes 
functionality  for  user  reference. 


SUBROUTINE 


FUNCTION 


ROTWASH 

CLOUD 

DEVTX 

FREAD 

GDVTX 


HAZARD 

HIPVEL 

HOMCLS 

HSVTX 

HWJVEL 


Main  Program  Driver  and  Initialization 
Calculates  Particle  Cloud  Boundaries 
Locates  Disk  Edge  Vortex  System  and 
Calculates  Induced  Velocity  Field 
Prompt/Validate  for  Floating  Point  Input  Data 
Locates  Ground  Vortex  System  and  Calculates 
Induced  Velocity  Field 

Driver  Subroutine  for  Hazard  Analysis 
Twin  Rotor  Overturning  Forces  and  Moments 
Home  Cursor  and  Clear  Screen 
Calculates  Induced  Velocity  Field  of  a 
Horseshoe  Vortex  System 

Single  Rotor  Overturning  Forces  and  Moments 


INKEY 

INPUT 

INPUTV 

INPUTX 

IOFNSH 


Menu  Input  Data  Control 

Rotorcraft  Characteristics  Input  Data  Menu 
Velocity  Profile  Status  Menu 
Ground/Disk  Vortex  Input  Data  Menu 
Close  Disk  I/O  Files 


IOINIT 

IPVEL 

IREAD 

LEGAL 

LOCATE 


File  I/O  Management  Menu 

Calculates  Interaction  Plane  Velocitv  Profile 
Prompt/Val idate  for  Integer  Input  Data 
Check  Validity  of  Input  Data  Selection  Codes 
Locate  Cursor  Position 


MOMENT 

PROPRM 

VLINE 

WALJET 

WJVEL 


Calculates  Personnel  Overturning  Forces  and 
Moments 

Calculates  Radial  Wall  Jet  Velocity  Profile 
Calculates  Induced  Velocity  from  a  Line 
Vortex  Field 

Defines  Initial  Wall  Jet  Position  and  Growth 
Parameters 

Calculates  Single  Rotor  Velocity  Profile 


PROGRAM  ROTWASH 


1  c 

2  C 

3 

4  C 

5  C 

6  C 

7  C 

8  C 

9  C 

10  C 

11  C 

12  C 

13  C 

14  C 

15  C 

16  C 

17  C 

18  C 

19  r 

20 

21  C 

22 

23 

24 

25 

26 

27 

28  C 

29 

30  C 

31 

32 

33 

34 

35 

36 

37 

38  C 

39 

40 

41 

42 

43 

44 

45 

46 

47  C 

48 

49  C 

50  C 

51  C 

52  C 

53  C 

54  C 


PROGRAM  ROTWASH 


ROTORCPAFT  DOWNWASH  HAZARD  ANALYSIS  PROGRAM 

EMA  /  COMPUTATIONAL  METHODOLOGY  ASSOCIATES 
SAMUEL  W.  FERGUSON  /  J.  DAVID  KOCUREK 

1  APRIL  1991 

PROGRAM  VERSION  2.0 

THIS  PROGRAM  WAS  DEVELOPED  USING  MICROSOFT  FORTRAN  77  FOR 
THE  DOS  OPERATING  SYSTEM.  CONSOLE  DISPLAY  CONTROL  IS 
PROVIDED  WITH  THIS  PROGRAM  IN  SEVERAL  SPECIAL  SUBROUTINES. 

*★★★'★*•*★****★★★*■*****'  ■*•*★**★■*★★★★•*■★*★■*•★★★★★★★★  k  ★*★★★★★★★★★★★ 


PARAMETER ( NUM  =  10) 

CHARACTER* 1  OKLIST(NUM) 

CHARACTER* 1  KE Y , KKEY , FLOW , VELHAZ , HAZTYP 
CHARACTER* 1  ICONT ( 5 ) 

CHARACTER* 1  TEMCHAR 
CHARACTER* 1 2  PTSFIL(4) 

CHARACTER* 50  COMM(2) 

DIMENSION  CONT ( 9 ) ,CONTV(7) ,CONTX(8) 

COMMON  /  CKEY/  KEY, KKEY 

COMMON  /CONSTS/  PI , RHO , FPSPKN , RHOD2 , DRC 

COMMON  /HELGEO/  H , DL , YSEP , WSPD , RADIUS , SHFTAN , DXO 

COMMON  /INPUTC/  ICONT , COMM , PTSFIL 

COMMON  /INPUTD/  CONT , CONTV , CONTX 

COMMON  /PROFIL/  RJ , ZB J , Z HJ , ZMJ , UMJ , ZB , ZH , ZM , UM , CU , CY 
COMMON  /  UNITS/  IOU1 , IOU4 , IOU5 , IOU6 , IOU7 , IOU8 , IGRAPH 

DATA  CONT/2. 0,32. 2, 12. 5, 130 00. 0,1 3. 0,37. 0,0. 0,1. 0,0.0/ 

DATA  ICONT/ 'V','N','N','W','L'/ 

DATA  CONTV/ 50 .0,1.0,10.0,50.0,10.0,100.0,0.0/ 

DATA  CONTX/7 33.0,7.0,18.0,200.0,0.0,20.0,100.0,0.0/ 

DATA  PTSFIL/ ' DFRC . PTS ' , 'DAIP.PTS ' , ' OTDFRC . PTS ' , 'OTDAIP.PTS ' / 
DATA  COMM/ 

1  'BELL  XV- 15  CHARACTERISTICS  ARE  ED  AS  INPUT  DATA  ’, 

2  'GROSS  WEIGHT  MIGHT  BE  ONE  OF  THE  OMMENT  STRINGS  '/ 

DATA  OKLIST  / ' W’ , V, 'I', ' i ' , 'G', ' q' , 'D' , 'd' , 'X', 'x  ’  / 


INITIALIZE  I/O  SYSTEM 
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55  C 

56  CALL  IOINIT 

57  C 

58  C  - 

59  C  INITIALIZE  MISCELLANEOUS  CONSTANTS 

60  C  - 

61  C 

62  PI  =  ACOS (-1.0) 

63  RHOSL  =  0.0023769 

64  FPSPKN  =  1.687 

65  DRC  =  0.01745329252 

66  C 

67  50  CONTINUE 

68  C 

69  KEY  =  *  ' 

70  KKEY=  '  1 

71  C 

72  C  - 

73  C  OBTAIN  INPUT  DATA  PARAMETERS  FROM  STATUS  SCREEN 

74  C  - 

75  C 

76  CALL  INPUT 


77 

c 

78 

c 

79 

c 

DEFINE 

PARAMETERS  OBTAINED  FROM  SUBROUTINE  INPUT 

80 

c 

81 

c 

82 

ROTORS 

=  CONT(l) 

83 

YYSEP 

=  CONT ( 2 ) 

84 

RADIUS 

=  CONT ( 3 ) 

85 

HELGW 

=  CONT ( 4 ) 

86 

DWNLD 

=  CONT ( 5 ) 

87 

HAGL 

=  CONT ( 6 ) 

88 

SHFTAN 

=  CONT ( 7 ) 

89 

SIGPR 

=  CONT ( 8 ) 

90 

WSPD 

=  CONT ( 9 ) 

91 

c 

92 

VELHAZ 

=  ICONT(l) 

93 

c 

94 

c 

95 

c 

ADJUST 

GEOMETRY  IF  SHAFT  ANGLE  >  0.0  DEGREES 

9  g 

c 

97 

c 

98 

RSHFT 

=  SHFTAN* DRC 

99 

CSHFTA 

=  COS ( RSHFT ) 

o 

o 

DXO 

=  HAGL*TAN( RSHFT) 

101  C 

102  C  - 

103  C  NON-DTMENSIONALIZE  SOME  OF  THE  INPUT  PARAMETERS 

104  C  - 

io:  C 

106  H  =  HAGL/R’IDIUS/CSHFTA 

107  YSEP  =  YYSEP/T . 0/RADIUS 

108  EFFGW  =  HELGW* (1.0  +  ( DWNLD/ 100.0) ) 
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109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 


DL  =  EFFGW/ ROTORS /PI /RADIUS*  *2 

RHO  =  SIGPR*RHOSL 

RHOD2  =  0 . 5*RHO 
C 

c  - 

C  SWITCHING  CALLS  HAZARD  PROGRAM  AND  ALLOWS  RETURN  (IF  DESIRED) 

C  TO  THE  MAINLINE  ROUTINE  TO  CHANGE  ROTORCRAFT  INPUT  PARAMETERS 

C  - 

c 

IF ( VELHAZ . EQ . ' H ' ) THEN 
C 

KKEY  =  ’H' 

CALL  HAZARD (HAZTYP) 

C 

IF ( KEY . EQ . ' X ' ) GOTO  999 
IF ( FLOW . EQ . ' X ' ) GOTO  999 
IF ( HAZTYP . EQ . ' X ' ) GOTO  999 
GOTO  50 
C 

ELSE 

C 

KKEY  =  'V' 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 

c 

c 

c 

c 


END  IF 


SELECT  FLOWFIELD  OPTION 


10  CONTINUE 


HOME  CURSOR  AND  CLEAR  SCREEN 


ICD  =  0 

CALL  HOMCLS ( ICD ) 
CALL  LOCATE (3,1) 


WRITE (IOU1 , 11 ) 


11 

FORMAT  (  2 OX,  SELECT  TYPE  OF  FLOW  TO 

BE  ESTIMATED',/// 

1 

2 OX,  'WALL  JET  PROFILE, 

TYPE 

<W> ' , /  , 

2 

2 OX,  ’ INTERACTION  PLANE  PROFILE, 

TYPE 

<!>',/, 

2 

20X, 'GROUND  VORTEX, 

TYPE 

<G> ' , / , 

2 

2 OX ,  ' DISK  VORTEX, 

TYPE 

<D> ' , / , 

3 

2 OX ,  'TO  EXIT  PROGRAM, 

TYPE 

<X>’ ,//) 

PROMPT  FOR,  OBTAIN,  AND  CHECK  FOR  LEGAL  INPUT 


40  CONTINUE 

WRITE (IOU1, 1 ( 23X,A, $ ) ' ) 
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163 

164  C 

165 

166  C 

167 

168  C 

169  C 

170  C 

171  C 

172  C 

173 

174 

175 

176 

177 

178  C 

179  C 

180  C 

181  C 

182  C 

183 

184 

185  C 

186  C 

187  C 

188  C 

189  C 

190  C 

191 

192 

193  C 

194 

195 
1S6 

197 

198 

199 

200 

201  C 

202 

203 

204 

205 

206 

207 

208 

209  C 

210 
211 
212 

213 

214 

215 

216 


1  '  ENTER  DATA  ENTRY  CODE  ==>  ’ 

READ ( IOU1 , ' ( A1 ) ' )  FLOW 

IF ( LEGAL ( FLOW, IOU1 , OKLIST , NUM ) ,EQ.l)GOTO  40 

MAKE  LEGAL  LOWERCASE  INPUTS  UPPER  CASE  BEFORE  BRANCHING 


IF ( FLOW . EQ . ' w 1 )  FLOW 
IF ( FLOW . EQ . 1 i ' )  FLOW 
IF ( FLOW . EQ .  '  g  '  )  FLOW 
IF ( FLOW . EQ . ' d ’ )  FLOW 
IF ( FLOW . EQ . ’ X ’ )  FLOW 


'W' 
'  I  ’ 
'G  ' 
'D' 
'  X' 


HOME  CURSOR  AND  CLEAR  SCREEN 


ICD  =  0 

CALL  HOMCLS ( ICD ) 


BRANCH  BASED  ON  CHOSEN  OPTION,  ALSO 

CHECK  NUMBER  OF  ROTORS  TO  LIMIT  SOME  OPTIONS 


I F ( FLOW . EQ . ' X ' ) GOTO  999 
I F ( FLOW . EQ . ' W ' ) GOTO  12 

IF ( FLOW . EQ . ' I ' ) THEN 

IF ( ROTORS. GT. 1.0 (GOTO  12 
WRITE (IOU1, ‘ (////, 14X,A,$)’( 

1  '  REQUIRES  TWO  ROTORS,  TYPE  <RETURN>  TO  CONTINUE  ’ 

READ( IOU1 , ' (A1 ) ' )  TEMCHAR 
GOTO  10 
ENDIF 

IF ( FLOW . EQ . 'G' (THEN 

IF (ROTORS.LT. 2.0 (GOTO  1000 
WRITE (IOU1, '(////, 14X, A, $)’) 

1  ’  REQUIRES  ONE  ROTOR,  TYPE  <RETURN>  TO  CONTINUE  ' 

READ ( I OU 1 , ' ( A 1 ) ' )  TEMCHAR 
GOTO  10 
ENDIF 

IF ( FLOW . EQ . ’D' (THEN 

IF (ROTORS.LT. 2.0 (GOTO  1000 
WRITE (IOU1, '(////, 14X, A, $)•) 

1  '  REQUIRES  ONE  ROTOR,  TYPE  <RETURN>  TO  CONTINUE  ' 

READ( IOU1 , ' (A1 ) ' )  TEMCHAR 
GOTO  10 
ENDIF 
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217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


GOTO  10 
12  CONTINUE 

★★★★★★★★★★★★★★★★★★★■ft* 

RADIAL  WALL  JET  FLOWS 

ACCELERATED  SLIPSTREAM  MEAN  VELOCITY 

UN  =  SQRT ( 2 . 0*DL/RHO) 

GROUND  EFFECT  CORRECTION 

AKG  =  1.0  -  0 . 9  *EXP ( -2 . 0*H ) 

MEAN  VELOCITY  AT  ROTOR  DISK  (RATIOED  TO  UN) 

UB  =  AKG/2.0 

FIND  INITIAL  RADIUS  OF  WALL  JET 

CALL  WALJET ( H , UB , UN , UMB ) 

500  CONTINUE 

IF(KEY.EQ. 'X' )GOTO  999 
IF(KEY.EQ. ’N ’ )GOTO  50 
IF ( FLOW . EQ . ’ I ' ) GOTO  700 

600  CONTINUE 

WALL  JET  REGION 

OBTAIN  INPUT  DATA  FOR  THE  WALL  JET  OPTION 

CALL  INPUTV ( FLOW ) 

RVZ  =  (CONTV(l)  -  DXO) /RADIUS 
DELZ  =  CONTV( 2) /RADIUS 
ZMAX  =  CONTVf 3) /RADIUS 


C 
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271  C  GENERATE  VELOCITY  PROFILE  AT  RVZ  IN  WALL  JET  REGION 

272  C  - 

273  C 

274  CALL  WJVEL ( H , UN , UMB , RVZ , RADIUS , WSPD , DELZ , ZMAX , DXO) 

275  C 

276  GOTO  500 

277  C 

278  700  CONTINUE 

279  C 

280  C  - 

281  C  INTERACTION  PLANE  UPWASH  DEFLECTION  ZONE 

282  C 

283  C  OBTAIN  INPUT  DATA  FOR  THE  I PLANE  OPTION 

284  C  - 

285  C 

286  CALL  INPUTV ( FLOW ) 

287  C 

288  XIP  =  (CONTV(l)  -  DXO) /RADIUS 

289  DELZ  =  CONTV( 2 ) /RADIUS 

290  ZMAX  =  CONTV( 3 ) /RADIUS 

291  C 

292  C  - 

293  C  GENERATE  VELOCITY  PROFILE  AT  XIP  IN  INTERACTION  PLANE 

294  C  - 

295  C 

296  CALL  IPVEL( H, UN, RADIUS, UMB, XIP, YSEP , WSPD , DELZ , ZMAX , DXO ) 

297  C 

298  GOTO  500 

299  C 

300  1100  CONTINUE 

301  C 

302  IF(KEY.EQ. 'X' )GOTO  999 

303  IF ( KEY . EQ . ' N ' ) GOTO  50 

304  C 

305  1000  CONTINUE 

306  C 

307  C  ********************** 

308  C  HORSESHOE  VORTEX  FLOWS 

309  c  ********************** 


310  C 

311  C 


312 

313 

c 

Q 

OBTAIN 

INPUT  DATA  FOR  VORTEX  OPTIONS 

314 

c 

315 

316 

c 

CALL  I 

NPUTX 

317 

OMEGAR 

= 

CONTX ( 1 ) 

318 

B 

= 

CONTX ( 2 ) 

319 

VF 

= 

CONTX ( 3 ) 

320 

XT 

= 

CONTX (4 ) /RADIUS 

321 

YT 

= 

CONTX ( 5 ) /RADIUS 

322 

DELZ 

= 

CONTX ( 6 ) /RADIUS 

323 

324 

c 

ZMAX 

CONTX ( " ) /RADIUS 
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325 

VF  =  VF*FPSPKN 

326 

AMU  =  VF/OMEGAR 

327 

CT  =  DL/RHO/ OMEGAR*  *  2 

328 

C 

329 

C 

330 

C 

ITERATE  TO  GET  INFLOW  RATIO 

331 

C 

332 

C 

333 

ALOLD  =  SQRT(CT/2 . 0) 

334 

C 

335 

DO  1300  ITER=1 / 100 

336 

C 

337 

ALNEW  =  CT/2. 0/SORT (ALOLD**2  +  AMU**2) 

338 

C 

339 

IF(ABS(ALNEW  -  ALOLD) .LE . 1 . 0E-05 }GOTO  1301 

340 

C 

341 

ALOLD  =  ALNEW 

342 

C 

343 

1300 

CONTINUE 

344 

c 

345 

c 

346 

c 

HOME  CURSOR  AND  CLEAR  SCREEN 

347 

c 

348 

c 

349 

ICD  =  0 

350 

CALL  HOMCLS ( ICD ) 

351 

CALL  LOCATE (5,1) 

352 

c 

35. 

WRITE ( IOU1 ,20) 

354 

20 

FORMAT (  '************************************’, 

355 

1  , ' ITERATIONS  EXCEEDED  FOR  INFLOW  RATIO 

355 

2  , ************************************* 

357 

C 

35b 

STOP  ’  ' 

35 

c 

3  6  ‘  • 

1301 

CONTINUE 

36 

C 

36; 

ALAMDA  =  ALNEW 

36  ■ 

AMUS  =  AMU /SQRT (CT/2 . 0 ) 

36  1 

GAMT  =  OM£GAR*RADIUS*2 . Q*PI*CT/B 

36  i 

CHI  =  AT AN (ALAMDA /AMU)/ 2.0 

36 ; 

GAMW  «  PI*RADIUS*OMEGAR**2*CT/VF/2 . 0 

36 

HOD  =  H/2.0 

365 

C 

365 

c 

37  C 

c 

HOME  CURSOR  AND  CLEAR  SCREEN 

371 

c 

372 

c 

373 

ICD  =  0 

374 

CALL  HOMCLS (ICD) 

375 

CALL  LOCATE (3,1) 

376 

c 

377 

I F ( FLOW . EO • ' D ’ ) GOTO  1200 

378 

c 
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379  C  - 

380  C  GROUND  VORTEX 

381  C  - 

382  C 

383  WRITE( IOU1 , 1001 )  HOD,AMUS,AMU 

384  1001  FORMAT (  18X, 'ROTOR  HEIGHT  ABOVE  GROUND  H/D  ',2X,F8.4,/ 

385  1  , 18X, 'ADVANCE  RATIO  MU-STAR  ',2X,F8.4,/ 

386  2  , 18X, 'ADVANCE  RATIO  MU  ' , 2X , F8 . 4 , / / ) 

387  C 

388  C  - 

389  C  THE  VALUE  INPUT  HERE  REQUIRES  USE  OF  THE  CHART  IN  FIGURE  18 

390  C  OF  THE  ACCOMPANYING  DOCUMENTATION  FOR  THE  GROUND  VORTEX 

391  c  - 

392  C 

393  30  CONTINUE 

394  C 

395  WRITE (IOU1, 31) 

396  31  FORMAT (  18X, 'ENTER  GROUND  VORTEX  STRENGTH  RATIO ',/ , 

397  1  25X, ' (SEE  FIGURE  18)  ==>  ’$) 

398  C 

399  READ ( IOU1 , * , ERR=3  0 )  GAMG 

400  C 

401  IF ( GAMG. LT. 0.0)  GAMG  =  0.0 

402  C 

403  C  - 

404  C  CONTINUE  ANALYSIS 

405  C  - 

406  C 

407  GAMG  =  GAMG*GAMT 

408  C 

409  CALL  GDVTX ( H , RADIUS , AMU ,CT ,GAMG , XT , YT , DELZ , ZMAX ) 

410  C 

411  GOTO  1100 

412  C 

413  1200  CONTINUE 

414  C 

415  C  - 

416  C  DISK  VORTEX 

417  C  - 

418  C 

419  CALL  DEVTX ( H , RADIUS , GAMW, CHI , XT , YT , DELZ , ZMAX ) 

420  C 

421  GOTO  1100 

422  C 

423  C  - 

424  C  NORMAL  PROGRAM  EXIT 

425  C  - 

426  C 

427  999  CONTINUE 

428  C 

429  END 

430  C 

431 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 


C 

C 

SUBROUTINE  CLOUD ( UN , UMB ) 

C 

Q  *★*★*★★■*★*★★*★*★*★*★***★*★★*★★★★★*★★★★★★★★★★*★**★★***★★*★★★*★★ 

C  SUBROUTINE  CLOUD 

C 

C  THIS  SUBROUTINE  MAKES  THE  CALCULATIONS  REQUIRED  IN  ESTIMATING 

C  THE  PARTICLE  CLOUD  BOUNDARIES  (  NO  DENSITIES  )  FOR  SINGLE  AND 

C  TWIN  ROTOR  CONFIGURATIONS 

c 

CHARACTER* 1  KEY,KKEY 
CHARACTER* 1  ICONT ( 5 ) 

CHARACTER* 12  PTSFIL(4) 

CHARACTER* 50  COMM(2) 

C 

COMMON  /  CKEY/  KEY , KKEY 
COMMON  / CLOUDK/  QSMAX 

COMMON  /CONSTS/  PI , RHO, FPSPKN , RHOD2 , DRC 

COMMON  /HELGEO/  H , DL , YSEP, WSPD , RADIUS , SHFTAN, DXO 

COMMON  /INPUTC/  ICONT , COMM, PTSFIL 

COMMON  /PROFIL/  RJ , 2BJ , ZH J , ZMJ , UMJ , ZB , ZH , ZM, UM, CU , CY 
COMMON  /  UNITS/  IOU1 , IOU4 , IOU5 , IOU6 , IOU7 , IOU8 , IGRAPH 
C 

C  if************************************************************* 

c 

c  - 

C  CLEAR  SCREEN  AND  HOME  CURSOR 

C  - 

C 

ICD  =  0 

CALL  HOMCLS ( ICD ) 

CALL  LOCATE (3,1) 

C 

C  - 

C  READ  IN  THE  TERRAIN  FACTOR  CONSTANT 

C  (SEE  FIGURE  30  OF  USER'S  GUIDE) 

C  - 

C 

10  CONTINUE 
C 

WRITE (IOU1, 20) 

20  FORMAT (  15X,'ENTER  TERRAIN  EROSION  FACTOR  ( -ND- )  ==>  ’,$) 

C 

READ ( IOU1 , * , ERR=1 0 )  XKT 
C 

C  - 

C  VALIDATE  REAL  INPUT  VALUE 

C  - 

C 

I F ( XKT . LE . 0 . 0 . OR . XKT . GT . 5  0  0 . 0 ) GOTO  1 0 
C 

c  - 
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55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 


C  DEFINE  CLOUD  BOUNDARY  CONSTANTS 

C  - 

c 

XKT  =  SQRT ( XKT ) 

C 

QSMX  =  RHOD2  * ( ( SQRT ( QSMAX ) *  UN ) *  *  2 ) 

ERC  =  -0.437 
XUM  =  ( UMB  *  UN ) *  *  2 
XCU  =  CU*CU 
Cl  =1.0 
C2  =2.2 
C 

C  - 

C  SINGLE  ROTOR  CLOUD  BOUNDARY  CALCULATIONS 

C  - 

c 

RCR  =  RADIUS* ( ( XKT/ (Cl  * RHOD2 *XUM*XCU ) ) **ERC ) 

RVR  =  0 . 785*RCR 

ZVR  =  0 . 329  *RCR 

RCVR  =  RCR  -  RVR 

AR  =  ( 2 . 0 /PI ) *ALOG ( ZVR /RCVR ) 

PHIR  =  ( PI /2 . 0 ) *ALOG ( RCVR) / ALOG ( ZVR/RCVR) 

AXLV  =  AR*( (-PI/2.0)  +  PHIR) 

XLV  =  EXP (AXLV) 

HCR  =  XLV  +  ZVR 
C 

c  - 

C  INITIALIZE  INTERACTION  PLANE  BOUNDARIES 

C  - 

c 

RCI  =0.0 
RVI  =  0.0 
ZVI  =0.0 
HCI  =  0.0 

IF ( YSEP . LE .0.1) GOTO  30 
C 

c  - 

C  INTERACTION  PLANE  CLOUD  BOUNDARY  CALCULATIONS 

C  - 

C 

RCI  =  RADIUS* (  ( XKT/ ( C2  *RHOD2*XUM*XCU) ) **ERC ) 

RVI  =  0 . 7  8  5  *  RC I 

ZVI  =  0 . 3  2  9  *  RC I 

RCVR  =  RCI  -  RVI 

AR  =  (2 . 0/PI ) *ALOG( ZVI/RCVR) 

PHIR  =  ( PI /2 . 0 ) *ALOG ( RCVR ) / ALOG ( ZVI/RCVR ) 

AXLV  =  AR* ( ( -PI / 2.0)  +  PHIR) 

XLV  =  EXP (AXLV) 

HCI  =  XLV  +  ZVI 
C 

30  CONTINUE 
C 

c  - 

C  PRINTOUT  OF  BOUNDARY  LIMITS 
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109  C  - 

110  C 

111  IF ( IOU6 . NE . IOU1 )  WRITE ( IOU6 , ' ( ' ' 1 ' ' ) ‘ ) 

112  C 

113  C  - 

114  C  WRITE  "40  FORMAT"  IF  OUTPUT  TO  GRAPHICS  FILE 

115  C  - 

116  C 

117  IF  < IOU6 . EQ. 6 )  WRITE ( IOU6 , 4 0 )  COMM( 1 ) , COMM( 2 ) 

118  40  FORMAT (  10X, A50 , / , 10X, A50 , // ) 

119  C 

120  WRITE (IOU6, 50) 

121  50  FORMAT (  //, 

122  1  20X,  SUMMARY  OF  CLOUD  BOUNDARIES',//, 

123  2  20X,  RC  AND  RV  ARE  FROM  ROTOR  CENTER  (FT)’,//, 

124  3  20X,  RC  '  , 7X,  RV  ’,7X,'  ZV  ’,7X,’  HC  ',/) 

125  C 

126  WRITE (IOU6, 60)  RCR, RVR, ZVR , HCR 

127  60  FORMAT (  1 3X , ' SR ' , F10 . 1 , 3F1 1 . 1 ) 

128  C 

129  WRITE (IOU6, 70)  RCI , RVI , 2VI , HCI 

130  70  FORMAT (  1 3X , ' IP ' , FI 0 . 1 , 3F1 1 . 1 ) 

131  C 

132  WRITE (IOU6, 80)  QSMX 

133  80  FORMAT (  /,12X,'  QSMAX  =',F7.1,’  PSF',//) 

134  C 

135  C  - 

136  C  DECIDE  NEXT  OPTION  WITH  INKEY 

137  C  - 

138  C 

139  CALL  INKEY 

140  C 

141  RETURN 

142  END 

143  C 

144 
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1  c 

2  C 

3 

4  C 

5  C 

6  C 

7  C 

8  C 

9  C 

10  c 

11  c 

12 

13 

14 

15 

16 

17  C 

18 

19 

20 
21 
22 

23 

24  C 

25  C 

26  C 

27  C 

28  C 

29  C 

30  C 

31 

32 

33 

34  C 

35 

36 

37 

38  C 

39  C 

40  C 
4  j.  C 

42  C 

43  C 

44 

45 

46 

47  C 

48 

49 

50 

51  C 

52  C 

53  C 

54  C 


SUBROUTINE  DEVTX ( H , RADIUS , GAMW , CHI , XT , YT , DELZ , ZMAX ) 

SUBROUTINE  DEVTX 

THIS  SUBROUTINE  LOCATES  THE  DISK  EDGE  VORTEX  SYSTEM,  AND 
DIRECTS  THE  CALCULATION  OF  ITS  INDUCED  VELOCITY  FIELD 


CHARACTER* 1  TEMCHAR 
CHARACTER* 1  KEY , KKEY 
CHARACTER* 1  ICONT{5) 

CHARACTER* 12  PTSFIL(4) 

CHARACTER* 50  COMM ( 2 ) 

COMMON  /  CKEY /  KEY, KKEY 

COMMON  /CHSVTX/  XL1 , YL1 , ZL1 , XL2 , YL2 , ZL2 , XL3 , YL3 , ZL3 , 
1  XR1 , YR1 , ZR1 , XR2 , YR2 , ZR2 , XR3 , YR3 , ZR3 

COMMON  /CONSTS/  PI , RHO , FPSPKN , RHOD2 , DRC 
COMMON  /INPUTC/  ICONT , COMM, PTSFIL 

COMMON  /  UNITS/  IOU1 , IOU4 , IOU5 , IOU6 , IOU7 , IOU8 , IGRAPH 


ASSUME  HORSESHOE  SHAPE  -  ASSIGN  LEFT  AND  RIGHT  CORNERS 


XL1  =0.0 
YL1  =  -1.0 
ZL1  =  H 

XR1  =0.0 
YR1  =  1.0 
ZR1  =  H 


SET  UP  DIRECTION  POINTERS  FOR  TRAILER  ELEMENTS 
POINT  2  IS  AT  GROUND  IMPINGEMENT 


XL 2  =  XL1  +  H/TAN (CHI ) 

YL2  =  YL1 

ZL2  =  0.0 

XR2  =  XR1  +  H/TAN (CHI ) 

YP2  =  YR1 

ZR2  =  0.0 


POINT  THREE  EXTENDS  TRAILER  PARALLEL  TO  GROUND 


DEVTX- 1 


SUBROUTINE  DEVTX 


55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 


C 


XL  3 

= 

XL2  + 

1.0 

YL3 

= 

YL2 

ZL 

= 

ZL2 

XR3 

SB 

XR2  + 

1.0 

YR3 

= 

YR2 

ZR3 

= 

ZR2 

c 

XT  =  XT*RADIUS 
YT  =  YT*RADIUS 
C 

c  - 

C  CLEAR  SCREEN  AND  HOME  CURSOR 

C  - 

c 

ICD  =  0 

CALL  HOMCLS ( ICD ) 

CALL  LOCATE (4,1) 

C 

C  - 

C  WRITE  OUTPUT  HEADER 

C  - 

c 

IF ( IOU6 . NE . IOU1 )  WRITE ( IOU6 , ' ( ' ' 1 ' ' ) ' ) 

C 

WRITE ( IOU6 , 1000  )  XT ,  YT 

1000  FORMAT (  21X,'DISK  VORTEX  VELOCITY  PROFILE  DATA', 

1  /// , 14X,  ‘  X-LOCATION  (XT)  =  ' ,2X,F8.2,2X,  ' FT  1 , 

2  / , 1 4X ,  Y -LOCATION  (YT)  =  ,2X,F8.2,2X, ' FT ' , / ) 

C 

WRITE (IOU6, 1001)  GAMW 

1001  FORMAT (  15X,  'VORTEX  CIRCULATION  =  '  , FI 0 . 2 , 2X ,  ' FT* *2 /SEC ' ) 
C 

CHID  =  CHI* 180.0/PI 
C 

WRITE(IOU6, 1005)  CHID 

1005  FORMAT (  15X,  ‘SETTLING  ANGLE  =  '  , FI  0 . 2 , 2X ,  ' DEG ’  , //// ) 

C 

WRITE (IOU1, ' (23X,A,$) ' ) 

1  '  PRESS  <RETURN>  TO  CONTINUE 

C 

READ (IOU1, ' (Al) ' )  TEMCHAR 
C 

CALL  HOMCLS (ICD) 

C 

IF ( IOU6 . EQ . 6 )  WRITE ( IOU6 ,93 )  COMM ( 1 ) , COMM ( 2 ) 

93  FORMAT (  14X,A50,/,10X,A50, //I 
C 

WRITE ( IOU6 ,1100) 

1100  FORMAT (  12X, 'HEIGHT ', 8X, 'MEAN  VELOCITY ', 5X ,' MEAN  Q',/, 

1  ' 0 ' , 12X, ' (FT) ' , 8X, ' (FPS)  , 6X , ’ ( KN ) ' , 5X, ' (PSF) ' , / ) 

C 

c  - 
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109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 
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135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 


C  SET  UP  SWEEP  OF  Z  AT  SPECIFIED  X,Y 

C  DELZ  AND  ZMAX  COME  FROM  A  MAINLINE  STATUS  MENU 

C  - 

C 

XT  =  XT/RADIUS 

YT  =  YT/RADIUS 

NPTS  =  IFIX ( ZMAX /DELZ )  +  2 
LINES  =  0 
C 

DO  200  1=1, NPTS 
C 

LINES  =  LINES  +  4 

ZT  =  (I  -  1 ) *DELZ 

C 

CALL  HSVTX(XT, YT, ZT,VXF,VYF, VZF,GAMW, RADIUS) 

C 

ZZ  =  ZT*RADIUS 

VTF  =  SQRT ( VXF*  *  2  +  VYF**2  +  VZF**2) 

VXK  =  VXF / FPSPKN 
VYK  =  VYF/FPSPKN 
VZK  =  VZF/FPSPKN 
VTK  =  VTF/FPSPKN 
C 

QX  =  RHOD2  *VXF*  *  2 
QY  =  RHOD2*VYF**2 
QZ  =  RHOD2*VZF**2 
QT  =  RHOD2*VTF**2 

c 

c  - 

C  KEEP  OUTPUT  PAGE  LENGTH  TO  SIZE  OF  SCREEN 

C  - 

C 

IF ( IOU6 . EQ . IOU 1 ) THEN 

IF ( LINES . LE . 1 2 ) GOTO  100 
LINES  =  4 
CALL  INKEY 

IF ( KEY . ME . ' C ' ) GOTO  999 
WRITE (IOU6, 1100) 

ENDIF 

C 

100  CONTINUE 
C 

c  - 

C  REPORT  X  COMPONENT  OF  VELOCITY 

C  _ 

C 

WRITE (IOU6, 1101 )  2  Z , VXF , VXK , OX 
1101  FORMAT (  9X,F8.2,2X,  X' , 3F 1 0 . 3 ) 

C 

c  - 

C  REPORT  Y  COMPONENT  OF  VELOCITY 

C  - 

C 

WRITE ( IOU6 ,1102)  VYF , VYK , QY 
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163 

164 

165 

1102 

C 

c 

FORMAT (  1 9  X , ‘ Y ' , 3 F 1 0 . 3 ) 

166 

C 

REPORT  Z  COMPONENT  OF  VELOCITY 

167 

c 

168 

c 

169 

WRITE ( IOU6 , 1103 )  VZF,VZK,QZ 

170 

1103 

FORMAT (  1 9  X ,  ' Z '  , 3F 1 0 . 3 ) 

171 

c 

172 

173 

c 

c 

REPORT  TOTAL  VELOCITY 

H4 

c 

175 

c 

176 

WRITE ( IOU6 , 1 104 )  VTF , VTK , QT 

177 

1104 

FORMAT (  19X, ' T ‘ , 3F1 0.3) 

178 

c 

179 

200 

CONTINUE 

180 

C 

181 

CALL  INKEY 

182 

c 

183 

999 

CONTINUE 

184 

C 

185 

RETURN 

186 

END 

187 

188 

C 

DSVTX-4 


SUBROUTINE  FREAD 


1  c 

2  C 

3  SUBROUTINE  FREAD ( IGU1 , PROMPT , VALUE , CONST ) 

4  C 

^  £  **★*★★***★*★★★*★★★★****★*★★★★★★***★***★★★★★■*★★ 

6  C  SUBROUTINE  FREAD  PROMPTS  USER  FOR  A  FLOATING 

7  C  POINT  DATA  ENTRY  AND  CHECKS  VALIDITY  OF  ENTRY 

g  c  ********************************************** 

9  C 

10  PARAMETER (LAST=50) 

11  C 

12  CHARACTER* 50  PROMPT , SHOWIT 

13  CHARACTER* 15  ENTRY, BLANK 

14  C 

15  DATA  BLANK  / ’  ' < 

16  C 

17  C  ******  +  **********************************■-**** 

18  C 

19  c  - 

20  C  PROMPT  USER  FOR  SCALED  FLOATING  POINT  ENTRY. 

21  C  FIND  POSITION  OF  LAST  NON-BLANK  CHARACTER  IN  PROMPT, 

22  C  THEN  STORE  RIGHT  JUSTIFIED  IN  SHOWIT. 

23  C  - 

24  C 

25  N  =  LAST  -r  1 

26  C 

27  10  IF ( N . EQ . 1 ) GOTO  20 

28  C 

29  N  =  N  -  1 

30  C 

31  IF ( PROMPT (N: N) .EQ.  ')GOTO  10 

32  C 

33  20  JS  =  LAST  -  N 

34  C 

35  WRITE ( SHOWIT, '( 50A1 )’ )  ('  ' , J= 1 , J S ) , ( PROMPT ( I ; I ) , 1=1 , N ) 

36  C 

37  c  - 

38  C  NOW  ASK  USER  FOR  DATA  ENTRY 

39  c  - 

40  C 

41  30  WRITE ( IOU 1, ’ (/,1X,A,G13.6) ' )  SHOWIT , VALUE*CONST 

42  C 

43  WRITE ( IOU1 , ' ( / , 8X , A , $ ) ') 

44  1  '  ENTER  NEW  VALUE  OR  <RETURN>  TO  LEAVE  AS  IS  =  -->  ' 

45  C 

46  READ ( IOU1 , ' ( A) ' )  ENTRY 

47  C 

48  IF ( ENTRY. EQ. BLANK) RETURN 

49  C 

50  READ (ENTRY, ' (BN, F15.0) ' ,ERR=30)  TEMP 

51  C 

52  C  - 

53  C  CONSTANT  CAN  BE  USED  TO  SCALE  OR 

54  C  OR  CONVERT  UNITS  OF  AH  INPUT  VALUE 


FREAD- 1 


SUBROUTINE  FREAD 


VALUE  =  TEMP/CONST 

RETURN 

END 


FREAD-2 


SUBROUTINE  GDVTX 


1  c 

2  C 

3 

4  C 

5  C 

6  C 

7  C 

8  C 

9  C 

10  c 

11  c 

12  C 

13  C 

14  C 

15  C 

16 

17 

18 

19 

20 

21  C 

22 

23 

24 

25 

26 

27 

28  C 

29  C 

30  C 

31 

32 

33 

34 

35 

36 

37 

38  C 

39  C 

40  C 

41  C 

42  C 

43 

44 

45 

46  C 

47 

48 

49 

50  C 

51  C 

52  C 

53  C 

54  C 


SUBROUTINE  GDVTX ( H , RADIUS , AMU , CT , GAMG , XT , YT , DELZ , ZMAX) 

SUBROUTINE  GDVTX 

THIS  SUBROUTINE  LOCATES  THE  GROUND  VORTEX  BASED  ON  THE 
EXPERIMENTS  BY  SUN  AND  CURTIS  (PRINCETON  UNIV. ) ,  AND  THEN 
DIRECTS  THE  CALCULATION  OF  ITS  INDUCED  VELOCITY  FIELD 

THE  OUTPUT  FROM  THIS  SUBROUTINE  SHOULD  BE  USED  CAREFULLY 
FOR  GROSS  ESTIMATION  PURPOSES  ONLY 

CHARACTER* 1  ICONT ( 5 ) 

CHARACTER* 1  TEMCHAR 
CHARACTER* 1  KEY , KKEY 
CHARACTER* 12  PTSFIL(4) 

CHARACTER*50  COMM(2) 

COMMON  /  CKEY /  KEY, KKEY 

COMMON  / CHSVTX/  XL1 , YL1 , ZL1 , XL2 , YL2 , ZL2 , XL3 , YL3 , ZL3 , 

1  XR1 , YR1 , ZR1 , XR2 , YR2 , ZR2 , XR3 , YR3 , ZR3 

COMMON  /CONSTS/  PI , RHO , FPSPKN , RHOD2 , DRC 
COMMON  / INPUTC/  ICONT , COMM , PTSFIL 

COMMON  /  UNITS/  IOU1 , IOU4 , IOU5 , IOU6 , IOU7 , IOU8 , IGRAPH 
★  *****★★*★*★*★★★*★★*★★**★*★★★★***★★*★★★★■**★★★★★★*★*★★★★★★★ 


HOD 

= 

H/2.0 

Cl 

= 

1.0  +  1 . 2086*HOD**0 . 43 

C2 

= 

-0.2786*  HOD  **0.6757 

ZGV 

= 

-10. 0*AMU  +  0.6 

XGV 

= 

-(Cl  +  C2  * ( AMU/CT ) ) **2 

XXGV 

= 

XGV* RAD I US 

ZZGV 

ZGV*RADIUS 

ASSUME  HORSESHOE  SHAPE  -  ASSIGN  LEFT  AND  RIGHT  CORNERS 


XL 1  =  XGV 
YL1  =  -1.0 
ZL1  =  ZGV 

XRl  =  XGV 
YR1  =1.0 
ZR1  =  ZGV 


SET  UP  DIRECTION  POINTERS  FOR  TRAILER  ELEMENTS 


GDVTX- 1 


SUBROUTINE  GDVTX 


55  XL 2  =  XL 1  +1.0 

56  YL2  =  YL1 

57  ZL2  =  ZL1 

58  C 

59  XR2  =  XR1  +1.0 

60  YR2  =  YR1 

61  ZR2  =  ZR1 

62  C 

63  XL 3  =  XL 2  +1.0 

64  YL3  =  YL2 

65  ZL3  =  ZL2 

66  C 

67  XR3  -  XR2  +1.0 

68  YR3  =  YR2 

69  ZR3  =  ZR2 

70  C 

71  C  - 

72  C  WRITE  OUT  GROUND  VORTEX  POSITION/STRENGTH  DATA 

73  c  - 

74  C 

75  IF ( IOU6 . NE . IOU1 )  WRITE ( IOU6 1 ) 

76  C 

77  WRITE (IOU6, 1001 )  XXGV,ZZGV 

78  1001  FORMAT (  //, 24X GROUND  VORTEX  CORE  POSITION',//, 

79  1  18X, 'X-LOCATION  (XXGV)  =  ' , IX, F8 . 2 , 2X, ’ FT  * , / 

80  2  18X, ' Y-LOCATION  (ZZGV)  =  ' , IX, F8 . 2 , 2X, ' FT ' ) 

81  C 

82  WRITE  (T.OU6, 1002)  GAMG 

83  1002  FORMAT (  /, 1 8X ,' GROUND  VORTEX  CIRCULATION  =  ’ , IX, 

84  1  F8 . 2 , 2X , ' FT* *2 / SEC ' , / / ) 

85  C 

86  WRITE(IOUl, ’ (23X,A,$) ' ) 

87  1  ’  PRESS  <RETURN>  TO  CONTINUE  ' 

88  C 

89  READ ( IOU1 , ' ( A1 ) ' )  TEMCHAR 

90  C 

91  ICD  =  0 

92  CALL  HOMCLS ( ICD ) 

93  C 

94  C  - 

95  C  WRITE  OUTPUT  HEADER 

96  C  - 

97  C 

98  IF ( IOU6 . EQ . 6 )  WRITE ( IOU6 , 9 3 )  COMM ( 1 ) , COMM ( 2 ) 

99  93  FORMAT (  1 4X , A5 0 , / , 1  OX , A50 , / / ) 

100  C 

101  WRITE (IOU6, 1100) 

102  1100  FORMAT {  1 2 X HEIGHT 8X MEAN  VELOCITY 5X MEAN  Q',/, 

103  1  ' 0 ' , 1 2X , ' ( FT ) ' , 8X , ' ( FPS ) ' , 6X , ' ( KN ) ' , 5X , ’ ( PSF ) ' , / ) 

104  C 

105  C  - 

106  C  SET  UP  SWEEP  OF  Z  AT  SPECIFIED  X,Y 

107  C  DELZ  AND  ZMAX  COME  FROM  MAINLINE  STATUS  SCREEN 

108  C  - 
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109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 


C 

NPTS  =  IF1X ( ZMAX/DELZ )  +  2 
LINES  =  0 
C 

DO  200  1=1, NPTS 
C 

LINES  =  LINES  +  4 
ZT  =  (I  -  1 ) *DELZ 
C 

CALL  HSVTX ( XT , YT , Z  T , VXF , V YF , VZ  F , GAMG , RADIUS ) 
C 

ZZ  =  ZT*RADIUS 

VTF  =  SQRT ( VXF**2  +  VYF**2  +  VZF**2) 

VXK  =  VXF/FPSPKN 
VYK  =  VYF/FPSPKN 
VZK  =  VZF/FPSPKN 
VTK  =  VTF/FPSPKN 
C 

QX  =  RHOD2*VXF**2 
QY  =  RHOD2*VYF**2 
QZ  =  RHOD2*VZF**2 
QT  =  RHOD2*VTF**2 
C 

IF  (IOU6.EQ.IOUl)THEN 
IF ( LINES . LE .12) GOTO  100 
LINES  =  4 
CALL  INKEY 

IF (KEY.NE . 'C ' )GOTO  999 
WRITE (IOU6, 1100) 

ENDIF 

C 

100  CONTINUE 
C 

C  - 

C  REPORT  X  COMPONENT  OF  VELOCITY 

C  - 

C 

WRITE (IOU6, 1101)  Z  Z , VXF , VXK , QX 

1101  FORMAT (  9X,F8.2,2X, 'X' , 3F1 0 . 3 ) 

C 

c  - 

C  REPORT  Y  COMPONENT  OF  VELOCITY 

C  - - 

C 

WRITE ( IOU6 ,1102)  VYF,VYK,QY 

1102  FORMAT (  19X ,  ' Y '  , 3F 1 0 . 3 ) 

C 

c  - 

C  REPORT  Z  COMPONENT  OF  VELOCITY 

C  - 

C 

WRITE ( IOU6 ,1103)  VZF,VZK,QZ 

1103  FORMAT (  1 9X ,  ' Z  '  , 3F1 0 . 3 ) 

C 


GDVTX-3 
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163  C  - 

164  C  REPORT  TOTAL  VELOCITY 

165  C  - 

166  C 

167  WRITE ( IOU6 ,1104)  VTF, VTK,QT 

168  1104  FORMAT (  19X, ' T ' , 3F10 . 3 ) 

169  C 

170  200  CONTINUE 

171  C 

172  CALL  INKEY 

173  C 

174  999  CONTINUE 

175  C 

176  RETURN 

177  END 

178  C 

179 
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SUBROUTINE  HAZARD 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 
4  4 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 


C 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 

c 

c 


c 

c 


SUBROUTINE  HAZARD ( HAZTYP ) 

SUBROUTINE  HAZARD  IS  THE  MAINLINE  DRIVER  FOR  THE 
CALCULATION  OF  SPECIFIC  HAZARDS 

PARAMETER (NUM1  =  6) 

PARAMETER (NUM2  =  13) 

PARAMETER (NUM3  =  4) 

PARAMETER (NUM4  =  4) 

CHARACTER* 1  OKLSTl(NUMl) 

CHARACTER* 1  OKLST2(NUM2) 

CHARACTER* 1  OKLST3(NUM3) 

CHARACTER* 1  OKLST4(NUM4) 

CHARACTER* 1  KEY, KKEY , HAZTYP , HUMTYP 
CHARACTER* 1  CHDOL,CVALUE 
CHARACTER* 1  ICONT(5) 

CHARACTER* 12  PTSFIL(4) 

CHARACTER* 12  TMPFIL 
CHARACTER* 50  COMM(2) 

CHARACTER* 50  PROMPT 

DIMENSION  CONT (  9 )  , CONTV (7 ) ,CONTX(8) 

COMMON  /  CKEY/  KEY, KKEY 

COMMON  /CONSTS/  PI , RHO , FPSPKN , RHOD2 , DRC 
COMMON  /HELGEO/  H , DL , YSEP, WSPD, RADIUS , SHFTAN, DXO 
COMMON  / INPUTC/  ICONT , COMM, PTSFIL 
COMMON  / INPUTD/  CONT , CONTV , CONTX 

COMMON  /PROFIL/  RJ , ZB J , ZH J , ZMJ , UMJ , ZB , ZH , ZM, UM,CU,CY 
COMMON  /  UNITS/  IOU1 , IOU4 , IOU5 , IOU6 , IOU7 , IOU8 , ICRAPH 

★  *★★*★★*★*★★★***★*★****★**'*'★*■*★*★★★★*★★**'*★★*★★★★★★★★ 

DATA  OKLSTl  /  '  C  '  ,  '  c  '  ,  '  H'  ,  '  m'  ,  '  X'  ,  '  x.'  / 

DATA  OKLST2  /’  ' , ' ' a ' , ' B ' , ' b ' , ' C ' , ' c ' , ' D ' , ' d ' , 

1  'E','e','F','f'/ 

DATA  OKLST3 

DATA  OKLST4  / ' h' , ' 1' , ' S ‘ ' S ' / 

CLEAR  SCREEN  AND  HOME  CURSOR 

ICD  =  0 

CALL  HOMCLS(ICD) 

CALL  LOCATE <  4 , 1 ) 
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55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 
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81 
82 

83 

84 

85 

86 

87 
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89 
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92 
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94 

95 
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103 
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C  DETERMINE  THE  TYPE  OF  HAZARD  ANALYSIS 

C  OPTION  THAT  WILL  BE  EXECUTED 

C  - 

C 

WRITE ( IOU1 ,10) 

10  FORMAT  (  25X, '  SELECT  TYPE  OF  HAZARD',///, 

1  1 8 X, 'OVERTURNING  FORCE/MOMENT,  TYPE  <M> ' , / , 

2  18X,  'PARTICULATE  CLOUDS,  TYPE  <C>  ’  , / , 

3  18X, 'TO  EXIT  PROGRAM,  TYPE  <X>’,//) 

C 

C  - 

C  INQUIRE,  OBTAIN,  AND  CHECK  FOR  VALID  MENU  CODE 

C  - 

C 

11  CONTINUE 
C 

WRITE ( IOU1 , ' (23X,A,$) ' )  '  ENTER  HAZARD  CODE  ==>  ' 

C 

READ ( IOUl , ' ( A1 ) ' )  HAZTYP 
C 

IF ( LEGAL ( HAZTYP , IOUl , OKLST1 , NUMl ) .EQ.l)GOTO  11 
C 

C  - 

C  CORRECT  LOWER  CASE  LETTERS  TO  UPPER  CASE 

C  TO  USE  AS  VALID  FLAGS  IN  PARENT  SUBROUTINE 

C  - 

C 

IF ( HAZTYP. EQ. 'c' )  HAZTYP  =  'C' 

IF ( HAZTYP . EQ. 'm' )  HAZTYP  =  'M' 

IF ( HAZTYP. EQ. 'x' )  HAZTYP  =  'X' 

C 

C  - 

C  BRANCH  IF  EXIT  OPTION  CHOSEN 

C  - 

C 

IF ( HAZTYP . EQ . ' M ' ) GOTO  1 8 
IF( HAZTYP. EQ. 'C ' )GOTO  18 
IF ( HAZTYP . EQ . ' X ' ) GOTO  999 
C 

GOTO  11 
C 

18  CONTINUE 
C 

Q  ********************************* 

C  RADIAL  WALL  JET  FLOW  INFORMATION 

Q  ********************************* 

C 

c  - 

C  ACCELERATED  SLIPSTREAM  MEAN  VELOCITY 

C  - 

C 

UN  =  SQRT ( 2 . 0  *DL/RHO ) 

C 

c  - 
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109  C  GROUND  EFFECT  CORRECTION 

110  C  - 

111  C 

112  AKG  =  1.0  -  0 . 9 *EXP ( -2 . 0*H ) 

113  C 

114  C  - 

115  C  MEAN  VELOCITY  AT  ROTOR  DISK  ( RATIOED  TO  UN) 

116  C  - 

117  C 

118  UB  =  AKG/2.0 

119  C 

120  C  - 

121  C  FIND  INITIAL  RADIUS  OF  WALL  JET 

122  C  - 

123  C 

124  CALL  WALJET ( H , UB , UN , UMB ) 

125  C 

126  500  CONTINUE 

127  C 

128  IF(KEY.EQ. 'X' )GOTO  999 

129  IF ( KEY . EQ . ' N ' ) GOTO  999 

130  C 

131  c  - 

132  C  BRANCH  IF  CLOUD  OPTION  CHOSEN 

133  c  - 

134  C 

135  IF ( HAZTYP . EQ . ' C ' ) GOTO  800 

136  C 

137  c  *************************** 

138  C  OVERTURNING  FORCES /MOMENTS 

139  c  *************************** 

140  C 

141  c  - 

142  C  DETERMINE: 

143  C 

144  C  1.  THE  AXIS  ALONG  WHICH  THE  OVERTURNING 

145  C  FORCES /MOMENTS  WILL  BE  CALCULATED 

146  C 

147  C  2.  THE  SIZE  OF  THE  PERSON  AFFECTED 

148  C 

149  C  3.  THE  DISTANCES  AT  WHICH  THE  OVERTURNING 

150  C  FORCES /MOMENTS  WILL  BE  CALCULATED 

151  c  - 

152  C 

153  20  CONTINUE 

154  C 

155  ICD  =  0 

156  CALL  HOMCLS ( ICD ) 

157  CALL  LOCATE (2, 1 ) 

158  C 

159  WRITE ( IOU1 , 12 ) 

160  12  FORMAT (  20X,'  OVERTURNING  FORCE/MOMENT  DATA  MENU',//, 

161  1  1  OX,  'CODE  PARAMETER  VALUE', 

162  2  '  UNITS',/) 
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163  C 

164  C  - 

165  C  PRINT  OUT  MENU  VARIABLES  AS  BASED  ON  THE  WALL  JET 

166  C  OPTION  OR  INTERACTION  PLANE  OPTION  SWITCH  SETTING 

167  C  - 

168  C 

169  IF(ICONT(4) .  EQ. 'W' )THEN 

170  C 

171  WRITE (IOU1, 14)  ICONT (  4  )  , ICONT ( 5 ) , PTSFIL ( 3 ) , 

172  1  CONTV ( 4 ) , CONTV ( 5 ) , CONTV ( 6 ) 

173  C 

174  ELSE 

175  C 

176  WRITE (IOU1, 14)  ICONT ( 4 ), ICONT ( 5 ), PTSFIL ( 4 ) , 

177  1  CONTV ( 4 ) , CONTV {  5  )  , CONTV ( 6 ) 

178  C 

179  ENDIF 

180  C 

181  14  FORMAT (  11X,'A  <W>ALL  JET  OR  <I>I.  TERACTION  PLANE 5X,A2 ,/ , 

182  1  1 IX , ' B  <L>ARGE  OR  <S>MALL  PERSON  ’,5X,A2,/, 

183  2  11X,'C  DATA  OUTPUT  FILENAME  ',2X,A12,/, 

184  3  11X,'  ',/, 

185  4  1 IX ,  '  D  J.NITIAL  STATION  POSITION  ,5X,F7.2,4X,  1  FT '  ,  / , 

186  5  1 IX, ‘ E  HORIZONTAL  INCREMENT  ’ ,5X,F7.2,4X, ’ FT ' , / , 

187  6  1 IX, ' F  MAXIMUM  STATION  POSITION  1 , 5X,F7 . 2 , 4X, ' FT ’ , // ) 

188  C 

189  c  - 

190  C  PROMPT  FOR,  OBTAIN,  AND  CHECK  FOR  LEGAL  INPUT  DATA 

191  C  - 

192  C 

193  16  CONTINUE 

194  C 

195  WRITE ( IOU1 , ’ ( 8X , A , $ ) ' ) 

196  1'  ENTER  DATA  ENTRY  CODE  OR  <RETURN>  TO  CONTINUE  ==>  ' 

197  C 

198  READ ( IOU1 , ' { A 1 ) ’ )  CHDOL 

199  C 

200  IF ( LEGAL ( CHDOL, IOU1 , OKLST2 , NUM2 ) . EQ . 1 ) GOTO  16 

201  C 

202  C  - 

203  C  DIRECT  OPTIONS  BASED  ON  CHOICE  FOR  "CHDOL'' 

204  C  - — 

205  C 

206  IF (CHDOL . EQ .  ‘)GOTO  30 

207  C 

208  C  - 

209  C  CHOOSE  WALJET  OR  IPLANE  OPTION 

210  C  - 

211  C 

212  IF ( CHDOL . EQ . 'A' . OR . CHDOL . EQ . 'a' ) THEN 

213  C 

214  40  CONTINUE 

215  C 

216  WRITE ( IOU1 , ' (/,35X,A,1X,A2/) ' )  '  ANALYSIS  TYPE  =  ’,ICONT(4) 
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217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 
261 
262 
26  J 

264 

265 

266 

267 

268 

269 

270 


C 

c 


c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 

c 

r 


WRITE(IOUl, ' (37X,A,$) ' )  '  ENTER  NEW  CODE  ==> 

READ ( IOU1 , ' ( A1 ) ' )  CVALUE 

IF ( LEGAL (CVALUE , IOU1 , OKLST3 , NUM3 ) .EQ.l)GOTO  40 
ICONT ( 4 )  =  CVALUE 

IF ( ICONT ( 4 )  . EQ . ’ w ' )  ICONT(4)  =  'W' 

IF ( ICONT ( 4 ) . EQ . ' i ' )  ICONT(4)  =  'I' 

GOTO  20 

ENDIF 


CHOOSE  LARGE  OR  SMALL  PERSON 


IF(CHDOL.EQ. 'B’ . OR . CHDOL . EQ . 'b‘ ) THEN 
41  CONTINUE 

WRITE ( IOU 1, ' (/,35X,A,1X,A2/) ’ )  '  PERSON  TYPE  = 

WRITE ( IOU1 , ' ( 3  7X , A , $ ) ' )  '  ENTER  NEW  CODE  ==> 

READ( IOU1 , ' (A1 ) ' )  CVALUE 

IF (LEGAL (CVALUE, IOU 1 ,OKLST4 ,NUM4 ) .EQ.l)GOTO  41 
ICONT ( 5 )  =  CVALUE 

IF(ICONT(5) .EQ. '1 ' )  ICONT (5)  =  'L' 

IF< ICONT ( 5) .EQ. ' 3 ' )  ICONT) 5)  ~  'S’ 

GOTO  20 

ENDIF 


C  CHOOSE  GRAPHICS  FILENAME 

C  - 

C 

IF ( CHDOL. EQ. 'C' . OR . CHDOL . EQ . 'c' ) THEN 
C 

I F ( ICONT ( 4 ) .EQ. 'W* ) THEN 

WRITE (IOU1,  (/,25X,A,1X,A12/) ' ) 

1  '  FILENAME  =  ’,PTSFIL(3) 

ELSE 

WRITE ( IOU 1 , ' ( / ,25X,A, IX, A12/ ) ' ) 

1  '  FILENAME  =  ',PTSFIL(4) 

END  I F 
C 

WRITE ( IOU 1 , ' ( 20X, A, $ ) ' ) 

1  '  ENTER  NEW  FILENAME  ( xxxxxxxx . xxx )  ==> 
C 

READ ( IOU1 , ' ( A1 2 ) ' )  TMPFIL 
C 

IT ( ICONT ( 4 ) . EQ . ' W ' )  PTSFIL(3)  =  TMPFIL 
IF( ICONT( 4 ) .EQ. ’ I ' )  PTSFIL(4)  =  TMPFIL 


’  , ICONT ( 5 ) 
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271 

272 

273 

274 

275 

276 

277 

278 

279 

280 
281 
282 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 

296 

297 

298 

299 

300 

301 

302 
301 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 

322 

323 

324 


C 

C 

C 

C 

C 

C 

C 


c 


c 

c 

c 

c 

c 

c 

c 


c 


c 

c 

c 

c 

c 

c 


c 


c 


c 

c 


c 


c 


c 


GOTO  20 
ENDIF 


CHOOSE  INITIAL  STATION  POSITION 


IF (CHDOL . EQ . 'D' . OR . CHDOL . EQ . 'd* ) THEN 

PROMPT  =  'INITIAL  STATION  POSITION 
CALL  FREAD ( IOU1 , PROMPT , CONTV ( 4 ) , 1 . 0 

IF ( CONTV ( 4 ) . LT .0.0)  CONTV(4)  =0.0 
GOTO  20 

ENDIF 


CHOOSE  HORIZONTAL  INCREMENT 


IF ( CHDOL . EQ . 'E' . OR . CHDOL . EQ . 'e' ) THEN 

PROMPT  =  ' HORIZONTAL  INCREMENT  =  ' 
CALL  FREAD ( TOU1 , PROMPT , CONTV { 5 ) ,1.0 

IF (CONTV ( 5 ) .LT .0.0)  CONTV ( 5 )  =  0.0 
GOTO  20 

ENDIF 


CHOOSE  HORIZONTAL  INCREMENT 


IF (CHDOL . EQ . 'F' . OR . CHDOL . EQ . ' f ' ) THEN 

PROMPT  =  'MAXIMUM  STATION  POSITION 
CALL  FREAD ( IOU1 , PROMPT , CONTV ( 6 } , 1 . 0 

IF ( CONTV ( 6 ) . LT . CONTV (4)1  CONTV(6)  = 
GOTO  20 

ENDIF 

GOTO  20 

30  CONTINUE 

ICD  -  0 

CALL  HOMCLS ( ICD ) 

IF ( ICOHT ( 4 ) . EQ . ' I ' ) GOTO  700 


) 


) 

CONTV ( 4 ) 
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325  C 

326  600  CONTINUE 

327  C 

328  C  - 

329  C  WALL  JET  REGION 

330  C 

331  C  OBTAIN  DATA  FOR  THE  HWJVEL  OPTION 

332  C  - 

333  C 

334  RVZ  =  ( CONTV ( 4 )  -  DXO) /RADIUS 

335  DELH  =  CONTV(5) 

336  HMAX  =  CONTV ( 6 ) 

337  HUMTYP  =  ICONT(5) 

338  C 

339  c  - 

340  C  GENERATE  VELOCITY  PROFILE  AT  RVZ  IN  WALL  JET  REGION 

341  c  - 

342  C 

343  CALL  HWJVEL (H, UN, UMB, RVZ, RADIUS, WSPD, DELH, HMAX, HUMTYP, DXO) 

344  C 

345  GOTO  500 

346  C 

347  700  CONTINUE 

348  C 

349  C  - 

350  C  INTERACTION  PLANE  UPWASH  DEFLECTION  ZONE 

351  C 

352  C  OBTAIN  INPUT  DATA  FOR  THE  HIPVEL  OPTION 

353  c  - 

354  C 

355  XIP  =  ( CONTV ( 4 )  -  DXO)  ./RADIUS 

356  DELH  =  CONTV(5) 

357  HMAX  =  CONTV ( 6 ) 

358  HUMTYP  =  ICONT(5) 

359  C 

360  C  - 

361  C  GENERATE  VELOCITY  PROFILE  AT  XIP  IN  INTERACTION  PLANE 

362  C  - 

363  C 

364  CALL  HIPVEL(H, UN, RADIUS, UMB, XIP, YSEP, WSPD, DELH, HMAX, 

365  *  HUMTYP, DXO) 

366  C 

367  GOTO  500 

368  C 

369  800  CONTINUE 

370  C 

371  c  ■***  +  +  **•»******»****■.**,•«•*  +  *****  +  ******» 

372  C  CALCULATE  PARTICULATE  CLOUD  BOUNDARIES 

373  C  ***  +  *♦*********  +  ***♦**♦***■>♦  +  ■.  *  *  * 

374  C 

375  CALL  CLOUD (UN, UMB) 

376  C 

377  GOTO  500 

378  C 
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379 

380 

381 

382 

383 

384 

385 

386 

387 

388 


C  - 

C  NORMAL  PROGRAM  EXIT 

C  - 

C 

999  CONTINUE 
C 

RETURN 

END 

C 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 


C 

C 

SUBROUTINE  HIPVEL ( H , UN , RADIUS , UMB , XIP , YSEP , WSPD , DELH , 

*  HMAX , HUMTYP , DXO ) 

C 

c 

C  SUBROUTINE  HIPVEL  GENERATES  THE  VELOCITY  PROFILE  AND  THE  FORCES 

C  AND  OVERTURNING  MOMENTS  FOR  A  HUMAN  BEING  ALONG  THE  INTERACTION 

C  PLANE  FOR  THE  TWIN  ROTOR  CASE 

C 

c 

CHARACTER* 1  TEMCHAR 
CHARACTER* 1  KEY , KKEY , HUMTYP 
CHARACTER* 1  ICONT(5) 

CHARACTER* 12  PTSFIL(4) 

CHARACTER* 50  COMM(2) 

C 

COMMON  /  CKEY/  KEY, KKEY 
COMMON  /CONSTS/  PI , RHO , FPSPKN , RHOD2 , DRC 
COMMON  /INPUT C/  ICONT , COMM, PTSFIL 
COMMON  /PERSON/  QP(12),DSET 

COMMON  /PROFIL/  RJ , ZBJ , ZHJ , ZMJ , UMJ , ZB, ZH , ZM, UM, CU , CY 
COMMON  /  UNITS/  IOU1 , IOU4 , IOU5 , IOU6 , IOU7 , IOU8 , IGRAPH 
C 

c 

c 

ICD  =  0 

CALL  HOMCLS ( ICD ) 

C 

c  - 

C  INPUT  FOR  DELH  AND  HMAX  COMES  FROM  INPUTV  STATUS  MENU 

C  - 

c 

DSET  =  DELH 

IF ( DSET . EQ . 0 . ) DELH  =  HMAX 
C 

DELH  =  DELH/RADIUS 
HMAX  =  HMAX  /  RAD  I  US 

NHPTS  =  I F I X ( ( HMAX  -  XIP) /DELH)  +  1 
C 

IF ( DSET . EQ . 0 . ) GOTO  3  3 
C 

C  - 

C  WRITE  OUTPUT  HEADER  (FOR  PLOT  FILE,  SEE  BELOW) 

C  - 

c 

IF ( IOU6 .NE.IOU1 )  WRITE ( IOU6 , ' ( ’ ’ 1 * ’ ) ’ ) 

C 

IF ( IOU6 . EQ . 6 )  WRITE ( IOU6 ,93)  COMM ( 1 ) ,COMM( 2 ) 

93  FORMAT (  1  OX , A50 , / , 1  OX , A50 , / / ) 

C 
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55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 


WRITE (IOU6, 1001 ) 

1001  FORMAT (  12X, '  SUMMARY  OF  OVERTURNING  FORCES  AND  MOMENTS',//, 

1  19X, ’RADIUS' ,6X, ' TOTF ' ,6X, ' TOTM ' ,/, 

2  2 OX,  ' (FT)  ' ,7X,  ' (LB)  ' ,5X,  ' (FT-LB)  ’ ,/) 

33  CONTINUE 

C 

c  - 

c  WRITE  OUT  GRAPHICS  FILES  IF  SWITCH  IS  SET  BY  USER 

C  - - - 

c 


c 

c 

c 

c 

c 

c 


c 


c 


c 


c 


IF ( IGRAPH .  EQ .  1  ) THEN 


OPEN  GRAPHICS  FILE  AND  WRITE  FILE  HEADER 


OPEN ( IOU8 , FILE=PTSFIL ( 4 ) , STATUS= ' NEW ' ,ERR=2000) 

WRITE (IOU8, 83)  COMM ( 1 ) , COMM ( 2 ) 

83  FORMAT (  10X,A50,/,10X,A50,//) 

WRITE (IOU8, 80) 

80  FORMAT (  IX , ' TITLE= "TWIN  ROTOR  DAIP  DATA"’) 

WRITE ( IOU8 ,81) 

81  FORMAT (  IX, 'VARIABLES  =  DAIP , TOTF , TOTM ' ) 

WRITE (IOU8, 88) 

88  FORMAT (  IX, 'ZONE  T  =  "GW  =  xxxxx  LB,  WAGL  =  xx  FT",', 
*  ’  I— X ,  F=POINT ' ) 


ENDIF 


C 

C  - 

C  BEGIN  LOOP  INCREMENTING  THE  RADIAL  POINTS  AT  WHICH 

C  THE  OVERTURNING  MOMENT  CALCULATIONS  ARE  MADE 

C  - 

C 

DO  565  K  =  1 , NHPTS 
C 

c  - 

C  TF  IS  INTERACTION  PLANE  AMPLIFICATION  FACTOR 

C  - 

C 

TF  =  1.55  -  (0.55)*EXP(-1 . 35*XIP ) 

C 

c  - 

C  GET  PARAMETERS  AT  BASE  RADIUS  FOR  'BOUNDARY  LAYER' 

C  - 

c 


RIP0  =  SQRT ( XIP*  *  2  +  YSEP  *  *  2 ) 

C 

c  - 

C  ' PROPRM '  PROVIDES  THE  VELOCITY  PROFILE  PARAMETERS 
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109  C  OF  A  RADIAL  WALL  JET  (WITHOUT  INTERACTION  PLANE) 

110  C  - 

111  C 

112  CALL  PROPRM ( H , UMB , R I P  0 ) 

113  C 

114  ZIPB  =  ZB 

115  ZIPM  =  ZM 

116  ZIPH  =  ZH 

117  C 

118  RIPM  =  SQRT(XIP**2  +  (YSEP  +  ZIPM)**2) 

119  C 

120  CALL  PROPRM (H, UMB, RIPM) 

121  C 

122  UMM  =  UM 

123  C 

124  C  - 

125  C  OUTPUT  HEADER 

126  C  - 

127  C 

128  IF ( DSET . NF . 0 . ) GOTO  78 

129  C 

130  XXIP  =  RADIUS*XIP 

131  XIPOUT  =  XXIP  +  DXO 

132  C 

133  IF ( IOU6 .NE . IOU1 )  WRITE ( IOU6 1 '')' ) 

134  C 

135  IF ( IOU6 . EQ . 6 )  WRITE ( IOU6 , 93)  COMM ( 1 ) , COMM( 2 ) 

136  C 

137  WRITE(IOU6, 1000)  XIPOUT 

138  1000  FORMAT (  2X,’TWIN  ROTOR  INTERACTION  PLANE  VELOCITY  PROFILE’, 

139  1  'AT  DISTANCE  =  ’,F7.1,’  FT’,//) 

140  C 

141  WRITE ( IOU6 ,1002) 

142  1002  FORMAT (  3X ,’ HEIGHT ’, 6X , ’MEAN  VELOCITY ’, 7X, ’ PEAK  VELOCITY ’, 6X, 

143  1  ’MEAN  Q ’ , 4X, ’ PEAK  Q  ’  ,  /  , 

144  2  3X, ’ (FT) ’ ,7X, ’ (FPS) ’ ,6X, ’ (KN) ’ ,5X, ’ (FPS) ’ ,6X, ’ (KN) ’ ,5X, 

145  3  ’ (PSF) ’ ,5X, ’ (PSF) ’ ,/) 

146  78  CONTINUE 

147  C 

148  C  - 

149  C  ’AN’  IS  ACTUALLY  ’=  1. 0/7.0’ 

150  C  - 

151  C 

152  AN  =  0.142857142 

153  C 

154  DELZ  =  0.5/RADIUS 

155  NPTS  =  12 

156  C 

157  DO  500  I  =  1 , NPTS 

158  C 

159  ZIP  =  DELZ  * ( I  -  1)  +  (0.25/RADIUS) 

160  C 

161  C  - 

162  C  GET  MAX  WALL  JET  VELOCITY  AT  EFFECTIVE  RADIUS 
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163  C 

164  C 

165 

166  C 

167 

168  C 

169 

170 

171  C 

172  C 

173  C 

174  C 

175  C 

176 

177 

178 

179  C 

180  C 

181  C 

182  C 

183  C 

184 

185 

186  C 

187 

188  C 

189  C 

190  C 

191  C 

192  C 

193 

194 

195  C 

196  C 

197  C 

198  C 

199  C 

200 

201  C 

202 

203 

204 

205  C 

206  C 

207  C 

208  C 

209  C 

210  C 

211  C 

212 

213  C 

214 

215  C 

216 


RIP  =  SQRT ( XIP*  *2  +  ( YSEP  +  ZIP)**2) 

CALL  PROPRM( H , UMB , RIP ) 

VN  =  UN 
VZ  =  UM 

INTERACTION  PLANE  'BOUNDARY  LAYER' 

IF ( ZIP • LT . Z I PM) THEN 

VZ  =  UMM* ( Z IP/Z IPM) *  *AN 
ENDIF 

DEVELOPED  INTERACTION  PLANE  JET 

VH  =  TF*VZ  *XIP/RIP 
W  =  TF*VZ  * (YSEP  +  ZIP) /RIP 

ZZ  =  ZIP*RADIUS 

MEAN  HORIZONTAL  VELOCITIES  AND  DYNAMIC  PRESSURE 


VHMF  =  VH*UN 
VHMK  =  VHMF/FPSPKN 


PEAK  HORIZONTAL  VELOCITIES  AND  DYNAMIC  PRESSURE 


VMFD3I  =  (XIP*0 . 2444 )  +  0.8 

IF(VMFD3I .GT. 2.5)  VMFD3I  =  2.5 
VHPF  =  VMFD3I *VHMF 
VHPK  =  VHPF/FPSPKN 


THE  EFFECT  OF  WIND  IS  TO  ADD  (DOWNWIND  SIDE)  OR  SUBTRACT 
(UPWIND  SIDE)  ' XWK '  TIMES  THE  AMBIENT  WIND  VELOCITY  TO 
THE  HORIZONTAL  PROFILE  VELOCITY  (EMPIRICAL,  CH-5",E  BASED) 


XKW  =  (-0 .5*H)  +  2.5 
IF ( XKW . LT . 1 . 0) XKW  =  1.0 
WSPD2  =  WSPD*XKW 
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217  VHMK  =  VHMK  +  WSPD2 

218  VHMF  =  VHMK*FPSPKN 

219  VHPK  =  VHPK  +  WSPD2 

220  VHPF  =  VHPK*FPSPKN 

221  C 

222  C  - 

223  C  DYNAMIC  PRESSURE 

224  C  - 

225  C 

226  QHM  -  RHOD2*YHMF**2 

227  QP(I)  =  RHOD2*VHPF**2 

228  C 

229  IF ( DSET . NE . 0 . ) GOTO  77 

230  C 

231  C  - 

232  C  REPORT  HORIZONTAL  COMPONENTS 

233  C  - 

234  C 

235  WRITE ( IOU6 ,1003)  ZZ , VHMF , VHMK, VHPF , VHPK, QHM, QP< I ) 

236  1003  FORMAT  (  F8 . 2 , 2X , 6F 1 0 . 3 ) 

237  77  CONTINUE 

238  C 

239  500  CONTINUE 

240  C 

241  IF ( DSET . NE . 0 . ) GOTO  520 

242  WRITE (IOU1 ,73) 

243  73  FORMAT (  ) 

244  C 

245  WRITE ( IOU1 » ‘ ( 1 9X , A , $ ) ' ) 

246  1  TYPE  <RETURN>  TO  CONTINUE 

247  C 

248  READ ( IOU1 , ’ ( A1 ) ’ )  TEMCHAR 

249  C 

250  ICD  =  0 

251  CALL  HOMCLS ( ICD ) 

252  IF  ( IOU6  .NE  .  IOU1 )  WRITE  { IOU6  1 '')' ) 

253  C 

254  IF ( IOU6 .  EQ . 6 )  WRITE ( IOU6 , 9 3 )  COMM ( 1 ) ,COMM( 2 ) 

255  C 

256  WRITE(IOU6, 1007)  XIPOUT 

257  1007  FORMAT (  12X,'TWIN  ROTOR  FORCE  PROFILE  AT  DISTANCE  = 

258  1  F7.1, '  FT' ,//) 

259  WRITE (IOU6, 1008) 

260  1008  FORMAT (  2X HEIGHT 6X , ' PEAK  O' , 6X , ' FOVER ' , 7X , ' OVERM ' , 7X, 

261  1  'TOT  F ’ , 7X , 'TOT  M ' , / , 

262  2  3X, ' (FT) ‘ ,8X, ’ (PSF) ' ,7X, ’ (LB) ' ,6X, ' (FT-LB) ’ ,7X, 

263  3  ' (LB) ' ,6X, ’ (FT-LB) ’ ,/) 

264  520  CONTINUE 

265  C 

266  C  - 

267  C  CALL  SUBROUTINE  TO  CALCULATE  THE 

268  C  FORCES  AND  MOMENTS  ON  A  HUMAN  BEING 

269  C  - 

270  C 


HIPVEL-5 


SUBROUTINE  HIPVEL 


271 

CALL  MOMENT { NPTS , HUMTYP , TOTF , TOTM ) 

272 

C 

273 

IF ( DSET . EQ . 0 . ) GOTO  545 

274 

C 

275 

HH  =  XIP*RADIUS 

276 

HHOUT  =  HH  +  DXO 

277 

C 

278 

WRITE ( IOU6 ,1014)  HHOUT, TOTF, TOTM 

279 

1014 

FORMAT (  18X,F8 . 2 , 2F10 . 3 ) 

280 

C 

281 

IF ( IGRAPH . EQ . 1 ) THEN 

282 

c 

283 

WRITE (IOU8, 90)  HHOUT , TOTF , TOTM 

284 

90 

FORMAT (  1X,F7.2,1X,F7.2,1X,F8.2) 

285 

c 

286 

ENDIF 

287 

c 

288 

545 

CONTINUE 

289 

c 

290 

XIP  =  XIP  +  DELH 

291 

c 

292 

565 

CONTINUE 

293  C 

294  C 

295  C 

296  C 

297  C 

298 

299  C 

300 

301  C 

302 

303  C 

304 

305  C 

306 

307  C 

308  C 

309  C 

310  C 

311  C 

312  C 

313  2000  CONTINUE 

314  C 

315  CALL  HOMCLS ( 0 ) 

316  WRITE( IOU1 ,2001 ) 

317  2001  FORMAT (  ////,8X, 

318  1  '  ***  ERROR  ***  PLEASE  CHOOSE  A  NEW  OUTPUT  FILENAME', 

319  2  /// , 8X ,  '  TYPE  <RETURN>  TO  CONTINUE  ',$) 

320  READ ( IOU1 , ’ ( A1 ) ' )  TEMCHAR 

321  KEY  =  'P' 

322  C 

323  999  CONTINUE 

324  C 


CLOSE  AN  OPEN  GRAPHICS  FILE 

IF ( IGRAPH . EQ . 1 ) THEN 

CLOSE ( IOU8 , STATUS= ' KEEP ' ) 

ENDIF 

CALL  INKEY 
GOTO  999 

THE  ERROR  LOGIC  ALLOWS  FOR  THE  HANDLING  OF  FILE 
OPEN  ERRORS  BY  RETURNING  THE  USER  TO  A  MENU 
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325  RETURN 

326  END 

327  C 

328 
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SUBROUTINE  HOMCLS 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


SUBROUTINE  HOMCLS (CODE) 

********************************************************* 
SUBROUTINE  HOMCLS 

THIS  SUBROUTINE  HOMES  THE  CURSOR  AND  CLEARS  THE  TERMINAL- 
SCREEN  ( CODE=0 )  OR  HOMES  THE  CURSOR  ONLY  (CODE=l) 


COMMON  /  UNITS/  IOU1 , IOU4 , IOU5 , IOU6 , IOU7 , IOU8 , IGRAPH 

INTEGER* 4  CODE 

CHARACTER* 4  ED 
CHARACTER* 1  EED ( 4 ) 

EQUIVALENCE  ( ED , EED ( 1 )  ) 

CHARACTER* 3  EE 
CHARACTER* 1  EEE ( 3 ) 

EQUIVALENCE  ( EE , EEE ( 1 )  ) 

IF ( CODE. EQ. 1 ) GOTO  20 


HOME  CURSOR  AND  CLEAR  SCREEN 
ANSI  CONTROL  SEQUENCE:  ED  =  ESC[2J 


EED ( 1 )  =  CHAR (27) 
EED ( 2 )  =  CHAR (91) 
EED ( 3 )  =  CHAR( 50 ) 
EED ( 4 )  =  CHAR (74) 

WRITE ( IOU1 ,* )  ED 

20  CONTINUE 


HOME  CURSOR  ONLY 

ANSI  CONTROL  SEQUENCE:  EE  =  ESC[H 


EEE ( 1 )  =  CHAR (27) 
EEE ( 2  )  =  CHAR (91) 
EEE ( 3 )  =  CHAR( 72 ) 

WRITE ( IOU 1 , * )  EE 

RETURN 

END 
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55  C 

56 


HOMCLS-2 
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1  C 

2  C 

3  SUBROUTINE  HSVTX ( XT ,YT, ZT,VX,VY ,VZ , GAMMA, RADIUS ) 

4  C 

5  C  ****************************************************************** 

6  C  SUBROUTINE  HSVTX 

7  C 

8  C  THIS  SUBROUTINE  DIRECTS  THE  CALCULATION  OF  THE  INDUCED  VELOCITY 

9  C  FIELD  DUE  TO  A  HORSESHOE  VORTEX  SYSTEM  OF  UNIT  STRENGTH.  POINT  1 

10  C  (LEFT  =  LI,  RIGHT  =  R1 )  DEFINE  THE  EXTENT  OF  THE  BOUND  PORTION  OF 

11  C  THE  HORSESHOE.  THE  TRAILERS  START  AT  POINT  1  AND  EXTEND  THROUGH 

12  C  POINT  2,  AND  THEN  ON  TO  POINT  3.  THIS  ALLOWS  TWO  ELEMENTS  FOR 

13  C  EACH  TRAILER  SO  THAT  IT  CAN  'BEND'  TO  ACCOUNT  FOR  GROUND  CONTACT. 

^4  q  ****************************************************************** 

15  C 

16  COMMON  / CHSVTX/  XL1 , YLl , ZLl , XL2 , YL2 , ZL2 , XL3 , YL3 , ZL3 , 

17  1  XR1 , YR1 , 2R1 , XR2 , YR2 , ZR2 , XR3 , YR3 , ZR3 

18  COMMON  /CVLINE/  IFI , XA, YA, 2A, XB , YB , ZB , XC , YC, ZC , Ql , Q2 ,Q3 

19  C 

20  C  ****************************************************************** 

21  C 

22  C  - 

23  C  AT  SPECIFIED  (X,Y,Z)  TARGET  POINT  IN  VICINITY 

24  C  OF  ROTOR,  CALCULATE  THE  VECTOR  VELOCITY 


25 

c 

— 

26 

c 

27 

VX  « 

O 

o 

28 

VY  = 

o 

o 

29 

vz  « 

0.0 

30 

c 

31 

XC  = 

XT 

32 

►2 

O 

H 

YT 

33 

ZC  = 

ZT 

34 

c 

35 

C 

36 

c 

LEFT 

TRAILER  CONTRIBUTION,  POINT  1  TO  POINT  2 

37 

38 

c 

39 

IFI  = 

=  0 

40 

XA  = 

XL  1 

41 

YA  = 

YLl 

42 

ZA  * 

ZLl 

43 

XB  = 

XL2 

44 

YB  * 

YL2 

45 

ZB  = 

ZL2 

46 

c 

47 

CALL 

VLINE 

48 

c 

49 

VX  = 

VX  -  Ql 

50 

VY  = 

VY  -  Q2 

51 

VZ  * 

VZ  -  Q3 

52 

c 

r 

54  C  LEFT  TRAILER  IMAGE 


HSVTX- 1 


SUBROUTINE  HSVTX 


55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 


C 

c 


c 

c 


c 

c 

c 

c 

c 


c 

c 


c 

c 

c 

c 

c 


c 

c 


c 

c 

c 

c 

c 


ZA  =  -ZA 
ZB  =  -ZB 

CALL  VLINE 

VX  =  VX  +  01 
VY  =  VY  +  02 
VZ  =  VZ  +  Q3 


LEFT  TRAILER  CONTRIBUTION,  POINT  2  TO  POINT  3 


IFI 

=  1 

XA 

= 

XL2 

YA 

= 

YL2 

ZA 

= 

ZL2 

XB 

- 

XL  3 

YB 

= 

YL3 

Z3 

= 

ZL3 

CALL 

VLINE 

VX 

= 

VX  - 

Ql 

VY 

= 

VY  - 

<32 

VZ 

= 

VZ  - 

Q3 

LEFT  TRAILER  TMAGE 


ZA  =  -ZA 
ZB  =  -ZB 

CALL  VLINE 

VX  =  VX  +  Q1 
VY  =  VY  +  Q2 
VZ  =  VZ  +  Q3 


PANWISE  VORTEX  CONTRIBUTION 


IFI  =  0 
XA  =  XL  1 
YA  =  YL1 
ZA  =  ZL1 
XB  =  XR1 
YB  =  YR1 
ZB  =  ZR1 
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109 

CALL 

VLINE 

110 

C 

111 

VX  = 

VX  +  Ql 

112 

VY  = 

VY  +  Q2 

113 

VZ  = 

VZ  +  Q3 

114 

C 

115 

c 

— 

— 

116 

c 

SPANWISE  VORTEX  IMAGE 

117 

c 

— 

— 

118 

c 

119 

ZA  = 

-ZA 

120 

ZB  = 

-ZB 

121 

c 

122 

CALL 

VLINE 

123 

c 

124 

VX  = 

VX  -  Ql 

125 

VY  = 

VY  -  Q2 

126 

VZ  = 

VZ  -  Q3 

127 

c 

128 

c 

129 

c 

RIGHT  TRAILER  CONTRIBUTION,  POINT  1  TO  POINT  2 

130 

Q 

131 

c 

132 

IFI 

=  0 

133 

XA  = 

XR1 

134 

YA  = 

YRl 

135 

ZA  = 

ZR1 

136 

XB  = 

XR2 

137 

YB  = 

YR2 

138 

ZB  = 

ZR2 

139 

c 

140 

CALL 

VLINE 

141 

c 

142 

VX  = 

VX  +  Ql 

143 

VY  = 

VY  +  Q2 

144 

VZ  = 

VZ  +  Q3 

145 

c 

146 

c 

147 

c 

RIGHT  TRAILER  IMAGE 

148 

c 

— 

149 

c 

150 

ZA  = 

-ZA 

151 

ZB  = 

-ZB 

152 

c 

153 

CALL 

VLINE 

154 

c 

155 

VX  = 

VX  -  Ql 

156 

VY  = 

VY  -  Q2 

157 

VZ  = 

VZ  -  Q3 

158 

c 

1 

r 

160  C  RIGHT  TRAILER  CONTRIBUTION,  POINT  2  TO  POINT  3 

161  C  - 

162  C 
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163 

IFI  =  1 

164 

XA  =  XR2 

165 

YA  =  YR2 

166 

ZA  =  ZR2 

167 

XB  =  XR3 

168 

YB  =  YR3 

169 

ZB  =  ZR3 

170 

C 

171 

CALL  VLINE 

172 

C 

173 

VX  =  VX  +  Q1 

174 

VY  =  VY  +  Q2 

175 

VZ  =  VZ  +  Q3 

176 

C 

177 

c 

178 

c 

RIGHT  TRAILER  IMAGE 

179 

c 

180 

c 

181 

ZA  =  -2A 

182 

ZB  =  -ZB 

183 

c 

184 

CALL  VLINE 

185 

c 

186 

VX  =  VX  -  Q1 

187 

VY  =  VY  -  Q2 

188 

VZ  =  V2  -  Q3 

189 

c 

190 

c 

191 

c 

DIMENS IONALIZE 

192 

c 

193 

c 

194 

GDR=  GAMMA/RADIUS 

195 

VX  =  VX*GDR 

196 

VY  =  VY*GDR 

197 

VZ  =  VZ*GDR 

198 

c 

199 

RETURN 

200 

END 

201 

c 
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SUBROUTINE  HWJVEL 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 


C 

C 

SUBROUTINE  HWJVEL ( H , UN , UMB , RVZ , RADIUS , WSPD , DELH , HMAX , 

*  HUMTYP , DXO ) 

C 

C  SUBROUTINE  HWJVEL  GENERATES  THE  VELOCITY  PROFILE  AND 

C  THE  FORCES  AND  OVERTURNING  MOMENTS  FOR  A  HUMAN  BEING 

C  AT  A  GIVEN  RADIUS 

C 

CHARACTER* 1  ICONT(5) 

CHARACTER* 1  TEMCHAR 
CHARACTER* 1  KEY ,KKEY, HUMTYP 
CHARACTER* 1 2  PTSFIL<4) 

CHARACTER* 50  COMM(2) 

C 

COMMON  /  CKEY/  KEY , KKEY 
COMMON  /CONSTS/  PI , RHO , FPSPKN , RHOD2 , DRC 
COMMON  / INPUTC/  ICONT , COMM , PTSFIL 
COMMON  /PERSON/  QP(12),DSET 

COMMON  /PROFIL/  RJ ,  ZBJ , ZHJ, ZMJ, UMJ , ZB, ZH , ZM,UM,CU,CY 
COMMON  /  UNITS/  IOU1 , IOU4 , IOU5 , IOU6 , IOU7 , IOU8 , IGRAPH 
C 

c  ***************************************************** 

C 

ICD  =  0 

CALL  HOMCLS ( ICD) 

C 

c  - 

C  INPUT  FOR  DELH  AND  HMAX  COMES  FROM  INPUTV  STATUS  MENU 

C  - 

C 

DSET  =  DELH 
C 

IF ( DSET . EQ . 0 . )  DELH  =  HMAX 
C 

DELH  =  DELH/RADIUS 
HMAX  =  H MAX/RADIUS 

NHPTS  =  IFIX ( ( HMAX  -  RVZ) /DELH)  +  1 
C 

IF ( DSET . EQ . 0 . ) GOTO  50 

IF ( IOU6 .NE.IOU1 )  WRITE (IOU6, ' ( ’ 1 1 ' ' ) ' ) 

C 

IF { IOU6 . EQ . 6 )  WRITE (IOU6, 93)  COMM ( 1 ) , COMM ( 2 ) 

93  FORMAT (  1  OX , A50 , / , 1  OX , A50 , / / ) 

C 

WRITE ( IOU6 ,1001) 

1001  FORMAT (  12X,  SUMMARY  OF  OVERTURNING  FORCES  AND  MOMENTS’,//, 

1  19X, 'RADIUS' ,6X, ’ TOTF ' ,6X, ' TOTM ' ,/, 

2  2 OX, ’ (FT) ' ,7X, ' (LB) ' ,5X, ' (FT-LB) ' ,/) 

50  CONTINUE 

C 

c  - 
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55  C  WRITE  OUT  GRAPHICS  FILE  IF  SWITCH  IS  SET  BY  USER 

56  C  - 

57  C 

58  IF ( IGRAPH . EQ . 1 ) THEN 

59  C 

60  C  - 

61  C  OPEN  GRAPHICS  FILE  AND  WRITE  FILE  HEADER 

62  C  - 

63  C 

64  OPEN ( IOU8 , FILE=PTSFIL ( 3 ) , STATUS= ’ NEW ’ ,ERR=2000) 

65  C 

66  WRITE (IOU8, 83)  COMM( 1 ) , COMM( 2 ) 

67  83  FORMAT (  10X,A50 , / , 10X, A50 , // ) 

68  C 

69  WRITE (IOU8, 80) 

70  80  FORMAT (  IX,  'TITLE=''SINGLE  ROTOR  DFRC  DATA"') 

71  C 

72  WRITE (IOU8, 81 ) 

73  81  FORMAT (  IX, 'VARIABLES  =  DFRC , TOTF , TOTM ' ) 

74  C 

75  WRITE (IOU8, 88) 

76  88  FORMAT (  IX, 'ZONE  T  =  "GW  =  xxxxx  LB,  WAGL  =  XX  FT",', 

77  *  '  I=x,  F=POINT ' ) 

78  C 

79  ENDIF 

80  C 

81  C  - 

82  C  BEGIN  LOOP  INCREMENTING  THE  RADIAL  POINTS  AT  WHICH 

83  C  THE  OVERTURNING  MOMENT  CALCULATIONS  ARE  MADE 

84  C  - 

85  C 

86  DO  565  K  =  1,NHPTS 

87  C 

88  C  - 

89  C  ' PROPRM '  PROVIDES  THE  VELOCITY  PROFILE  PARAMETERS 

90  C  OF  A  RADIAL  WALL  JET  (WITHOUT  INTERACTION  PLANE) 

91  c  - 

92  C 

93  CALL  PROPRM ( H, UMB, RVZ ) 

94  C 

95  ZETAM  =  ZM/ZB 

96  C 

97  c  - 

98  C  BOUNDARY  LAYER  REGION  EXPONENT 

99  c  - 

100  C 

101  C  - 

102  C  'AN'  IS  ACTUALLY  '=  1.  0/7.0’ 

103  C  - 

104  C 

105  AN  =  0.142857142 

106  C 

107  C  - 

108  C  SHEAR  LAYER  REGION  EXPONENT,  TO  MEET  EDGE  CONDITIONS 
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109  C  (FROM  FIGURE  7,  USAAVLABS  TECHNICAL  REPORT  68-52,  JULY  1968) 

110  C  - 

111  C 

112  ALPW  «  ALOG( 1 . 0  -  1.0/ SQRT (2.0))/ ALOG ( ( 2 H  -  ZM)/(ZB  -  ZM) ) 

113  C 

114  VN  =  UN 

115  VMN  =  UM 

116  C 

117  C  - 

118  C  CALCULATION  OF  THE  PEAK  VELOCITY  CORRECTION  FACTOR 

119  C  IS  MADE  IN  THE  NEXT  SECTION  AT  THE  3  FT  POSITION 

120  C  - 

121  C 

122  VZ3  =  3.0/RADIUS 

123  ZETA3  =  VZ3/ZB 

124  VZM3  =  0.0 

125  IF ( ZETA3 .GE . ZETAM) GOTO  51 

126  C 

127  C  - 

128  C  THE  3  FT  HEIGHT  IS  IN  THE  BOUNDARY  LAYER 

129  C  - 

130  C 

131  IF ( ZETAM . GT .0.0)  VZM3  =  ( ZETA3/ZETAM) **AN 

132  GOTO  52 

133  C 

134  51  CONTINUE 

135  C 

136  C  - 

137  C  THE  3  FT  HEIGHT  IS  IN  THE  SHEAR  LAYER 

138  C  - 

139  C 

140  IF ( VZ3 .LE . ZB) THEN 

141  VZM3  =  (1.0  -  ( ( ZETA3  -  ZETAM)/ (1.0  -  ZETAM) ) **ALPW) **2 

142  ENDIF 

143  C 

144  52  CONTINUE 

145  C 

146  C  - 

147  C  THE  DELTA  VELOCITY  TO  ADD  TO  THE  MEAN  VELOCITY  IS 

148  C  CALCULATED  AS  THE  CONSTANT  VMFD3  (FT/SEC) 

149  c  - 

150  C 

151  VZN3  =  VZM3  *VMN 

152  VMF3  =  VZN3*VN 

153  VMFD3S  =  <RVZ*0.2444)  +  0.4 

154  C 

155  IF(VMFD3S .GT. 1 . 5)  VMFD3S  =  1.5 

156  C 

157  VMFD3  =  VMFD3S*VMF3 

158  C 

159  IF ( DSET . NE . 0 . ) GOTO  78 

160  C 

161  RRVZ  =  RVZ*RADIUS 

162  RVZOUT  =  RRVZ  +  DXO 
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163  C 

164  IF ( IOU6 .NE . IOU1 )  WRITE ( IOU6 ,'('' 1 '•)' ) 

165  C 

166  IF ( IOU6 . EQ . 6 )  WRITE ( IOU6 , 9 3 )  COMM( 1 ) ,COMM( 2 ) 

167  C 

168  WRITE ( IOU6 ,  1000 )  RVZOUT 

169  1000  FORMAT (  10X,'SINGLE  ROTOR  VELOCITY  PROFILE  AT  RADIUS  = 

170  1  F7 . 1 , '  FT',//) 

171  WRITE (IOU6, 1005) 

172  1005  FORMAT (  2X, 'HEIGHT 5X, ‘MEAN  VELOCITY ', 7X, ' PEAK  VELOCITY 6X, 

173  1  'MEAN  Q’ ,4X, 'PEAK  Q' ,/, 

174  2  3X, ' (FT) ' ,5X, ' (FPS) ‘ ,6X, ' (KN) ' ,5X, ' (FPS) ' ,6X, ' (KN) ' ,5X, 

175  3  ' (PSF) ' ,5X, ' (PSF) ' ,/) 

176  78  CONTINUE 

177  C 

178  C 

179  C 

180  C 

181  C 

182  C 

183 

184 

185  C 

186 

187  C 

188 

189 

190  C 

191 

192  C 

193  C 

194  C 

195  C 

196  C 

197 

198 

199  C 

200 

201  C 

202  300  CONTINUE 

203  C 

204  C  - 

205  C  Z  IS  WITHIN  SHEAR  LAYER 

206  C  - 

207  C 

208  VZM  =  0.0 

209  C 

210  IF ( Z .LE . ZB ) THEN 

211  VZM  =  (1.0  -  ((ZETA  -  ZETAM)/(1.0  -  ZETAM) ) **ALPW) **2 

212  ENDIF 

213  C 

214  400  CONTINUE 

215  C 

216  VZN  =  VZM*VMN 


SET  UP  ABILITY  TO  CALCULATE  AT  0.5  FT. 
INCREMENTS  UP  THE  VELOCITY  PROFILE 

DELZ  =  0.5 /RADIUS 
NPTS  =12 

DO  500  I  =  1 , NPTS 

Z  =  DELZ* (I  -  1)  +  (0.25/RADIUS) 

ZETA  =  Z/ZB 

IF ( ZETA. GE. ZETAM) GOTO  300 

Z  IS  WITHIN  BOUNDARY  LAYER 

VZM  =0.0 

IF ( ZETAM.GT .0.0)  VZM  =  ( ZETA/ZETAM) **AN 
GOTO  400 
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217  C 

218  C  - 

219  C  DIMENSIONAL  HEIGHT 

220  C  - 

221  C 

222  ZZ  =  Z*RADIUS 

223  C 

224  C  - 

225  C  MEAN  VELOCITIES 

226  C  - 

227  C 

228  VMF  *  VZN*VN 

229  VMK  =  VMF/FPSPKN 

230  C 

231  C  - 

232  C  PEAK  VELOCITIES 

233  C  - 

234  C 

235  VPF  *  VMF  +  VMFD3 

236  VPK  =■  VPF/FPSPKN 

237  C 

238  C  - 

239  C  THE  EFFECT  OF  WIND  IS  TO  ADD  (DOWNWIND  SIDE)  OR  SUBTRACT 

240  C  (UPWIND  SIDE)  ' XWK '  TIMES  THE  AMBIENT  WIND  VELOCITY  TO 

241  C  THE  HORIZONTAL  PROFILE  VELOCITY  (EMPIRICAL,  CH-53E  BASED) 

242  C  - 

243  C 

244  XKW  =  (-0 .5*H)  +  2.5 

245  C 

246  IF ( XKW . LT .1.0)  XKW  =  1.0 

247  C 

248  WSPD2  =  WSPD*XKW 

249  VMK  =  VMK  +  WSPD2 

250  VMF  =  VMK*FPSPKN 

251  VPK  =  VPK  +  WSPD2 

252  VPF  =  VPK*FPSPKN 

253  C 

254  C  - 

255  C  DYNAMIC  PRESSURE 

256  C  - 

257  C 

258  QM  =  RHOD2*VMF**2 

259  QP ( I )  =  RHOD2  *VPF*  *2 

260  C 

261  IF ( DSET . NE . 0 . ) GOTO  77 

262  C 

263  WRITE ( IOU6 ,1002)  ZZ , VMF , VMK, VPF , VPK , QM, QP ( I ) 

264  1002  FORMAT (  F8.2,6F10.3) 

265  77  CONTINUE 

266  C 

267  500  CONTINUE 

268  C 

269  IF ( DSET .NE . 0 . )GOTO  520 

270  C 
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299 
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301 
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303 
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324 


WRITE ( IOU1 / 7  3 ) 

7  3  FORMAT (  ) 

C 

WRITE ( IOU1 ,  ’  ( 1 5X , A, $ )  '  ) 

1  *  TYPE  <RETURN>  TO  CONTINUE  ’ 

C 

READ( IOU1 , ' (A1 ) ' )  TEMCHAR 
C 

ICD  =  0 

CALL  HOMCLS ( ICD ) 

IF ( IOU6 . NE . IOU1 )  WRITE (IOU6, ' ( ' *1' 1 ) ' ) 

C 

IF ( IOU6 . EQ . 6 )  WRITE ( IOU6 ,93)  COMM ( 1 ) , COMM ( 2 ) 

C 

WRITE ( IOU6 , 1007 )  RV20UT 

1007  FORMAT (  12X, 'SINGLE  ROTOR  FORCE  PROFILE  AT  RADIUS  = 

1  F7.1, 1  FT' ,//) 

C 

WRITE (IOU6, 1008) 

1008  FORMAT (  2X HEIGHT ', 6 X PEAK  Q  ' , 6 X ,  ’ FOVER ' , 7 X ,  ’ OVERM ' , 7 X , 

1  'TOT  F’ ,7X, 'TOT  M ' , / , 

2  3X, ' (FT) ’ ,8X, ' ( PSF ) * ,7X, ' (LB) \6X,  ' (FT-LB) '  ,7X, 

3  '  (LB)  '  ,6X,  '  (FT-LB)  '  ,/) 

C 

520  CONTINUE 
C 

C  - 

C  CALL  SUBROUTINE  TO  CALCULATE  THE 

C  FORCES  AND  MOMENTS  ON  A  HUMAN  BEING 

C  - 

C 

CALL  MOMENT ( NPTS , HUMTYP , TOTF , TOTM) 

C 

IF ( DSET . EQ . 0 . ) GOTO  545 
C 

HH  =  RVZ  * RADIUS 
HHOUT  =  HH  +  DXO 
C 

WRITE ( IOU6, 1014)  HHOUT , TOTF , TOTM 
1014  FORMAT (  1  8X  ,  F8 . 2 , 2F1 0 . 3 ) 

C 

IF ( IGRAPH . EQ . 1 ) THEN 
C 

WRITE ( IOU8 , 90 )  HHOUT , TOTF , TOTM 
90  FORMAT (  IX , F7 . 2 , 1 X , F7 . 2 , IX, F8 . 2 ) 

C 

ENDIF 

C 

545  CONTINUE 
C 


RVZ  =  RVZ  +  DELH 
C 

565  CONTINUE 
C 
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325  C  - 

326  C  CLOSE  AN  OPEN  GRAPHICS  FILE 

327  C  - 

328  C 

329  IF(IGRAPH.EQ.1)THEN 

330  C 

331  CLOSE ( IOU8 , STATUS= ' KEEP* ) 

332  C 

333  END IF 

334  C 

335  CALL  INKEY 

336  C 

337  GOTO  999 

338  C 

339  C  - 

340  C  THE  ERROR  LOGIC  ALLOWS  FOR  THE  HANDLING  OF  FILE 

341  C  OPEN  ERRORS  BY  RETURNING  THE  USER  TO  A  MENU 

342  C  - 

343  C 

344  2000  CONTINUE 

345  C 

346  CALL  HOMCLS ( 0 ) 

347  WRITE ( IOU 1,2001) 

348  2001  FORMAT (  ////,8X, 

349  1  '  ***  ERROR  ***  PLEASE  CHOOSE  A  NEW  OUTPUT  FILENAME’, 

350  2  /// , 8X, '  TYPE  <RETURN>  TO  CONTINUE  ’,$) 

351  READ ( IOU1 , ' ( A1 ) ' )  TEMCHAR 

352  KEY  =  'P' 

353  C 

354  999  CONTINUE 

355  C 

356  RETURN 

357  END 

358  C 

359 
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1  c 

2  C 

3  SUBROUTINE  INKEY 

4  C 

5  C  ★  *  ★★★*★★★★★★★*★*★*★★★★★★★★★  ****-#  it 

6  C  SUBROUTINE  INKEY 

7  <3  ****************************************************** 

8  C 

9  PARAMETER ( NUM  =  8) 

10  C 

11  CHARACTER* 1  KEY , KKEY 

12  CHARACTER*!  OKLIST (NUM) 

13  C 

14  COMMON  /  CKEY /  KEY, KKEY 

15  COMMON  /  UNITS/  IOU1 , IOU4 , IOU5 , IOU6 , IOU7 , IOU8 , IGRAPH 

16  C 

17  c  ****************************************************** 

18  C 

19  DATA  OKLIST  / ' C '  ,  ' c ’  ,  ' P 1  ,  ’ p '  ,  ' N ' , ' n ' ,  *  X '  ,  *  X  ’  / 

20  C 

21  c  - 

22  C  INQUIRE,  OBTAIN,  AND  CHECK  FOR  VALID  MENU  OPTION 

23  C  - 

24  C 

25  10  CONTINUE 

26  C 

27  WRITE (IOUl, 20) 

28  20  FORMAT (  //,7X,'  TYPE  <C>ONTINUE,  NEXT  <P>OINT,  <N>EW  CASE,' 

29  1  , '  E<X>IT  ==>  ' ,$) 

30  READ ( I OU 1 , ' ( A 1 ) ' )  KEY 

31  C 

32  IF (LEGAL (KEY, IOUl , OKLIST, NUM) .EQ. 1 ) GOTO  10 

33  C 

34  c  - 

35  C  CORRECT  LOWER  CASE  LETTERS  TO  UPPER  CASE 

36  C  TO  USE  AS  VALID  FLAGS  IN  PARENT  SUBROUTINE 

37  c  - 

38  C 

39  IF ( KEY . EQ . '  c  ' )  KEY  =  'C' 

40  IF ( KEY . EQ • ’ p ' )  KEY  =  'P' 

41  IF ( KEY . EQ . 1 n ' )  KEY  =  'N' 

42  IF ( KEY . EQ . ' X ' )  KEY  =  'X' 

43  C 

44  c  - 

45  C  CLEAR  SCREEN  AND  HOME  CURSOR 

46  C  - 

47  C 

48  ICD  =  0 

49  CALL  HOMCLS ( ICD ) 

50  C 

5 1  RETURN 

52  END 

53  C 

54 
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27 

28 

29 

30 

31 

32 

33 

34 
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C 

c 

SUBROUTINE  INPUT 
C 

c  ******************************************************** 

C  SUBROUTINE  INPUT 

C 

C  THIS  SUBROUTINE  PRESENTS  THE  INPUT  STATUS  MENU 

C  AND  MANIPULATES  THE  DATA  FOR  PROGRAM  USE 

£  ♦it************'************************'****************** 

C 

PARAMETER (NUM1  =  19) 

PARAMETER (NUM2  =  9) 

PARAMETER (NUM3  =  4) 

PARAMETER (NUM4  =  4) 

C 

CHARACTER* 1  CHDOL 
CHARACTER* 1  CENTRY 
CHARACTER* 1  OKLSTl(NUMl) 

CHARACTER* 1  OKLST2(NUM2) 

CHARACTER* 1  OKLST3(NUM3) 

CHARACTER* 1  OKLST4(NUM4) 

CHARACTER* 1  ICONT ( 5 ) 

C 

CHARACTER* 12  PTSFIL(4) 

C 

CHARACTER* 50  COMM ( 2 ) , LENTRY 
CHARACTER* 50  PROMPT 
C 

DIMENSION  CONT<9) ,CONTV(7) ,CONTX(8) 

C 

COMMON  /INPUTC/  ICONT , COMM, PTSFIL 
COMMON  /INPUTD/  CONT , CONTV , CONTX 

COMMON  /  UNITS/  IOU1 , IOU4 , IOU5 , IOU6 , IOU7 , IOU8 , IGRAPH 
C 

Q  ******************************************************** 

C 

c  - 

C  SET  DATA  TO  CHECK  FOR  ILLEGAL  DATA  INPUT 

C  - 

C 

DATA  OKLST1  /'  ‘,'A',’a’,'B','b‘,'C*,'c',’D’,'d’, 

1  E|6ff,f|0»9»H,h,I|i/ 

DATA  OKLST2  /'  1 , ’ A1 , ' a' , ' B’ , 1  b ' , ' C * , ' c’ , '  D ' , ' d' / 

DATA  OKLST3  /,V,/,v,,,H',,h'/ 

DATA  OKLST4  /'Y',,y','N,r,n'/ 

C 

10  CONTINUE 
C 

C  - 

C  CLEAR  SCREEN  AND  HOME  CURSOR 

C  - 

c 

CALL  HOMCLS ( 0 ) 


INPUT-1 


SUBROUTINE  INPUT 


55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 


CALL  LOCATE (2,1) 

C 

c  - 

C  WRITE  FIRST  ENGINEERING  DATA  MENU 

C  - 

C 

IROTOR  =  IFIX(CONT ( 1 ) ) 

C 

WRITE ( IOU1 , ' ( 24X / A/ ) ' )  '  ROTWASH  USER  INPUT  DATA  MENU' 

WRITE ( IOU1 , 1 (T7,A,T25,A,T51,A,T61,A/) ' )  '  CODE ’ , *  PARAMETER', 
1  ’  VALUE','  UNITS' 

WRITE (IOU1, ' (8X,A,T50,I5,8X,A) ' ) 

1  '  A  NUMBER  OF  ROTORS  (1  OR  2)  ', IROTOR, ' -ND- ' 

WRITE ( IOUl , ' (8X,A,T50,F8.1,6X,A) ’ ) 

1  '  B  HUB  TO  HUB  ROTOR  SEPARATION ', CONT ( 2 ),’ FT ' 

WRITE (IOUl, ' (8X,A,T50,F8.1,6X,A) ' ) 

1  '  C  ROTOR  RADIUS  ' , CONT ( 3 ) , ' FT ’ 

WRITE (IOUl, ’ (8X,A,T50,F8.1,6X,A) ' ) 

1  '  D  GROSS  WEIGHT  ' , CONT ( 4 ) , ' LB ’ 

WRITE (IOUl, ' (8X,A,T50,F8.1,6X,A) ' ) 

1  '  E  FUSELAGE  DOWNLOAD  FACTOR  ' , CONT ( 5 ) , ' PCT ' 

WRITE (IOUl, ' (8X,A,T50,F8.1,6X,A) ' ) 

1  ’  F  ROTOR  HEIGHT  ABOVE  GROUND  ' , CONT ( 6 ) , ' FT ' 

WRITE (IOUl, ' (8X,A,T50,F8.1,6X,A) ' ) 

1  '  G  SHAFT  TILT  ANGLE  (<20  DEG)  ’ , CONT ( 7 ) , ' DEG ’ 

WRITE (IOUl, ' (8X,A,T50,F8.4,6X,A) ' ) 

1  '  H  AIR  DENSITY  RATIO  ,CONT ( 8 ) , ’ ND ' 

WRITE (IOUl, ' (8X,A,T50,F8.1,6X,A//) ' ) 

1  '  I  AMBIENT  WIND  (-10  TO  10  KT ) ' , CONT ( 9 ) , ' KT ' 

C 

C  - 

C  PROMPT  FOR,  OBTAIN,  AND  CHECK  FOR  LEGAL  INPUT  DATA 

C  - 

C 

20  WRITE ( IOUl , ' (8X,A,$) ' ) 

1 '  ENTER  DATA  ENTRY  CODE  OR  <RETURN>  TO  CONTINUE  ==>  ' 

C 

READ (IOUl, ' (Al) ' )  CHDOL 
C 

IF (LEGAL (CHDOL, IOUl ,OKLSTl ,NUM1 ) .EQ. 1 )GOTO  20 
C 

IF ( CHDOL . EQ .  ')GOTO  30 
C 

IF ( CHDOL. EQ.  'A'  . OR . CHDOL . EQ .  'a  1  ) THEN 
C 

PROMPT  =  'ROTORS  =  ' 

CALL  IREAD( IOUl , PROMPT, IROTOR) 

C 

IF ( IROTOR.LT. 1)  IROTOR  =  1 
IF( IROTOR. GT. 2 )  IROTOR  =  2 
CONT ( 1 )  =  FLOAT ( IROTOR) 

C 

GOTO  19 
C 


INPUT-2 


SUBROUTINE  INPUT 


109 

110  C 

111 

112  C 

113 

114 

115 

116  C 

117 

118  C 

119 

120  C 

121 
122 

123 

124  C 

125 

126  C 

127 

128  C 

129 

130 

131 

132  C 

133 

134  C 

135 

136  C 

137 

138 

139 

140  C 

141 

142  C 

143 

144  C 

145 

146 

147 

148  C 

149 

150  C 

151 

152  C 

153 

154 

155  C 

156 

157 

158 

159  C 

160 

161  C 

162 


ENDIF 

IF (CHDOL . EQ . 'B' . OR . CHDOL . EQ . 'b' ) THEN 

PROMPT  =  ‘ROTOR  SEPARATION  =  * 

CALL  FREAD(IOUl , PROMPT, CONT( 2) ,1.0) 
GOTO  10 

ENDIF 

IF (CHDOL . EQ . ' C ' . OR .CHDOL . EQ . ' c ' ) THEN 

PROMPT  =  ‘ROTOR  RADIUS  =  ’ 

CALL  FREAD ( IOU1 , PROMPT, CONT( 3) ,1.0) 
GOTO  10 

ENDIF 

IF(CHDOL.EQ. 'D' . OR . CHDOL . EQ . 'd‘ (THEN 

PROMPT  =  ‘GROSS  WEIGHT  =  ’ 

CALL  FREAD (IOU1, PROMPT, CONT( 4), 1.0) 
GOTO  10 

ENDIF 

IF ( CHDOL. EQ. ‘E‘ . OR . CHDOL . EQ . ‘e ’ ) THEN 

PROMPT  =  'DOWNLOAD  FACTOR  =  ’ 

CALL  FREAD  (  IOU  1 ,  PROMPT, CONT'(  5)  ,1.0) 
GOTO  10 

ENDIF 

IF(CHDOL . EQ. ‘F’ . OR . CHDOL . EQ . 'f ’ ) THEN 

PROMPT  =  'HEIGHT  ABOVE  GROUND  =  ’ 
CALL  FREAD ( IOU 1 , PROMPT , CONT ( 6 ) ,1.0) 
GOTO  10 

ENDIF 

IF(CHDOL.EQ. 'G‘ .OR .CHDOL . EQ . ’g‘ (THEN 

PROMPT  =  'SHAFT  TILT  =  ' 

CALL  FREAD ( IOU 1 , PROMPT, CONT ( 7) ,1.0) 

IF(CONT( 7 ) .LT.  0.0)  CONT<7)  =  0.0 

I F ( CONT ( 7 ) . GT .20.0)  CONT(7)  =  20.0 
GOTO  10 

ENDIF 

IF (CHDOL. EQ. 'H' . OR . CHDOL . EQ . 'h‘ ) THEN 
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prompt  =  'DENSITY  RATIO  =  ' 

CALL  FREAD ( IOU 1 , PROMPT , CONT ( 8 ) ,1.0) 
GOTO  10 

ENDIF 

IF { CHDOL . EQ. ' I ' . OR . CHDOL. EQ. ’ i ' ) ThEN 

PROMPT  =  'WIND  VELOCITY  =  ' 

CALL  FREAl>(  IOU1  ,  PROMPT , CONT ( 9  )  ,  1 . 0) 

IF (CONT ( 9 ) . LT .-10.0)  CONT ( 9 )  =  -10.0 
IF (CONT ( 9 )  .GT .  1 0 . 0 )  CONT(9)  =  10.0 

GOTO  10 

ENDIF 

GOTO  10 

30  CONTINUE 


CLEAR  SCREEN  AND  HOME  CURSOR 


CALL  HOMCLS(O) 
CALL  LOCATE (2,1) 


WRITr  SECOND  ENGINEERING  DATA  MENU 


WRITE (ICU1, 40) 

40  FORMAT (  '  ',18X,'  RCTWASH  PROGRAM  LOGIC /COMMENT  MENU',//, 

1  10X, 'CODE  PARAMETER  VALUE'/) 

WRITE (IOU1 ,50)  ICONT ( 1 ) , ICONT ( 2 ), COMM ( 1 ) ,COMM(2) 

50  FORMAT (  12X,'A  ANALYSIS  TYPE,  <V>  OR  <H> ' , 5X , A4 , / , 

1  1  2  X ,  '  B  GRAPHICS  FILE,  <Y>  OR  <N>  '  ,  5X ,  A4  ill , 

^  15X,  USER  INPUT  COMMENTS  (FOR  "PRN”  AND  "PLT”  OUTPUT) ’ , 

3  //  ,  1 7X,  ’  < -  50  SPACES  _ >•, 

4  /, 12X, 'C' , 4X , A50 , /, 

5  12X,  'D  '  , 4 X , A 5 0 , / / ) 


PROMPT  FOR,  OBTAIN,  AND  CHECK  FOR  LEGAL  INPUT  DATA 


60  WRITE ( IOU 1 , ' ( 8X , A, $ ) ' ) 

1  ENTER  CODE  FOR  DATA  INPUT  OR  <RETUkN>  TO  CONTINUE 
READ ( IOU1 , ' ( A1 ) ‘ )  CHDOL 
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237 
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253 

254 

255 

256 

257 

258 
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260 
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263 

264 

265 

266 

267 

268 
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C 

C 

c 

c 

c 

c 

c 

c 


c 

c 


c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 

c 

c 

c 

c 

c 

c 


IF (LEGAL (CHDOL , IOU1 , OKLST2 , NUM2 ) .EQ.l)GOTO  60 
IF ( CHDOL. EQ. 1  1 )GOTO  70 

CHOOSE  VELOCITY  OR  HAZARD  ANALYSIS  OPTION 


IF (CHDOL. EQ. 'A' .OR. CHDOL. EQ . ’a ’ ) THEN 
80  CONTINUE 

WRITE ( IOU1 ,  1  (  /  , 35X , A , IX , A2 / ) ' )  '  ANALYSIS  TYPE  = 

WRITE ( IOU1 , ' (40X/A,$) * )  '  ENTER  NEW  VALUE  ==> 
READ ( IOU1 , ' ( A1 ) ' )  CENTRY 

IF (LEGAL (CENTRY , IOU1 , OKLST3 ,NUM3 )  .EQ.l)GOTO  80 
ICONT(l)  =  CENTRY 

IF ( ICONT ( 1 ) . EQ . ’ v ' )  ICONT(l)  =  'V 
IF( ICONT ( 1 ) . EQ . ' h ' )  ICONT ( 1 )  =  ’H' 

GOTO  30 

ENDIF 


CHOOSE  OUTPUT  TO  A  GRAPHICS  FILE  OR  NOT 


IF (CHDOL . EQ . 'B' . OR . CHDOL . EQ . 'b’ ) THEN 
90  CONTINUE 


WRITE ( IOU1 , ' (/,35X,A,1X,A2/) ' )  '  GRAPHICS  FLAG  = 

WRITE ( 10U1 , ' ( 4  OX , A , $ ) ' )  ’  ENTER  NEW  VALUE  ==> 

READ ( IOU1 , ' ( A1 ) ' )  CENTRY 

IF ( LEGAL ( CENTRY, IOU1 ,OKLST4 , NUM4 ) .EQ. 1 )GOTO  90 


ICONT ( 2 )  =  CENTRY 


IF ( ICONT ( 2 ) .EQ. 
IF ( ICONT ( 2 ) .EQ. 
IF ( ICONT ( 2 ) .EQ. 
IF ( ICONT ( 2 ) .EQ. 
GOTO  30 


'y')  ICONT ( 2 )  = 
’n')  ICONT ( 2 )  = 
'  Y'  )  IGRAPH  =  1 
'  N ’  )  IGRAPH  =  0 


1  Y' 
’N  ' 


ENDIF 


CHOOSE  COMMENT  STRINGS  FOR  OUTPUT  DATA 


IF ( CHDOL. EQ. 'C' . OR . CHDOL . EQ . 'c ’ ) THEN 


'  , ICONT ( 1 ) 


'  , ICONT ( 2 ) 
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WRITE ( IOU1 , ' ( /, 5X,A,A50/ ) ' )  '  COMMENT  STRING  = 
WRITE ( IOU1 , ' { 5X , A, $ ) ' )  '  ENTER  STRING  ==>  ' 
READ( IOU1 , 1 ( A50 ) ' )  LENTRY 
C 

COMM ( 1 )  =  LENTRY 
GOTO  30 
C 

ENDIF 

C 

IF { CHDOL . EQ . ' D ’ . OR . CHDOL . EQ . ' d ’ ) THEN 
C 

WRITE(IOUl, ' (/ ,5X,A,A50/) 1 )  1  COMMENT  STRING  = 
WRITE ( IOU1 , ' ( 5X , A, $ ) 1 )  ’  ENTER  STRING  ==>  ' 
READ ( IOU1 , ’ ( A50 ) ' )  LENTRY 
C 

COMM ( 2 )  =  LENTRY 
GOTO  30 
C 

ENDIF 

C 

GOTO  30 
C 

70  CONTINUE 
C 


RETURN 

END 


'  , COMM( 1 ) 


'  , COMM ( 2 ) 
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54 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 


c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  INPUTV (FLOW) 

★  ★★★★★★★★★★♦★★★★♦★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★•AT** 

SUBROUTINE  INPUTV 

THIS  SUBROUTINE  PRESENTS  THE  INPUT  STATUS  MENU  AND  MANIPULATES 
DATA  FOR  FOR  THE  WALJET  AND  I PLANE  PROGRAM  OPTIONS 

★  ★♦♦★★★★★★★★★★★★★♦★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★I** 


PARAMETER (NUN  =  9) 

CHARACTER* 1  OKLIST(NUM) 

CHARACTER* 1  CHDOL,FLOW 
CHARACTER* 1  ICONT(5) 

CHARACTER* 12  TMPFIL 
CHARACTER* 12  PTSFIL(4) 

CHARACTER* 50  COMM(2) 

CHARACTER* 50  PROMPT 

DIMENSION  CONT ( 9 ) ,CONTV(7) ,CONTX(8) 

COMMON  /INPUTC/  ICONT , COMM ,PTSFIL 
COMMON  /INPUTD/  CONT , CONTV ,CONTX 

COMMON  /  UNITS/  IOU1 , IOU4 , IOU5 , IOU6 , IOU7 , IOU8 , IGRAPH 


DATA  OKLIST  /'  * , ' A ' , ' a ’ , ’ B ’ , ’ b ’ , ’ C 1 , ' C ’ , ' D ' , ' d ' / 


CLEAR  SCREEN  AND  HOME  CURSOR 


ICD  =  0 

20  CONTINUE 

CALL  HOMCLS ( ICD ) 

CALL  LOCATE (1,1) 

WRITE (IOU1, 12) 

12  FORMAT (  22X, '  VELOCITY  PROFILE  STATUS  MENU’,///, 

1  1  OX,  ’CODE  PARAMETER  VALUE’, 

2  3X, 'UNITS ', /) 


PRINT  OUT  MENU  VARIABLES  AS  BASED  ON  THE  WALL  JET 
OPTION  OR  INTERACTION  PLANE  OPTION  SWITCH  SETTING 


IF (FLOW . EQ. 'W' ) THEN 
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55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 


C 

c 

c 

c 


c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


c 


c 

c 

c 


c 


c 

c 

c 


c 


WRITE (IOU1, 14) (CONTV(I) ,1=1,3) ,PTSFIL(1) 
ELSE 


WRITE (IOU1, 14) (CONTV(I) ,1=1,3) ,PTSFIL(2) 
ENDIF 


14  FORMAT (  12X,  'A 

1  4X, ' FT ' , / , 

2  12X, ' B 

3  12X, 'C 

4  12X,  '  D 


PROFILE  STATION  POSITION  ',5X,F7.2, 


VERTICAL  INCREMENT 
MAXIMUM  PROFILE  HEIGHT 
DATA  OUTPUT  FILENAME 


'  ,5X,F7.2,4X, 'FT' ,/, 

1 ,5X,F7.2,4X, ’FT* ,//, 
'  , 7X, A12 , // ) 


PROMPT  FOR,  OBTAIN,  AND  CHECK  FOR  LEGAL  INPUT  DATA 


10  CONTINUE 

WRITE (IOU1, 1 ( 8X, A, $ ) * ) 

1  '  ENTER  DATA  ENTRY  CODE  OR  <RETURN>  TO  CONTINUE  ==> 

READ ( IOU1 , ' ( A1 ) ' )  CHDOL 

IF ( LEGAL ( CHDOL , IOU1 , OKLIST , NUM) .EQ.l)GOTO  10 

IF(CHDOL.EQ.  ' )GOTO  30 

IF ( CHDOL . EQ . ' A ' . OR . CHDOL .  EQ . ' a ' ) THEN 

PROMPT  =  'PROFILE  STATION  POSITION  =  ’ 

CALL  FREAD ( I OU 1 , PROMPT , CONTV ( 1 ) ,1.0) 

IF (CONTV ( 1 ) . LT .0.0)  CONTV(l)  =  0.0 
GOTO  20 

ENDIF 

IF ( CHDOL . EQ . ' B ' . OR . CHDOL . EQ . ’ b ' ) THEN 

PROMPT  =  'VERTICAL  INCREMENT  =  ' 

CALL  FREAD ( IOU1 , PROMPT , CONTV ( 2 ) , 1 . 0 ) 

IF ( CONTV ( 2 ) . LT .0.0)  CONTV(2)  =  0.0 
GOTO  20 

ENDIF 

IF ( CHDOL. EQ. 'C' . OR .CHDOL . EQ . 'c' ) THEN 

PROMPT  =  'MAXIMUM  PROFILE  HEIGHT  =  ' 

CALL  FREAD ( IOU1 , PROMPT , CONTV ( 3 ) , 1 . 0) 
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IF ( CONTV ( 3 ) . LT .0.0)  CONTV(3)  =  0. 

110 

GOTO  20 

111 

C 

112 

ENDIF 

113 

C 

114 

C 

115 

C 

CHOOSE  GRAPHICS  FILENAME 

116 

C 

117 

C 

118 

IF(CHDOL.EQ. ‘D‘ . OR . CHDOL . EQ . ‘d1 ) THEN 

119 

C 

120 

IF ( FLOW . EQ . 'W' ) THEN 

121 

WRITE (IOU1, ’ (/,25X,A,1X,A12/) ’ ) 

122 

1  1  FILENAME  =  ’.PTSFILfl) 

123 

ELSE 

124 

WRITE (IOU1, 1 (/,25X,A,1X,A12/) ’ ) 

125 

1  '  FILENAME  =  ’,PTSFIL(2) 

126 

ENDIF 

127 

C 

128 

WRITE (IOU1, ’ (20X,A,$) ' ) 

129 

1  '  ENTER  NEW  FILENAME  ( XXXXXXXX . XXX) 

130 

c 

131 

READ ( IOU1 , ' ( A1 2 ) ' )  TMPFIL 

132 

c 

133 

IF (FLOW. EQ. *W )  PTSFIL ( 1 )  =  TMPFIL 

134 

IF(FLOW. EQ . ' I ' )  PTSFIL( 2 )  =  TMPFIL 

135 

GOTO  20 

136 

c 

137 

ENDIF 

138 

c 

139 

GOTO  20 

140 

c 

141 

30  CONTINUE 

142 

c 

143 

ICD  =  0 

144 

CALL  HOMCLS ( ICD ) 

145 

c 

146 

RETURN 

147 

END 

148 

c 

149 

INPUTV 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


c 

c 


c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


c 


c 


SUBROUTINE  INPUTX 

***************************************************** 
SUBROUTINE  INPUTX  PRESENTS  THE  INPUT  STATUS  MENU  AND 
MANIPULATES  DATA  FOR  PROGRAM  USE  WITH  THE  GROUND  AND 
DISC  VORTEX  OPTIONS 


PARAMETER ( NUM  =  15) 

CHARt*<~TER*  1  OKLIST  (NUM ) 

CHARACTER* 1  CHDOL 
CHARACTER* 1  ICONT{5) 

CHARACTER* 12  PTSFIL(4) 

CHARACTER* 50  COMM<2) 

CHARACTER* 50  PROMPT 

DIMENSION  CONT ( 9 ) ,CONTV(7) ,CONTX(8) 

COMMON  / INPUTC/  ICONT , COMM , PTSFIL 
COMMON  /INPUTD/  CONT ,CONTV ,CONTX 

COMMON  /  UNITS/  IOU1 , IOU4 , IOU5 , IOU6 , IOU7 , IOU8 , IGRAPH 


DATA  OKLIST  /'  •,,A','a',,B,,,b',,C','c,,,D','<i', 

1  'E' , 'e' , 'F' , 'f ' , 'G' , 'g'/ 


CLEAR  SCREEN  AND  HOME  CURSOR 


ICD  =  0 
20  CONTINUE 


CALL  HOMCLS ( ICD ) 
CALL  LOCATE (1,1) 


WRITE (IOU1 ,10) 

10  FORMAT (  1 8X, '  GROUND/DISK 

1  15X, '  (FOR  SINGLE 

1  8X , 'CODE 

2  4X, 'UNITS' ,/) 


VORTEX  INPUT  DATA  MENU',/, 

MAIN  ROTOR  HELICOPTERS  ONLY)',///, 
PARAMETER  VALUE ’ , 


WRITE (IOU1, 11) (CONTX(I) ,1=1,7) 

11  FORMAT (  10X, 'A  ROTOR  TIP  SPEED 

1  4X , ' FPS ' , / , 

2  1  OX ,  ' B  NUMBER  OF  ROTOR  BLADES 

3  10X,'C  TRANSLATIONAL  SPEED 

4  10X, ' D  XT  POSITION 

5  10X, ' E  YT  POSITION 


’  ,5X,F8.2, 


'  ,5X,F8.2,4X,  '-ND-' ,/, 
’  ,5X,F8.2,4X,  'FPS' ,/, 

'  , 5X,F8 . 2 , 4X,  'FT' ,/, 

'  ,5X,F8.2,4X,  'FT'  ,/, 
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56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 


6  1  OX ,  ‘ F  ZT  CALCULATION  INCREMENT  ,5X,F8.2,4X,  ' FT ' , / , 

7  1  OX ,  ' G  MAXIMUM  CALCULATION  HEIGHT' ,5X,F8. 2, 4X, 'FT' ,///) 

C 

c  - 

C  INQUIRE,  OBTAIN,  AND  CHECK  FOR  VALID  MENU  OPTION 

C  - 

c 

40  CONTINUE 
C 

WRITE ( IOU1 , ' ( 8X,A, $ ) ' ) 

1  ’  ENTER  DATA  ENTRY  CODE  OR  <RETURN>  TO  CONTINUE  ==>  ‘ 

C 

READ ( IOU1 , ‘ { A1 ) ' )  CHDOL 
C 

IF (LEGAL (CHDOL, IOU1 , OKLIST, NUM) .EQ.lJGOTO  40 
C 

IF(CHDOL.EQ.  ')GOTO  30 
C 

IF ( CHDOL . EQ . ' A ' . OR , CHDOL . EQ . ' a ' ) THEN 
C 

PROMPT  =  ‘ROTOR  TIP  SPEED  =  ' 

CALL  FREAD ( I OU1, PROMPT ,CONTX(l) ,1.0) 

C 

IF (CONTX ( 1 ) . LT .0.0)  CONTX(l)  =  0.0 
GOTO  20 
C 

ENDIF 

C 

IF ( CHDOL. EQ. 'B' . OR . CHDOL . EQ . 'b' ) THEN 
C 

PROMPT  =  'NUMBER  OF  ROTOR  BLADES  =  ' 

CALL  FREAD(IOUl, PROMPT, CONTX ( 2 ) ,1.0) 

C 

IF ( CONTX ( 2 ) . LT .2.0)  CONTX(2)  =2.0 
GOTO  20 
C 

ENDIF 

C 

IF (CHDOL. EQ. 'C' . OR . CHDOL . EQ . 'c' ) THEN 
C 

PROMPT  =  'TRANSLATIONAL  SPEED  =  ’ 

CALL  FREAD ( IOU 1 , PROMPT , CONTX ( 3 ) ,1.0) 

C 

IF ( CONTX ( 3 ) . LT .0.0)  CONTX ( 3 )  =0.0 
GOTO  20 
C 

ENDIF 

C 

IF(CHDOL . EQ . 'D' . OR . CHDOL . EQ .  d' ) THEN 
C 

PROMPT  =  'XT  POSITION  =  ' 

CALL  FRE AD ( IOU 1 , PROMPT, CONTX ( 4 ) , 1 . 0 ) 

GOTO  20 
C 
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ENDIF 

110 

C 

111 

IF(CHDOL.EQ. 'E' .OR.CHDOL.EQ. 'e' ) THEN 

112 

C 

113 

PROMPT  =  'YT  POSITION  =  * 

114 

CALL  FREAD { IOU1 , PROMPT , CONTX ( 5 ) , 1 . 0 ) 

115 

GOTO  20 

116 

C 

117 

ENDIF 

118 

C 

119 

IF (CHDOL. EQ. 'F‘ .OR.CHDOL.EQ.  ' f  * ) THEN 

120 

C 

121 

PROMPT  =  *ZT  CALCULATION  INCREMENT  = 

122 

CALL  FREAD ( IOU 1 , PROMPT , CONTX ( 6 ) ,1.0) 

123 

C 

124 

IF  ( CONTX ( 6 ) . LT . 0 . 0 )  CONTX(6)  =  0.0 

125 

GOTO  20 

126 

C 

127 

ENDIF 

128 

C 

129 

IF ( CHDOL . EQ . ' G ’ . OR . CHDOL . EQ . ' g ‘ ) THEN 

130 

C 

131 

PROMPT  =*  'MAXIMUM  CALCULATION  HEIGHT 

132 

CALL  FREAD { I OU 1 , PROMPT , CONTX ( 7 ) ,1.0) 

133 

C 

134 

IF (CONTX ( 7 )  . LT .0.0)  CONTX(7)  =  0.0 

135 

GOTO  20 

136 

C 

137 

ENDIF 

138 

C 

139 

GOTO  20 

140 

c 

141 

30  CONTINUE 

142 

c 

143 

ICD  =  0 

144 

CALL  HOMCLS ( ICD ) 

145 

c 

146 

RETURN 

147 

END 

148 

149 

c 
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1  C 

2  C 

3 

4  C 

5  C 

6  C 

7  C 

8  C 

9 

10  C 

11  C 

12  C 

13  C 

14  C 

15  C 

16  C 

17 

18 

19  C 

20 
21 

22  C 


SUBROUTINE  IOFNSH 

*************************************************** 
SUBROUTINE  IOFNSH  CLOSES  FILES  OPENED  FOR  DISK  I/O 

COMMON  /  UNITS/  IOU1 , IOU4 , IOU5 , IOU6 , IOU7 , IOU8 , IGRAPH 

♦inn************************************************ 


KEEP  STANDARD  I/O  FILES  IF  ON  DISK,  ELSE  DELETE 


IF( IOU5 .EQ. 5 )  CLOSE ( IOU5 ,STATUS= 'KEEP  1 ) 
IF ( IOU6 . EQ . 6 )  CLOSE ( IOU6 , STATUS= ' KEEP’ ) 

RETURN 

END 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 


C 

C 

SUBROUTINE  IOINIT 
C 

C  SUBROUTINE  IOINIT  DISPLAYS  THE  OPENING  BANNER 

C  AND  OPENS  THE  FILES  FOR  DISK  I/O,  FILENAMES 

C  ARE  PROMPTED  FROM  THE  TERMINAL 
£  ★  **★*★★★★**★★*★**★*★**•*■★★★**★*★•**★★★★★★★★*★★★*★★★★★★* 

C 

CHARACTER* 3  IPFILE , OPFILE 
CHARACTER* 1  TEMCHAR 
C 

COMMON  /  UNITS/  IOU1 , IOU4 , IOU5 , IOU6 , IOU7 , IOU8 , IGRAPH 
C 

c 

c  - 

C  ASSIGN  DEFAULT  VALUES  TO  I/O  UNIT  POINTERS 

C  - 

c 

IOU1  =  0 

IOU4  =  0 

IOU5  =  0 

IOU6  =  6 

IOU7  =  0 

IOU8  =  8 

IGRAPH  =  0 

C 

C  - 

C  HOME  CURSOR  AND  CLEAR  THE  DISPLAY 

C  - 

C 

ICD  =  0 

CALL  HOMCLS ( ICD ) 

C 

C  - 

C  DISPLAY  BANNER 

C  - 

C 

CALL  LOCATE (5,1) 

C 

c  - 

C  ORIGINAL  ROTHAZ  PROGRAM  WAS  VERSION  1.0 
C  FIRST  VAX  MENU  VERSION  OF  ROTHAZ  WAS  VERSION  1.1 

C  PROGRAM  ROTWASH  REPLACES  ROTHAZ  AT  VERSION  2.0 

C  - 

C 

WRITE ( IOU1 , 10) 

10  FORMAT (  27X, 'ROTWASH  PROGRAM',// 

1  /, 17X, 'ROTORCRAFT  DOWNWASH  HAZARD  ANALYSIS',/ 

2  /, 15X, ' EMA/COMPUTATIONAL  METHODOLOGY  ASSOCIATES', 

3  / , 1 5X , ' *  *  *  PROGRAM  VERSION  2.0,  APRIL  1991  ***' 

4  ///////) 
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55  C 

56  WRITE (IOU1, ’ (A, $) ' ) 

57  1  '  PRESS  <RETURN> 

58  READ ( IOU1 , ' ( A1 ) ’ )  TEMCHAR 

59  C 

60  C  - 

61  C  HOME  CURSOR  AND  CLEAR  THE  DISPLAY 

62  C  - 

63  C 

64  ICD  =  0 

6  5  CALL  HOMCLS ( ICD ) 

66  CALL  LOCATE (3,1) 

67  C 

68  C  - 

69  C  PROMPT  FOR  I/O  FILES  AND  READ  USER  RESPONSE 

70  C  - 

71  C 

72  WRITE (IOU1, 12) 


73 

12  FORMAT ( 

14X,  ' 

I/O  CAN 

BE  DIRECTED  TO  FILES  OR  DEVICES 

74 

1 

19X,  ' 

VALID  DEVICES  ARE  AS  FOLLOWS:’,//, 

75 

2 

22X,  ' 

<CON> 

==>  CONSOLE',/, 

76 

2 

22X,  ' 

<PRN> 

==>  PRINTER',/, 

77 

4 

22X,  ' 

<PLT> 

==>  GRAPHICS  FILE ',///) 

78  C 

79  16  CONTINUE 

80  C 

81  WRITE ( IOU1 ,13) 

82  13  FORMAT (  16X, ’  ENTER  INPUT  FILE/DEV  NAME  ==>  ',$) 

83  11  FORMAT (A) 

84  READ ( IOU1 ,11)  IPFILE 

85  C 

86  C  - 

87  C  CHECK  FOR  DATA  ENTRY  ERROR 

88  C  - 

89  C 

90  IF ( IPFILE. NE. 'CON' . AND . IPFILE . NE . 'con' ) THEN 

91  WRITE ( IOU1 ,15) 

92  15  FORMAT ( / , 1 6 X ,  **  INPUT  ERROR,  PLEASE  REENTER  **',//) 

93  GOTO  16 

94  ENDIF 

95  C 

96  17  CONTINUE 

97  C 

98  WRITE ( IOU1 , 14) 

99  14  FORMAT (  16X,'  ENTER  OUTPUT  FILE/DEV  NAME  ==>  ',$) 

100  READ ( IOU1 ,11)  OPFILE 

101  C 

102  C  - 

103  C  CHECK  FOR  DATA  ENTRY  ERROR 

104  C  - 

105  C 

106  IF (OPFILE . NE . 'CON' .AND .OPFILE .NE . 'con' .AND. 

107  1  OPFILE. NE. 'PRN' .AND. OPFILE. NE. 'prn' .AND. 

108  2  OPFILE. NE. ' PLT ’ .AND. OPFILE. NE. 'pit' (THEN 
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109 

110 
111 

112  C 

113  C 

114  C 

115  C 

116  C 

117 

118 

119 

120  C 

121 
122 

123 

124  C 

125  C 

126  C 

127  C 

128  C 

129  C 

130  C 

131  C 

132 

133 

134  C 

135  C 

136  C 

137  C 

138  C 

139 

140 

141  C 

142 

143 

144  C 

145 


WRITE ( IOU1 ,15) 
GOTO  17 
END  IF 


REDIRECT  INPUT/OUTPUT  FILES  IF  REQUESTED 


IF ( IPFILE .EQ .  'CON  1 .OR . IPFILE . EQ.  ' con ‘ ) IOU5  =  IOU1 
IF(OPFILE . EQ. 'CON ' . OR . OPFILE . EQ . ’ con ' ) IOU6  =  IOU1 
IF ( OPFILE . EQ . ' PRN ' . OR . OPFILE . EQ . ' prn ' ) IOU6  =  6 

IF (OPFILE. EQ. 'PLT‘ .OR. OPFILE. EQ. 'pit' ) IGRAPH  =  1 
IF(OPFILE.EQ. 'PLT’ .OR. OPFILE. EQ. 'pit' ) IOU6  =  IOU1 
IF(OPFILE.EQ. 'PLT1 .OR. OPFILE. EQ. 'pit* ) IOU8  =  8 


OPEN  STANDARD  I/O  FILES 

IOU5  =  STD.  INPUT 
IOU6  *  STD.  OUTPUT 


IF ( IOU5 .EQ . 5 )  OPEN ( IOU5 , FILE=IPFILE, STATUS= ' OLD ’ ) 
IF ( IOU6 . EQ . 6 )  OPEN ( IOU6 , FILE=OPFILE , STATUS= ' NEW ' ) 


HOME  CURSOR  AND  CLEAR  THE  DISPLAY 


ICD  *  0 

CALL  HOMCLS ( ICD ) 

RETURN 

END 
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1  C 

2  C 

3 

4  C 

5  C 

6  C 

7  C 

8  C 

9  C 
10 

11 

12 

13 

14 

15  C 

16 

17 

18 

19  C 

20 
21 
22 

23 

24 

25  C 

26  C 

27  C 

28  C 

29  C 

30  C 

31  C 

32 

33  C 

34  C 

35  C 

36  C 

37  C 

38 

39  C 

40  C 

41  C 

42  C 

43  C 

44  C 

45 

46  C 

47 

48 

49 

50  C 

51 

52  C 

53 

54  C 


SUBROUTINE  IPVEL ( H , UN , RADIUS , UMB , XIP , YSEP , WSPD , DELZ , ZMAX, DXO) 


SUBROUTINE  IPVEL  GENERATES  THE  VELOCITY  PROFILE  V(X,Z)  AT 
XVZ  ALONG  THE  INTERACTION  PLANE  FOR  THE  TWO  ROTOR  CASE 


CHARACTER* 1  TEMCHAR 
CHARACTER* 1  ICONT(5) 

CHARACTER* 1  KEY, KKEY 
CHARACTER* 1 2  PTSFIL<4) 

CHARACTER* 50  COMM(2) 

DIMENSION  ZZ ( 40 ) ,VHMF(40) ,VHMK(40) ,VHPF(40) ,VHPK(40) , 

1  QHM (40)  ,  QHP  (40)  ,  WMF  (40)  ,WMK(40)  ,WPF(40)  , 

2  WPK  (40),  QVM  (40),  QVP  (40) 

COMMON  /  CKEY /  KEY, KKEY 

COMMON  /CONSTS/  PI , RHO , FPSPKN , RHOD2 , DRC 

COMMON  /INPUTC/  ICONT , COMM , PTSFIL 

COMMON  /PROFIL/  RJ , ZBJ , ZHJ , ZMJ , UMJ , ZB , ZH , ZM, UM,CU , CY 
COMMON  /  UNITS/  IOU1 , IOU4 , IOU5 , IOU6 , IOU7 , IOU8 , IGRAPH 


TF  IS  THE  INTERACTION  PLANE  AMPLIFICATION  FACTOR 


TF  =  1.55  -  (0.55)*EXP(-1.35*XIP) 


OBTAIN  PAPAMETERS  AT  BASE  RADIUS  FOR  THE  'BOUNDARY  LAYER' 


RIP0  a  SQRT(XIP**2  +  YSEP**2 ) 


’ PROPRM '  PROVIDES  THE  VELOCITY  PROFILE  PARAMETERS 
OF  A  RADIAL  WALL  JET  (WITHOUT  INTERACTION  PLANE) 


CALL  PROPRM ( H , UMB , RIP0 ) 

ZIPB  =  ZB 
Z I PM  =  ZM 
ZIPH  =  ZH 

RIPM  =  SQRT( XIP* *2  +  (YSEP  +  ZIPM)**2) 
CALL  PROPRM (H, UMB, RIPM) 
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55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 
7  5 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 


UMM  =  UM 
C 

c  - 

C  INCREMENT  AND  HEIGHT  ARE  DEL2  AND  ZMAX 

C  - 

c 

NPTS  =  IFIX ( ZMAX/DELZ )  +  2 
I F ( NPTS . GT .40)  NPTS  =40 
C 

c  - 

C  DIMENS IONALIZE  VELOCITY  PROFILE  PARAMETERS 

C  - 

c 

XXIP  =  RADIUS*XIP 
ZZB  =  ZIPB*RADIUS 
ZZH  =  ZI?H*RADIUS 
ZZM  =  ZIPM*RADIUS 
C 

c  - 

C  OUTPUT  THE  VELOCITY  AND  DYNAMIC  PRESSURE  PROFILE  HEADF° 

C  - - 

c 

ICD  =  0 

CALL  HOMCLS(ICD) 

C 

IF ( IOU6 .NE . IOU1 )  WRITE (IOU6,  '  (  1  '1 '  '  )  1  ) 

C 

IF ( IOU6 . EQ . 6 )  WRITE (IOU6, 93)  COMM ( 1 ) , COMM ( 2 ) 

93  FORMAT (  1  OX , A5 0 , / , 1 CX , A5 0 , / / ) 

C 

XIPOUT  =  XXIP  +  DXO 
C 

WRITE(IOU6, 1000)  XIPOUT 

1000  FORMAT (  2X,'TWIN  ROTOR  INTERACTION  PLANE  VELOCITY  PROFILE ' , 

1  *  AT  STATION  =  ' , F7 . 1 , *  FT',//) 

C 

WRITE (IOU6, 1002) 

1002  FORMAT (  2X ,  *  HE IGHT '  , 8 X ,  ' MEAN  VELOCITY 7X PEAK  VELOCITY 6X, 

1  'MEAN  Q ' , 4X , ' PEAK  Q ' , / , 

2  3X,  ' (FT)  ' ,8X,  '  (FPS)  , 6 X ,  ’ ( KN )  ' , 5X ,  ’ ( FPS )  ’ , 6X ,  ' ( KN )  ' , 5X , 

3  ’  (PSF)  '  ,5X,  '  (PSF)  '  ,/) 

C 

LINES  =  0 
C 

c  - 

C  'AN'  IS  ACTUALLY  '=  1.0/7. O' 

C  - 

c 

AN  =  0.142857142 
C 

c  - 

C  CALCULATE  THE  VELOCITY  PROFILE  POINTS  FOR  OUTPUT 

C  'NPTS'  IS  THE  NUMBER  OF  VERTICAL  STATION  POINTS 

C  - 
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109  C 

110 

111  C 

112 

113 

114  C 

115  C 

116  C 

117  C 

118  C 

119 

120  C 

121 

122  C 

123 

124 

125  C 

126  C 

127  C 

128  C 

129  C 

130 

131  C 

132  C 

133  C 

134  C 

135  C 

136 

137 

138  C 

139 

140  C 

141  C 

142  C 

143  C 

144  C 

145 

146 

147 

148 

149  C 

150  C 

151  C 

152  C 

153  C 

154 

155  C 

156 

157  C 

158 

159 

160 
161 

162  C 


DO  500  I  =  1 , NPTS 

LI.-^S  =  LINES  +  2 

ZIP  =  DELZ*FLOAT( I  -  1) 


GET  MAX  WALL  JET  VELOCITY  AT  EFFECTIVE  RADIUS 


RIP  =  SQRT(XIP**2  +  ( YSEP  +  ZIP)**2) 

CALL  PROPRM ( H , UMB , RIP ) 

VN  =  UN 
VZ  =  UM 

INTERACTION  PLANE  'BOUNDARY  LAYER' 

IF ( ZIP . LT . Z IPM)  VZ  =  UMM* ( ZIP/ZIPM) **AN 

DEVELOPED  INTERACTION  PLANE  JET 

VH  =  TF*VZ  *XIP/RIP 
W  =  TF*VZ* ( YSEP  +  ZIP) /RIP 

ZZ(I)  =  ZIP*RADIUS 

MEAN  VELOCITIES  AND  DYNAMIC  PRESSURE 


VHMF(I)  =  VH*UN 
WMF  ( I )  =  W*UN 
VHMK(I)  =  VHMF ( I ) /FPSPKN 
WMK(I)  =  WMF{  I ) /FPSPKN 


PEAK  VELOCITIES  AND  DYNAMIC  PRESSURE 


VMFD3I  =  ( XIP*  0 . 2  4  4  4  )  +  0.8 


I F ( VMFD  3 1 . GT .2.5)  VMFD3I  =  2.5 


VHPF ( I ) 
WPF(I ) 
VHPK ( I ) 
Wi>K  ( I ) 


VMFD  3 1 *VHMF ( I ) 
VMFD3I  *  WMF  (  I ) 
VHPF ( I ) /FPSPKN 
WPF  (  I  )  /FPSPKN 
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163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 
181 
182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 


C 

C 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


c 

c 


c 

c 

450 

c 

c 

c 

c 

c 


1003 

c 

c 

c 

c 

c 


IF ( VHMF ( I ) . EQ . 0 . ) GOTO  55 


THE  EFFECT  OF  WIND  IS  TO  ADD  (DOWNWIND  SIDE)  OR  SUBTRACT 
(UPWIND  SIDE)  ' XWK '  TIMES  THE  AMBIENT  WIND  VELOCITY  TO 
THE  HORIZONTAL  PROFILE  VELOCITY  (EMPIRICAL,  CH-53E  BASED) 


XKW  =  ( -0 . 5*H )  +  2.5 
IF ( XKW.LT .1.0)  XKW  =1.0 


WSPD2 
VHMK ( I ) 
VHMF ( I ) 
VHPK(I) 
VHPF(I) 


=  WSPD*XKW 
=  VHMK ( I )  +  WSPD2 
=  VHMK( I ) *FPSPKN 
=  VHPK(I)  +  WSPD2 
=  VHPK ( I ) *FPSPKN 


CONTINUE 


DYNAMIC  PRESSURE 


QHM(I)  =  RHOD2  *VHMF ( I ) *  *  2 
QVM(I)  =  RHOD2  *  WMF  ( I )  *  *  2 
QHP(I)  =  RHOD2*VHPF ( I ) **2 
QVP(I)  =  RHOD2 *WPF  ( I )  * *2 

IF ( IOU6 . EQ . IOU1 ) THEN 

IF (LINES.LT. 12 )GOTO  450 
LINES  =  2 
CALL  INKEY 

IF ( KEY . NE . ' C 1 ) GOTO  999 
WRITE ( IOU6 ,1002) 

ENDIF 

CONTINUE 


REPORT  HORIZONTAL  COMPONENTS  OF  VELOCITY  PROFILE 


WRITE ( IOU6 ,1003)  ZZ(I), VHMF ( I ) , VHMK ( I ) , VHPF ( I )  , 
VHPK ( I ) ,QHM(I) , QHP ( I ) 

FORMAT  (  F8.2,2X, 'H' , 6F1 0 . 3 ) 


REPORT  VERTICAL  COMPONENTS  OF  VELOCITY  PROFILE 


WRITE  ( IOU6 , 100  4  )  WMF(I)  ,WMK(I)  ,WPF(I)  ,WPK(I)  , 
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217  *  QVM ( I )  , QVP ( I ) 

218  1004  FORMAT  (  10X, ‘ V ‘ , 6F10 . 3 ) 

219  C 

220  500  CONTINUE 

221  C 

222  C  - 

223  C  WRITE  OUT  GRAPHICS  FILES  IF  SWITCH  IS  SET  BY  USER 

224  C  - 

225  C 

226  IF ( IGRAPH .EQ . 1 ) THEN 

227  C 

228  C  - 

229  C  OPEN  GRAPHICS  FILE 

230  C  - 

231  C 

232  OPEN ( IOU8 ,FILE=PTSFIL ( 2 ) , STATUS3 ' NEW ' ,ERR=2000) 

233  C 

234  WRITE (IOU8, 89)  COMM( 1 ) , COMM( 2 ) 

235  89  FORMAT (  10X, A50 , / , 10X, A50 , // ) 

236  C 

237  WRITE (IOU8, 80)  XIPOUT 

238  80  FORMAT (  IX, ' TITLE="TWIN  ROTOR,  DAIP  =',F6.2,'  FT,’ 

239  *  '  GW  =  xxxxx  LB,  WAGL  =  xx  FT" * ) 

240  C 

241  C  - 

242  C  PRINT  OUT  MEAN  VELOCITY,  PEAK  VELOCITY,  AND  PEAK 

243  C  DYNAMIC  PRESSURE  PROFILES  VERSUS  PROFILE  HEIGHT  (AGL) 

244  C  - 

245  C 

246  WRITE (IOU8, 88) 

247  88  FORMAT (  IX, 'VARIABLES  =  X,HT') 

248  C 

249  WRITE (IOU8, 81) 

250  81  FORMAT (  IX, 'ZONE  T  =  "MEAN  PROFILE,  KTS" ,  I=XX,  F=POINT ’ ) 

251  DO  82  I  =  1 ,NPTS 

252  WRITE (IOU8, 83)  VHMK ( I ) , Z Z ( I ) 

253  83  FORMAT (  IX, F6 . 1 , IX , F5 . 1 ) 

254  82  CONTINUE 

255  C 

256  WRITE (IOU8, 84) 

257  84  FORMAT (  IX, 'ZONE  T  =  "PEAK  PROFILE,  KTS",  I=XX,  F=POINT ' ) 

258  DO  85  I  =  1 , NPTS 

259  WRITE (IOU8, 83)  VHPK( I ) ,ZZ (I ) 

260  85  CONTINUE 

261  C 

262  WRITE (IOU8, 86) 

263  86  FORMAT (  IX, 'ZONE  T  =  "PEAK  Q,  PSF" ,  I=xx,  F=POINT ' ) 

264  DO  87  I  =  1 , NPTS 

265  WRITE ( IOU8 , 83 )  QHP(I),ZZ(I) 

266  87  CONTINUE 

267  C 

268  C  - 

269  C  CLOSE  GRAPHICS  FILE 

270  C  - 
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271 

272 

273 

274 

275 

276 

277 

278 

279 

280 
281 
282 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 

296 

297 

298 

299 

300 

301 

302 

303 


C 


C 

C 

C 

C 

C 

C 

C 

C 

C 


CLOSE ( IOU8 , STATUS3 *  KEEP ’ ) 

END  IF 

CALL  INKEY 
GOTO  999 

THE  ERROR  LOGIC  ALLOWS  FOR  THE  HANDLING  OF  FILE 
OPEN  ERRORS  BY  RETURNING  THE  USER  TO  A  MENU 


2000  CONTINUE 
C 


C 


CALL  HOMCLS ( 0 ) 


WRITE ( IOU1 , 2001 ) 
2001  FORMAT (  ////,8X, 

1  ’  ***  ERROR 

2  // , 8X, ' 

C 


***  PLEASE  CHOOSE 
TYPE  <R£TURN>  TO 


A  NEW  OUTPUT  FILENAME ’ , 
CONTINUE  ' , $ ) 


READ ( IOU1 , ' ( Al ) * )  TEMCHAR 
C 


KEY  =  ’P' 
C 

999  CONTINUE 
C 


RETURN 

END 
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1  C 

2  C 

3  SUBROUTINE  IREAD ( IOU1 , PROMPT , IVALUE ) 

4  C 

5  C 

6  C  SUBROUTINE  IREAD  PROMPTS  USER  FOR  AN  INTEGER 

7  C  DATA  ENTRY  AND  CHECKS  VALIDITY  OF  ENTRY 

8  c  ********************************************** 

9  C 

10  PARAMETER (LAST=50) 

11  C 

12  CHARACTER* 5 0  PROMPT , SHOWIT 

13  CHARACTER* 15  ENTRY, BLANK 

14  C 

15  DATA  BLANK  /'  '/ 

16  C 

17  c  ********************************************** 

18  C 

19  C  - 

20  C  PROMPT  USER  FOR  INTEGER  ENTRY.  FIND  POSITION 

21  C  OF  LAST  NON-BLANK  CHARACTER  IN  PROMPT, 

22  C  THEN  STORE  RIGHT  JUSTIFIED  IN  SHOWIT 

23  C  - 

24  C 

25  N  =  LAST  +  1 

26  C 

27  10  IF ( N . EQ . 1 ) GOTO  20 

28  C 

29  N  =  N  -  1 

30  C 

31  IF ( PROMPT ( N : N ) . EQ .  ' )GOTO  10 

32  C 

33  20  JS  =  LAST  -  N 

34  C 

35  WRITE ( SHOWIT, ' (50A1) ' )  ('  ' , J=1 , JS ) , ( PROMPT ( I s I ) , 1=1 ,N) 

36  C 

37  c  - 

38  C  NOW  ASK  USER  FOR  DATA  ENTRY 

39  c  - 

40  C 

41  30  WRITE ( IOU1 ,  ' ( / , IX , A , 1 3 ) ‘ )  SHOWIT, IVALUE 

42  C 

43  WRITE ( IOU 1,'(/,8X,A,$)') 

44  1  ‘  ENTER  NEW  VALUE  OR  <RETURN>  TO  LEAVE  AS  IS  ==>  ' 

45  C 

46  READ( IOU1 , ' (A) ' )  ENTRY 

47  C 

48  IF ( ENTRY. EQ. BLANK (RETURN 

49  C 

50  READ ( ENTRY , ’  ( BN , 1 7 )  '  , ERR=  3  0 )  ITEMP 

51  C 

52  IVALUE  =  ITEMP 

53  C 

54  RETURN 
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55 

56  C 

57 


END 
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1  C 

2  C 

3 

4  C 

5  C 

6  C 

7  C 

8  C 

9  C 

10  C 

11  C 

12 

13  C 

14  C 

15  C 

16 

17  C 

18 

19 

20 

21  C 

22 

23  C 

24 

25 

26  C 

27 

28 

29  C 

30 


FUNCT I ON  LEGAL ( C HDOL , I OU 1 , OKL 1ST, NUM ) 

**************************************** 
FUNCTION  LEGAL  DETERMINES  IF  THE  VALUE 
FOR  CHDOL  IS  A  VALID  INPUT.  THIS  VALUE 
IS  CHECKED  AGAINST  THE  LIST  OF  LEGAL 
VALUE  IN  ARRAY  OKLIST(NUM) 
**************************************** 

CHARACTER* 1  CHDOL, OKLIST( NUM) 

**************************************** 

LEGAL  =  0 

DO  10  1*1, NUM 

IF (CHDOL .EQ .OKLIST ( I ) )  RETURN 
10  CONTINUE 

LEGAL  *  1 

WRITE ( IOU1 , ' (/ ,T9,A,A1,A/I 1 )  ’  ***  CHDOL, 

1  ’  IS  NOT  A  VALID  INPUT  ***' 

RETURN 

END 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 


C 

C 

c 

c 

c 

c 

c 

c 


c 


c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


c 


c 


c 


c 


c 


SUBROUTINE  LOCATE ( IROW , ICOL ) 

♦  ★★★★♦★★♦★★★★★★★★♦★★★★★♦★★★★★★★★♦★★♦♦★♦★♦♦★♦♦★♦♦★♦Hr** 

SUBROUTINE  LOCATE  LOCATES  THE  CURSOR  POSITION 
***ir**irirA  x  ***•*★★*****★★★★★★★★★★*★**★★★★**★★*★'*★★•**★*★ 

COMMON  /  UNITS/  IOU1 ,  IOU4  ,  IOU5  ,  IOU6  ,  IOU7  ,  IOU8  ,  IGRA01! 

CHARACTER* 8  CUP 
CHARACTER* 1  ECUP ( 8 ) 

EQUIVALENCE  ( CUP , ECUP ( 1 ) ) 

CHARACTER* 10  FMT 
CHARACTER* 1  EFMT{10) 

EQUIVALENCE  ( FMT , EFMT ( 1 ) ) 

DATA  FMT  /  ' ( "  " , A?  \ ) '  / 

♦it*************************************************** 


ANSI  CONTROL  SEQUENCE:  CUP  =  ESC [' ROW COLUMN ' H 


IRl  =  IROW/ 10 

IR2  =  IROW  -  IRl *  1 0 

IC1  =  ICOL/ 1 0 

IC2  =  ICOL  -  IC1*10 

ECUP ( 1 )  =  CHAR (27) 

ECUP ( 2 )  =  CHAR ( 9 1 ) 

IPOS  =  3 

IF( IRl .GT. 0 )THEN 

ECUP (IPOS)  =  CHAR (IRl  4  48) 

IPOS  =  IPOS  +  1 

ENDIF 

ECUP (IPOS)  =  CHAR ( I R2  4  48) 
IPOS  =  IPOS  4  1 

ECUP (IPOS)  =  CHAR (59) 

IPOS  =  IPOS  4  1 

IF ( IC 1 .GT . 0 )THEN 

ECUP (IPOS)  =  CHAR ( IC 1  4  48) 

IPOS  =  IPOS  4  1 

ENDIF 

ECUP ( IPOS )  =  CHAR ( IC2  4  48) 
IPOS  =  IPOS  4  1 
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55 

ECUP(IPOS)  =  CHAR (72) 

56 

C 

57 

EFMT ( 7 )  =  CHAR( IPOS  +  48) 

58 

C 

59 

WRITE ( IOU1 ,* )  CUP 

60 

C 

61 

RETURN 

62 

END 

63 

C 

64 
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1  c 

2  C 

3  SUBROUTINE  MOMENT ( NPTS , HUMTYP , TOTF , TOTM) 

4  C 

5  c  ***************************************************** 

6  C  SUBROUTINE  MOMENT  CALCULATES  THE  TOTAL  OVERTURNING 

7  C  FORCE  AND  MOMENT  ON  A  MAN  OR  YOUNG  PERSON  AND  PRINTS 

8  C  OUT  THE  RESULTS. 

9  c  ***************************************************** 

10  c 

11  CHARACTER* 1  HUMTYP 

12  C 

13  COMMON  /PERSON/  QP(12),DSET 

14  COMMON  /  UNITS/  IOU1 , IOU4 , IOU5 , IOU6 , IOU7 , IOU8 , IGRAPH 

15  C 

16  q  ***************************************************** 

17  C 

18  C  - 

19  C  IN  THE  SUBROUTINE  TO  CALCULATE  FORCES  AND  MOMENTS: 

20  C 

21  C  HCDP"  IS  THE  COEFFICIENT  OF  DRAG  OF  THE  PERSON 

22  C  "WIDTHP"  IS  THE  WIDTH  OF  THE  PERSON  WHERE: 

23  C  WIDTHP  IS  'L'  TYPE  IF  1.1  FT 

24  C  WIDTHP  IS  'S'  TYPE  IF  0.7  FT 

25  C  - 

26  C 

27  WIDTHP  =1.1 

28  CDP  =1.1 

29  C 

30  C  - 

31  C  INITIALIZE  INTEGRATION  STEP  SIZE  AND  ZERO  SUMMATION  VARIABLES 

32  C  - 

33  C 

34  DELZZ  =  0.5 

35  TOTM  =0.0 

36  TOTF  =0.0 

37  C 

38  C  - 

39  C  CHOOSE  HUMAN  SIZE  TYPE 

40  C  - 

41  C 

42  IF ( HUMTYP . EQ .’S') THEN 

43  WIDTHP  =0.7 

44  NPTS  =  8 

45  END  IF 

46  C 

47  C  - 

48  C  INTEGRATE  DYNAMIC  PRESSURE  PROFILE  OVER 

49  C  THE  HEIGHT  OF  THE  PERSON  CHOSEN 

50  C  - 

51  C 

52  DO  10  J  =  1 , NPTS 

53  C 

54  FOVER  =  QP( J) *DELZZ*WIDTHP*CDP 
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ZZ  =  0 . 5* ( J  -  1)  +  0.25 
OVERM  =  FOVER* ZZ 
TOTF  =  TOTF  +  FOVER 
TOTM  =  TOTM  +  OVERM 


PRINT  OUT  RESULTS 


IF ( DSET . NE . 0 . ) GOTO  10 

WRITE (IOU6, 20)  ZZ , OP ( J) , FOVER, OVERM, TOTF, TOTM 
20  FORMAT (  F8 . 2 , 5 ( 2X, F10 . 3 ) ) 

10  CONTINUE 

RETURN 

END 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 


C 

C 

SUBROUTINE  PROPRM ( H , UMB , RVZ ) 

C 

£  ★  *★*★★★****★**★★★*★*★★★★*****★**★*★★★★★★★★★**••*•**★★**★★**★★★★★★★'*★★ 

C  SUBROUTINE  PROPRM 

C 

C  THIS  SUBROUTINE  CALCULATES  THE  VELOCITY  PROFILE  V(R,Z)  PARAMETERS 

C  OF  THE  RADIAL  WALL  JET  FOR  THE  NON-INTERACTING  ROTOR  CASE 

Q  ★★★★★★★★★★★★★*★★★★★*★*★★★★★★★★★★*★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★ 

C 

COMMON  /PROFIL/  RJ , ZB J , ZH J , ZMJ , UMJ, ZB , ZH , ZM, UM, CU , CY 
C 

c 

IF ( RVZ . GE . RJ ) GOTO  600 
C 

c  - 

C  RVZ  .LT.  RJ  ==>  TRANSITION  REGION 

C  - 

c 

UM  =  UMJ*  RVZ 
IF  ( RVZ  .  GT  .1.0)  UM  =  UMJ 
C 

ZH  =  ZHJ*RVZ 
IF ( RVZ .GT .1.0)  ZH  =  ZHJ 
C 

ZM  =  ZMJ*RVZ 
IF ( RVZ . GT .1.0)  ZM  =  ZMJ 
C 

c  - 

C  BOUNDARY  GROWTH  IN  TRANSITION  REGION 

C  - 

C 

ZB0  =1.5 

IF(H.LT. 1.5)  ZB0  =  H 
C 

ZH0  =  ZBO/2.8 

ZH  =  (ZH0  -  ZHJ) /RJ**2* (RJ  -  RVZ)**2  +  ZHJ 
ZB  =  2 . 8*ZH 
C 

GOTO  700 
C 

c  - 

C  RVZ  .GE.  RJ  ==>  DEVELOPED  WALL  JET  REGION 

C  - 

c 

600  CONTINUE 
C 

UM  =  CU*RVZ  **(-1.143) *UMB 
ZH  =  CY*RVZ** (1.028) 

ZB  =  2 . 8*ZH 
ZM  =  0 . 1944*ZH 
C 
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55 

700  CONTINUE 

56 

C 

57 

RETURN 

58 

END 

59 

C 
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:  c 

2  C 

3 

4  C 

5  C 

6  C 

7  C 

8  C 

9  C 

10  C 

11  C 

12  C 

13  C 

14  C 

15  C 

16  C 

17  C 

18  C 

19  C 

20 
21 

22  C 

23  C 

24  C 

25 

26 

27 

28  C 

29 

30 

31 

32  C 

33 

34  C 

35  C 

36  C 

37  C 

38  C 

39 

40 

41  C 

42  C 

43  C 

44  C 

45  C 

46 

47 

48 

49  C 

50  C 

51  C 

52  C 

53  C 

54 


subroutine  vline 

************************  ■-************************************ 
SUBROUTINE  VLINE 

THIS  SUBROUTINE  APPLIES  THE  BIOT-SAVORT  LAW  TO 
CALCULATE  THE  VELOCITY  INDUCED  BY  A  LINE  VORTEX 

XA,  YA, ZA  =  STARTING  POINT  OF  VORTEX 

XB ,  YB , Z B  =  ENDING  POINT,  OR  DIRECTION  POINTER 

XC ,  YC , ZC  =  TARGET  POINT  WHERE  VELOCITY  IS  INDUCED 

IFI  =  0  VORTEX  IS  FINITE,  FROM  POINT  A  TO  POINT  B 

IFI  =  1  VORTEX  IS  SEMI-INFINITE  FROM  POINT  A  THROUGH  B 

*************************************************************** 


COMMON  /CONSTS/  PI , RHO, FPSPKN , RHOD2 , DRC 

COMMON  /CVLINE /  IFI , XA , YA , Z A , XB , YB , ZB , XC , YC , ZC , Q1 , Q2 , Q3 


A  =  ( XA-XC ) *  *2  +  ( YA- YC ) *  *  2  +  (ZA-ZC)**2 

B  =  2 . 0* (  (XA-XB)  +  (XC-XA)  +  ( YA-YB ) * ( YC-YA,  +  ( ZA-ZB ) * ( ZC-SA)  ) 
C  =  ( XA-XB 1**2  +  ( YA-YB ) **2  +  (ZA-ZB)**2 

Cl  =  ( YC-YB ) *ZA  f  ( YA-YC ) *ZB  +  (Yb-YA)*ZC 

C2  =  ( ZC-ZB ) *XA  +  ( ZA-ZC ) *XB  +  (ZB-ZA)*XC 

C3  =  ( XC-XB ) * YA  +  ( XA-XC )*YB  +  (XB-XA)*YC 

Q  =  4 . 0 *A*C  -  B* * 2 


CHECK  FOR  COLINEAR  TARGET  POINT 


QB  =  0.0 

IF ( ABS (Q)  .LT.1.0E-06  JGOTO  100 


FINITE  LENGTH  VORTEX 


IF ( IFI . EQ . 0 ) THEN 

QB  =  1 . 0/Q* ( ( 2 . 0*C  +  B ) / SQRT ( A  +  B  +  C)  -  B/SQRT ( A) ) /2 . 0 /PI 
END  IF 


SEMI-INFINITE  VORTEX 


IF(IFI.EQ. 1 (THEN 
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55  QB  =  1 . 0/Q* ( 2 . 0  *SQRT ( C )  -  B/SQRT( A) ) /2 . 0/PI 

56  ENDIF 

57  C 

58  100  CONTINUE 

59  C 

60  C  - 

61  C  VELOCITY  COMPONENTS 

62  C  - 

63  C 

64  Ql  =  C1*QB 

65  Q2  =  C2*QB 

66  Q3  =  C3*QB 

67  C 

68  RETURN 

69  END 

70  C 
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1  c 

2  C 

3  SUBROUTINE  WAL JET ( H , UB , UN , UMB ) 

4  C 

5  £  ****★★**★*★***★★**★★*★★★★*★★*★★**★*★*★★★★★★★*★★★★★★**★★★*★**■*'**** 

6  C  SUBROUTINE  WALJET 

7  C 

8  C  THIS  SUBROUTINE  DEFINES  THE  START  OF  THE  WALL  JET  AND  ITS  GROWTH 

9  £  *★************★**★**★★★★★★★**★★★*★★★★★★★★★★★*★**★★★*★★★★★★★■*■★***★ 

10  C 

11  COMMON  / CLOUDK/  QSMAX 

12  COMMON  /PROFIL/  RJ , Z BJ , ZHJ , ZMJ , UMJ , ZB , ZH , ZM, UM, CU , CY 

13  COMMON  /  UNITS/  IOU1 , IOU4 , IOU5 , IOU6 , IOU7 , IOU8 , IGRAPH 

14  C 

16  C 

17  C  - 

18  C  INITIALIZATION  OF  EXPONENTS 

19  c  - 

20  C 

21  EXU  =  -1.143 

22  EXY  =  1.028 

23  EXM  =  1.0  +  EXU  +  EXY 

24  CQ  =  0.52 

25  CZM  =  0.1944 

26  CZB  =  2.8 

27  C 

28  C  - 

29  C  ITERATE  TO  FIND  INITIAL  RADIUS  OF  WALL  JET,  RJ 

30  C  - 

31  C 

32  TOL  =  1 . 0E-05 

3  3  R J  =  2 . 0 

34  C 

35  QSMAX0  =1.0 

36  C 

37  DO  100  1=1,20 

38  C 

39  c  - 

40  C  EQUIVALENT  JET  LENGTH  T/R 

41  c  - 

42  C 

43  TR  =  H  +  ( R J  -  1.0) 

44  TDE  =  0 . 7  07  *TR 

45  C 

46  c  - 

47  C  QSMAX  CURVEFIT  TO  FIG.  8,  USAAVLABS  TECHNICAL 

48  C  REPORT  68-52,  JULY  1968 

49  C  - 

50  C 

51  IF( TDE . LE .4.0)  QSMAX  =  QSMAX 0  +  (0.6  -  QSMAX 0 ) / 1 6 . 0*TDE* *2 

52  IF ( TDE . GT .4.0)  QSMAX  =  2 . 4 /TDE 

53  UM  =  SQRT( QSMAX) 

54  C 
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55  RJNEW  =  2.508078*(UB/UM)**(0.486) 

56  C 

57  IF (ABS (RJNEW  -  RJ ) . LE . TOL ) GOTO  200 

58  RJ  =  RJNEW 

59  C 

60  100  CONTINUE 

61  C 

62  WRITE (IOU1, 10) 

63  10  FORMAT (  •***********************************************',/ 

64  1  ,  1  ITERATIONS  EXCEEDED  FOR  WALL  JET  INITIAL  RADIUS’,/ 

65  2  ,•****************★******************************•) 

66  C 

67  STOP  ’  ’ 

68  C 

69  200  CONTINUE 

70  C 

71  RJ  =  RJNEW 

72  C 

73  C  - 

74  C  VELOCITY  GROWTH  FUNCTION  CONSTANTS 

75  C  - 

76  C 

77  UMB  =  ( (0.3586*RJ**EXM*(UM*UN)*(UB*UN)**(0.14) )**(0.88) )/UN 

78  ZHJ  =  0.654/ (UM/UMB) **2/RJ 

79  CU  =  UM/UMB^RJ*  *  ( -EXU ) 

80  CY  =  ZHJ*RJ**(-EXY) 

81  C 

82  C  - 

83  C  MAX  VELOCITY  AND  BOUNDARY  PARAMETERS  AT  RJ 

84  C  - 

85  C 

86  UMJ  =  CU*RJ** (EXU) *UMB 

87  ZHJ  =  CY*RJ*  * ( EXY ) 

88  ZMJ  =  CZM*  ZHJ 

89  ZBJ  =  CZB*ZHJ 

90  C 

9 1  RETURN 

92  END 

93  C 

94 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 


C 

C 

SUBROUTINE  WJVEL ( H , UN , UMB , RVZ , RADIUS , WSPD , DELZ , ZMAX , DXO ) 
C 

£  ★  ★**★*★**★*★****★★***★*★★***★★★*★★★★★*’★★★★★*★★★★★★★★★★★ 

C  SUBROUTINE  WJVEL  GENERATES  THE  VELOCITY  PROFILE  V(R,Z) 

C  AT  RVZ  FOR  THE  NON-INTERACTING  ROTOR  CASE 

C  ********************************************************* 

C 

CHARACTER* 1  TEMCHAR 
CHARACTER* 1  ICONT ( 5 ) 

CHARACTER* 1  KEY, KKEY 
CHARACTER* 12  PTSFIL(4) 

CHARACTER* 50  COMM(2) 

C 

DIMENSION  ZZ (40) ,VMF(40) ,VMK(40) ,VPF(40) , 

1  VPK (40), QM (40), QP (40) 

C 

COMMON  /  CKEY /  KEY, KKEY 

COMMON  /CONSTS/  PI , RHO , FPSPKN , RHOD2 , DRC 

COMMON  /INPUTC/  ICONT , COMM , PTSFIL 

COMMON  /PROFIL/  RJ, ZBJ, ZH J, ZMJ , UMJ , ZB, ZH , ZM, UM, CU, CY 
COMMON  /  UNITS/  IOU1 , IOU4 , IOU5 , IOU6 , IOU7 , IOU8 , IGRAPH 
C 

Q  **★★**********★****★***★**★★*★★★*★★★★★★*★*★*•**★★★★★★*★★ 

c 

c  - 

C  ' PROPRM '  PROVIDES  THE  VELOCITY  PROFILE  PARAMETERS 

C  OF  A  RADIAL  WALL  JET  (WITHOUT  INTERACTION  PLANE) 

C  - 

C 

CALL  PROPRM (H, UMB, RVZ) 

C 

c  - 

C  DIMENS IONALIZE  VELOCITY  PROFILE  PARAMETERS 

C  - 

c 

RRVZ  =  RADIUS*RVZ 
ZZB  =  ZB*RADIUS 
ZZH  =  ZH*RADIUS 
ZZM  =  ZM*RADIUS 
ZETAM  =  ZM/ZB 
C 

C  - 

C  OUTPUT  THE  VELOCITY  AND  DYNAMIC  PRESSURE  PROFILE  HEADER 

C  - 

C 

ICD  =  0 

CALL  HOMCLS ( ICD ) 

IF ( IOU6 . NE . IOU1 >  WRITE ( IOU6 ,  '  (  '  ' 1 ’  '  )  '  ) 

C 

IF ( IOU6 . EQ . 6 )  WRITE ( IOU6 ,93)  COMM( 1 ) ,COMM( 2 ) 

93  FORMAT (  1  OX , A50 , / , 1  OX , A50 , // ) 

C 
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55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 


RVZOUT  =  RRVZ  +  DXO 
C 

WRITE (IOU6, 1000)  RVZOUT 

1000  FORMAT (  9X, ' SINGLE  ROTOR  VELOCITY  PROFILE  AT  RADIUS  = 

1  F7 . 1 , ’  FT',/) 

C 

WRITER IOU6 , 1001 )  ZZB , ZZH , ZZM 

1001  FORMAT (  15X,  1  PROFILE  BOUNDARY  HEIGHT  =  ',F7.2,'  FT',/ 

1  , 15X, '  HALF-VEL. HEIGHT  =  *,F7.2,'  FT',/ 

2  , 15X,  MAX-VEL  HEIGHT  =  ',F7.2,'  FT',/) 

C 

C  - 

C  INCREMENTS  AND  HEIGHT  ARE  FROM  DELZ  AND  ZMAX 

C  - 

c 

NPTS  =  IFIX ( ZMAX /DELZ )  +  1 
C 

c  - 

C  BOUNDARY  LAYER  REGION  EXPONENT 

C  'AN'  IS  ACTUALLY  =  1. 0/7.0 

C  - 

c 

AN  =  0.142857142 
C 

c  - 

C  SHEAR  LAYER  REGION  EXPONENT,  TO  MEET  EDGE  CONDITIONS 

C  (FROM  FIGURE  7,  USAAVLABS  TECHNICS  REPORT  68-52,  JULY  1968) 

C  - 

C 

ALPW  =  ALOG (1.0  -  1 .0/SQRT(2 .0) )/ALOG( (ZH  -  ZM)/(ZB  -  ZM) ) 

C 

VN  =  UN 
VMM  =  UM 
C 

c  - 

c  CALCULATION  OF  THE  PEAK  VELOCITY  CORRECTION  FACTOR 

C  IS  MADE  IN  THE  NEXT  SECTION  AT  THE  3  FT  POSITION 

C  - - 

C 

VZ3  =  3.0/RADIUS 
ZETA3  =  VZ  3/ZB 
VZM3  =  0.0 
C 

IF ( ZETA3 .GE . ZETAM) GOTO  51 
C 

c  - 

C  THE  3  FT  HEIGHT  IS  IN  THE  BOUNDARY  LAYER 

C  - 

C 

IF ( ZETAM . GT .0.0)  VZM3  =  ( ZETA3 /ZETAM )* *AN 
C 

GOTO  52 
51  CONTINUE 


C 
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109  C  - 

110  C  THE  3  FT  HEIGHT  IS  IN  THE  SHEAR  LAYER 

111  C  - 

112  C 

113  IF ( V23 .LE . ZB ) THEN 

114  VZM3  =  (1.0  -  ( ( ZETA3  -  ZETAM)/(1.0  -  ZETAM) ) **ALPW) **2 

115  ENDIF 

116  C 

117  52  CONTINUE 

118  C 

119  C  - 

L20  C  THE  VELOCITY  TO  ADD  TO  THE  MEAN  VELOCITY  TO  OBTAIN  THE  PEAK 

L21  C  VELOCITY  IS  CALCULATED  AS  THE  CONSTANT  VMFD3  (FT/SEC) 

122  C  - 

123  C 

124  VZN3  =  VZM3*VMN 

125  VMF3  =  VZN3*VN 

126  VMFD3S  =  (RVZ*0.2444)  +  0.4 

127  IF ( VMFD3S . GT .1.5)  VMFD3S  =1.5 

128  VMFD3  =  VMFD3S* VMF  3 

129  C 

130  C  - 

131  C  OUTPUT  THE  VELOCITY  AND  DYNAMIC  PRESSURE  PROFILE  HEADER 

132  C  - 

133  C 

134  WRITE ( IOU6 ,1005) 

135  1005  FORMAT (  2X ,' HEIGHT ', 5X, ' MEAN  VELOCITY  1 , 7X PEAK  VELOCITY ', 6X, 

136  1  'MEAN  Q' ,4X, 'PEAK  Q' ,/, 

137  2  3Xf ' (FT) ’ ,5X, ' (FPS) ' ,6X, ’ (KN) ' ,5X, ' (FPS) ’ ,6X, ' (KN) ’ ,5X, 

138  3  ’ (PSF) ' ,5X, ' (PSF) ’ ,/) 

139  C 

140  C  - 

141  C  CALCULATE  THE  VELOCITY  PROFILE  POINTS  FOR  OUTPUT 

142  C  - 

143  C 

144  LINES  =  0 

145  C 

146  DO  500  I  =  1 , NPTS 

147  C 

148  LINES  =  LINES  +  1 

149  Z  =  DELZ*FLOAT( I  -  1) 

150  ZETA  =  Z/ZB 

151  IF ( ZETA. GE. ZETAM) GOTO  300 

152  C 

153  C  - 

154  C  Z  IS  WITHIN  BOUNDARY  LAYER 

155  C  - 

156  C 

157  VZM  =  0.0 

158  IF ( ZETAM . GT .0.0)  VZM  =  ( ZETA/ZETAM ) * *AN 

159  GOTO  400 

160  C 

161  300  CONTINUE 

162  C 
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163  C 

164  C 

165  C 

166  C 

167 

168  C 

169 

170 

171 

172  C 

173  400 

174  C 

175 

176  C 

177  C 

178  C 

179  C 

180  C 

181 

182  C 

183  C 

184  C 

185  C 

186  C 

187 

188 

189  C 

190  C 

191  C 

192  C 

193  C 

194 

195 

196 

197  C 

198 

199  C 

200  C 

201  C 

202  C 

203  C 

204  C 

205  C 

206 

207 

208 

209 

210 
211 
212 

213  C 

214  55 

215  C 

216  C 


2  IS  WITHIN  SHEAR  LAYER 

VZM  =0.0 

IF(  Z .LE . ZB ) THEN 

VZM  =  (1.0  -  ((ZETA  -  ZETAM ) / ( 1 . 0  -  ZETAM) ) “ALPW) “2 
ENDIF 

CONTINUE 

VZN  =  VZM*VMN 

DIMENSIONAL  HEIGHT 

ZZ(I)  =  Z  *RADIUS 

MEAN  VELOCITIES 


VMF(I)  =  VZN*VN 
VMK(I)  =  VMF ( I ) /FPSPKN 


PEAK  VELOCITIES 


VPF(I)  =  VMF ( I )  +  VMFD3 

IF ( VMF ( I ) . EC . 0 . )  VPF(I)  =  0.0 

VPK(I)  =  VPF( I } /FPSPKN 

IF ( VPK ( I ) • EQ .0.0) GOTO  55 


THE  EFFECT  OF  WIND  IS  TO  ADD  (DOWNWIND  SIDE)  OR  SUBTRACT 
(UPWIND  SIDE)  ' XWK '  TIMES  THE  AMBIENT  WIND  VELOCITY  TO 
THE  HORIZONTAL  PROFILE  VELOCITY  (EMPIRICAL,  CH-53E  BASED) 


XKW  =  ( -0 . 5*H )  +  2.5 
IF ( XKW . LT .1.0) XKW  =  1.0 
WSPD2  =  WSPD*XKW 
VMK(I)  =  VMK(I)  +  WSPD2 
VMF ( I )  =  VMK( I ) ‘FPSPKN 
VPK ( I )  =  VPK ( I )  +  WSPD2 
VPF(I)  =  VPK( I) ‘FPSPKN 

CONTINUE 
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217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 


C  DYNAMIC  PRESSURE 

C  - 

c 

QM(I)  =  RHOD2*VMF(I)**2 
QP{I)  =  RHOD2*VPF (I) **2 
C 

IF  ( IOU6 . EQ . IOU1 ) THEN 
IF (LINES .LT . 12 )GOTO  450 
LINES  =  1 
CALL  INKEY 

IF(XEY.NE. 'C ' )GOTO  999 
WRITE (IOU6, 1005) 

END  IF 
C 

450  CONTINUE 
C 

c  - 

c  OUTPUT  THE  VELOCITY  AND  DYNAMIC  PRESSURE  PROFILES 

C  - 

C 

WRITE (IOU6, 1002)  ZZ(I) ,VMF(I) ,VMK(I) ,VPF(I) , 

*  VPK(I) ,QM(I) ,QP(I) 

1002  FORMAT (  F8 . 2 , 6F1 0 . 3 ) 

C 

500  CONTINUE 


C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 

c 

c 

c 

c 


c 


WRITE  OUT  GRAPHICS  FILES  IF  SWITCH  IS  SET  BY  USER 


IF ( I GRAPH . EQ . 1 ) THEN 


OPEN  GRAPHICS  FILE 


OPEN ( IOU8 , FILE=PTSFIL ( 1 ) , STATUS= ' NEW ' ,ERR=2000) 

WRITE (IOU8, 89)  COMM( 1 ) , COMM ( 2 ) 

89  FORMAT (  1  OX, A50 , / , 10X , A50 , // ) 

WRITE (IOU8, 80)  RVZOUT 

80  FORMAT (  IX, ' TITLE=" SINGLE  ROTOR,  DFRC  =',F6.2,'  FT,' 

*  '  GW  =  xxxxx  LB,  WAGL  =  XX  FT"’) 


PRINT  OUT  MEAN  VELOCITY,  PEAK  VELOCITY,  AND  PEAK 
DYNAMIC  PRESSURE  PROFILES  VERSUS  PROFILE  HEIGHT  (AGL) 


WRITE (IOU8, 88) 

88  FORMAT (  IX, 'VARIABLES  =  X,HT') 
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271 

WRITE (IOU8, 81) 

272 

81 

FORMAT (  IX , ' ZONE  T  =  "MEAN  PROFILE,  KTS " ,  I=xx,  F=POINT ' ) 

273 

C 

274 

DO  82  I  =  1 , NPTS 

275 

WRITE (IOU8, 83)  VMK(I),Z2(I) 

276 

83 

FORMAT (  IX , F6 . 1 , 1 X , F5 . 1 ) 

277 

82 

CONTINUE 

278 

C 

279 

WRITE (IOU8, 84) 

280 

84 

FORMAT (  IX, 'ZONE  T  =  "PEAK  PROFILE,  KTS",  I=xx,  F=POINT ’ ) 

281 

C 

282 

DO  85  I  =  1 , NPTS 

283 

WRITE (IOU8, 83)  VPK(I),ZZ(I) 

284 

85 

CONTINUE 

285 

C 

286 

WRITE (IOU8, 86) 

287 

86 

FORMAT (  IX, 'ZONE  T  =  "PEAK  Q,  PSF",  I=xx,  F=POINT ' ) 

288 

C 

289 

DO  87  I  =  1 , NPTS 

290 

WRITE ( IOU8 ,83)  QP(I),ZZ(I) 

291 

87 

CONTINUE 

292 

C 

293 

C 

294 

C 

CLOSE  GRAPHICS  FILE 

295 

C 

296 

C 

297 

CLOSE ( IOU8 , STATUS3 ' KEEP ' ) 

298 

C 

299 

ENDIF 

300 

C 

301 

CALL  INKEY 

302 

C 

303 

GOTO  999 

304 

C 

305 

c 

306 

c 

THE  ERROR  LOGIC  ALLOWS  FOR  THE  HANDLING  OF  FILE 

307 

c 

OPEN  ERRORS  BY  RETURNING  THE  USER  TO  A  MENU 

308 

c 

309 

c 

310 

2000 

CONTINUE 

311 

c 

312 

CALL  HOMCLS ( 0 ) 

313 

WRITE ( IOU1 ,2001 ) 

314 

2001 

FORMAT (  ////, 8X, 

315 

1 

***  ERROR  ***  PLEASE  CHOOSE  A  NEW  OUTPUT  FILENAME' 

316 

2 

!  / / , 8X , '  TYPE  <RETURN>  TO  CONTINUE  ’,$) 

317 

READ ( I OU 1 , ' ( A 1 ) ' )  TEMCHAR 

318 

KEY  =  'P' 

319 

c 

320 

999 

CONTINUE 

321 

c 

322 

RETURN 

323 

END 

324 

c 
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